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Abstract 


The  problem  of  flow  separation  from  a  low  pressure  turbine  blade  was  investigated.  The  op¬ 
erating  conditions  under  which  the  separation  occurred  were  documented  through  measurement  of 
surface  pressure  coefficients,  boundary  layer  velocity  and  turbulence  profiles,  total  pressure  loss 
coefficient  and  wake  velocity  momentum  deficit.  Three  different  means  for  reducing  the  losses 
associated  with  the  flow  separation  were  also  investigated.  A  boundary  layer  trip,  dimples,  and 
V-grooves  were  studied  as  passive  means  requiring  no  additional  energy  to  reduce  the  separation 
losses.  The  boundary  layer  trip  was  only  successful  for  an  inlet  and  axial  chord  Reynolds  number  of 
50k  with  a  reduction  in  loss  coefficient  of  58.2%.  Three  sets  of  dimples  were  tested  with  the  place¬ 
ment  of  each  at  axial  chord  locations  of  50%,  55%,  and  65%.  The  dimples  provided  reductions  in 
the  loss  coefficient  for  Reynolds  numbers  of  50k,  100k,  and  200k  ranging  from  5.1%  (Re  =  100k, 
freestream  turbulence  level  of  4%)  to  51.7%  (Re  =  50k,  freestream  turbulence  level  of  4%).  Two 
sets  of  V-grooves  were  tested  with  axial  chord  start  locations  of  55%  and  60%.  The  V-grooves  pro¬ 
vided  smaller  reductions  in  loss  coefficient  than  the  dimples .  Boundary  layer  profiles,  total  pressure 
loss  coefficients,  and  wake  velocity  momentum  deficits  are  presented  for  the  three  passive  modifi¬ 
cations. 
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Flow  Separation  Prevention  on  a  Turbine  Blade  in 
Cascade  at  Low  Reynolds  Number 


Chapter  1  -  Introduction 

The  recent  Air  Force  interest  in  the  use  of  Unmanned  Aerial  \fehicles  (UAV)  for  reconnaissance 
has  led  to  new  requirements  for  turbomachinery  engine  designers.  The  low  Reynolds  number  effects 
of  high  altitude,  low  temperatures  and  sustained  low  speed  cruise  velocities  for  smaller  engines  has 
shown  inefficiencies  due  to  flow  separation  losses  with  the  low  pressure  turbine  (LPT).  For  example, 
UjW  aircraft  such  as  the  Teledyne  Ryan  Tier  11+  Global  Hawk  using  the  Allison  AE3007-H  3-stage 
low  pressure  turbine  (LPT)  suffer  from  this  loss  of  efficiency  (MacArthur  [1])  while  flying  at  its 
operational  altitude  of  19.812  km  (65,000  feet)  at  a  standard  cmise  velocity  of  176.4  m/s  (343  Kts 
TAS).  The  author  proposes  that  surface  modifications  to  LPT  blades  can  reduce  these  losses  by 
maintaining  attached  flow  across  the  surface  of  the  blades.  With  the  blade  efficiency  increased, 
loiter  time,  range,  and  performance  will  all  improve. 

To  solve  the  loss  in  efficiency  problem  a  technique  has  been  developed  that  promotes  attached 
flow  over  the  suction  surface  of  the  blade  at  low  Reynolds  numbers  but  does  not  inhibit  the  already 
acceptable  flow  field  present  with  high  Reynolds  numbers.  To  accomplish  this,  passive  techniques 
requiring  no  additional  power  from  the  engine  or  additional  weight  have  been  studied.  These  passive 
techniques  provide  a  solution  to  one  of  the  most  recently  discovered  problems  with  turbomachinery. 

1.1  Low  Pressure  Turbine 

Investigation  of  LPT  efficiency  has  remained  largely  neglected  until  recent  years  despite  the 
many  improvements  and  research  conducted  over  the  last  six  decades  for  other  turbomachinery 
components.  The  compressor  has  been  improved  by  increasing  its  efficiency  and  overall  pressure 
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ratio  as  well  as  the  individual  stage  pressure  ratios.  The  high  pressure  turbine  has  been  improved  by 
increasing  its  allowable  maximum  temperature  through  cooling  techniques.  The  lower  temperature 
characteristics  of  the  low  pressure  turbine  do  not  require  active  cooling  in  modem  engines,  and 
hence  more  focus  is  given  to  improving  the  LPT  blades’  aerodynamic  performance. 

Pratt  and  Whitney  introduced  the  Pak-B  blade  in  an  attempt  to  improve  the  low  Reynolds 
number  efficiency  of  the  LPT.  The  Pak-B  blade  shape  is  a  modified  version  of  the  commercially  used 
Pak  airfoil.  The  Pak-B  airfoil  was  intentionally  designed  for  lower  Mach  number  operation.  While 
in  operation,  each  of  the  blades  in  the  low  pressure  turbine  shed  wakes  which  travel  downstream. 
When  these  wakes  come  into  contact  with  a  blade  in  the  next  row,  the  wake  passes  over  the  surface. 
This  region  inside  the  wake  has  a  higher  ffeestream  turbulence  level  than  the  air  upstream  and 
downstream  of  it.  The  turbulence  level  inside  the  wake  may  reach  up  to  20%  (Halstead  [2,3]).  While 
this  higher  turbulence  level  provides  a  beneficial  effect  with  boundary  layer  attachment,  the  duration 
is  only  intermittent.  Between  wakes,  the  flow  over  the  suction  surface  of  the  blade  is  laminar 
with  separation  occurring  for  the  low  Reynolds  numbers  typical  of  LbW  missions.  Experiments 
by  Murawski  et  al.  [4]  and  Qiu  and  Simon  [5]  have  shown  that  losses  associated  with  separation 
decrease  the  effectiveness  of  the  Pak-B  blades  (see  Figure  1)  at  low  Reynolds  numbers. 


Figure  1.  Unsteady  wake  passing  over  Pak-B  blade  surface  (mounting  pin  locations  also  shown). 
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In  addition  to  the  work  of  Murawski  et  al.  [4],  ongoing  efforts  by  Simon  and  Volino  [6]  fo¬ 
cusing  on  the  low  pressure  Pak-B  turbine  blade  have  also  shown  separation  losses  for  applications 
at  high  altitudes  and  low  temperatures.  These  lower  Reynolds  number  effects  have  been  receiving 
increasing  attention  and  the  problem  is  now  understood  as  one  of  separation  and  transition  over  the 
suction  surface  of  the  blade  due  to  regions  of  strong  pressure  gradient  and  diffusion.  Hourmouziadis 

[7]  has  shown  that  for  these  low  Reynolds  number  environments  profile  losses  rise  sharply  and  ef¬ 
ficiency  is  decreased. 

1.2  Previous  Low  Pressure  Turbine  Experiments 

To  document  the  existing  baseline  efficiency  of  the  Pak-B  blade,  linear  cascade  tunnels  have 
been  used  to  simulate  the  engine  operating  environment.  A  linear  cascade  tunnel  contains  blades 
arranged  in  a  row,  as  if  the  low  pressure  blade  spool  were  cut  and  unrolled  onto  a  flat  surface.  Many 
experimenters  have  chosen  to  use  linear  cascades  for  their  LPT  experiments  including  Rivir  [8], 
Murawski  et  al.  [4],  and  Qiu  and  Simon  [9],  Murawski  et  al.  [4]  and  Qiu  and  Simon  [9]  have  docu¬ 
mented  the  existing  separation  that  exists  on  the  Pak-B  LPT  blade  at  low  Reynolds  numbers.  Rivir 

[8]  investigated  the  effect  of  turbulent  length  scale  on  transition  location  for  the  Langston  airfoil 
previously  described  by  Baughn  et  al.  [10],  Murawski  et  al.  [4]  and  Qiu  and  Simon  [9]  verified  the 
basic  low  Reynolds  number  separation  problem  of  the  Pak-B  airfoil  with  measurements  of  velocity, 
surface  pressure,  separation  location,  and  boundary  layer  thickness.  Murawski  et  al.  [4]  found  the 
Pak-B  to  have  suction  surface  separation  extending  from  the  75%  to  90%  chord  at  a  low  inlet  and 
axial  chord  Reynolds  number  of  53,000.  Qiu  and  Simon  [9]  documented  the  reattachment  of  turbu¬ 
lent  flow  downstream  of  the  90%  chord  position.  The  boundary  layer  thickness  shows  a  significant 
increase  with  flow  separation,  which  decreases  overall  performance.  To  better  understand  the  effect 
of  separation  bubbles,  a  brief  explanation  is  given. 
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Werle  [11]  classified  separation  that  occurs  on  turbine  blades  into  four  categories  that  include 
leading  edge  bubbles,  transonic  shock  induced  separation,  pressure  side  bubbles,  and  trailing  edge 
wake  separation.  Recent  tests  by  Murawski  et  al.  [4]  and  Qiu  and  Simon  [9]  have  shown  that  oper¬ 
ation  at  low  Reynolds  numbers  can  lead  to  an  additional  form  of  separation  where  the  laminar  flow 
detaches  due  to  the  strong  streamwise  pressure  gradient.  Once  the  laminar  flow  separates  from  the 
suction  side  of  the  blade,  transition  can  occur.  The  separated  flow  transition  which  occurs  on  the 
Pak-B  blade  is  one  of  three  modes  of  transition  that  have  been  defined  by  Mayle  [12]  and  Halstead 
[13].  These  three  modes  of  boundary  layer  transition  are  natural  transition,  bypass  transition  and 
separated-flow  transition.  The  control  of  this  particular  form  of  separated-flow  transition  and  lam¬ 
inar  separation  over  the  Pak-B  blade  is  the  focus  of  this  research. 

Most  past  experiments  have  focused  on  the  airflow  over  blades  of  existing  design,  with  almost 
no  comparable  experiments  on  slightly  modified  blades  to  improve  upon  the  performance.  One  ex¬ 
ception  to  this  is  Chen  etal.  [14]  who  attempted  to  decrease  separation  losses  for  compressor  airfoils 
using  riblets  on  the  suction  side  of  the  airfoil  -  the  experiment  was  unsuccessful.  The  poor  perfor¬ 
mance  of  the  riblets  on  convex  surfaces  excluded  their  use  in  this  experiment.  The  modification  of 
turbine  blades  with  passive  surface  features  to  reduce  separation  has  never  been  attempted  before. 
The  possibilities  for  passive  surface  modifications  goes  beyond  just  low  pressure  turbines.  Passive 
techniques  could  be  used  on  other  aircraft  parts  such  as  propellers  (Occhipinti  [15])  or  even  auto¬ 
mobiles.  This  research  documents  successful  integration  of  surface  modifications  with  the  Pak-B 
airfoil  which  reduced  the  losses  associated  with  separated  flow. 

1.3  Current  Research  Focus 

This  project  focused  on  the  control  of  low  Reynolds  number  separation  effects  on  low  pres¬ 
sure  turbine  blades  of  the  Pratt  and  Whitney  Pak-B  blade  design.  In  particular,  measurements  of 
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separation  and  loss  coefficient  were  made  for  the  unmodified  Pak-B  blades.  Data  was  taken  for 
freestream  turbulence  (Tu)  levels  of  1%  and  4%  to  model  various  engine  operating  conditions  for 
low  inlet  Reynolds  numbers  (Re)  of  50,000,  100,000  and  200,000.  These  Reynolds  numbers  are 
typical  of  the  conditions  a  U/V  low  pressure  turbine  blade  will  operate  under  during  a  mission.  The 
turbulence  intensities  are  typical  of  the  regions  of  flow  over  the  blade  between  wakes. 

The  modifications  to  the  Pak-B  blade  included  the  use  of  a  trip  wire,  dimples,  and  longitudi¬ 
nal  V-grooves  as  a  means  of  reducing  separation  losses.  The  benefits  of  dimples  in  a  cylinder  was 
shown  by  Bearman  and  Harvey  [16, 17],  while  application  of  V-grooves  to  the  surface  of  a  con¬ 
vex  ramp  was  tested  by  Lin  et  al.  [18].  An  investigation  into  the  use  of  a  simple  wire  to  trip  the 
boundary  layer  from  a  laminar  state  to  a  turbulent  one  was  also  investigated.  More  information  on 
the  above  surface  modifications  are  presented  in  Chapters  7,  8,  and  9.  A  typical  test  sequence  in¬ 
cluded  blade  surface  modifications  made  to  an  existing  central  blade  of  the  cascade  and  subsequent 
re-examination  (varied  Reynolds  number  and  turbulence  level)  for  the  improvements  in  reducing 
separation  and  improving  loss  coefficient.  The  overall  purpose  of  this  project  was  to  demonstrate 
that  a  turbine  blade  with  separation  problems  at  low  Reynolds  numbers,  if  properly  modified,  will 
experience  reduced  flow  separation  and  improved  blade  performance. 


1.4  Determination  of  Improved  Performance 

To  better  judge  the  merit  of  the  research  a  numerical  parameter  is  required  for  comparison. 
Improvements  to  aircraft  performance  such  as  range,  fuel  economy  and  engine  thermal  efficiency 
are  factors  which  can  be  easily  compared  to  a  baseline  case.  Aircraft  range  can  be  estimated  from 
the  Breguet  range  formula  (Anderson  [19]): 


Uo£l,Wo 
CD  Wi 


(1) 
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where  rj0  is  the  overall  engine  efficiency  and  c  is  the  specific  fuel  consumption.  Also  in  equation  1 
is  the  ratio  of  lift  (Cl)  to  drag  (Cx>)  coefficients  and  the  ratio  of  the  initial  weight  (Wo)  of  the  aircraft 


to  the  final  weight  (WO.  Equation  1  can  be  re-arranged  in  terms  of  the  weight  of  fuel  expended: 


AWfuel  =  Wi 


exp 


(RcCd\ 

Uo  CL)  . 


(2) 


where  AWfuei  =  Wo  —  Wi.  Examination  of  equation  1  and  equation  2  shows  two  definite  advan¬ 


tages  to  improving  the  overall  engine  efficiency  -  increased  range  or  decreased  fuel  consumption. 
The  thermal  efficiency  is  expressed  by  the  increased  kinetic  eneigy  imparted  to  the  gas  per  unit 


of  heat  eneigy  supplied  by  the  fuel.  For  a  stationary  engine  (e.g.  in  a  test  cell)  this  simplifies  to 

_(*) 


nth 


fhv 


(3) 


where  eg  is  the  nozzle  exit  velocity,  /  is  the  fuel-air  ratio  and  hv  is  the  heating  value  of  the  fuel.  It 
will  be  shown  later  that  reduced  total  pressure  losses  result  in  an  increase  in  eg  and  hence  an  increase 
in  thermal  efficiency. 

The  engine  efficiency  can  be  determined  with  the  aid  of  a  T-S  diagram.  A  T-S  diagram  depicts 
the  change  in  entropy  with  the  change  in  temperature  for  a  given  engine  cycle.  For  example,  the 
ideal  jet  engine  can  be  depicted  with  the  Brayton  cycle  shown  in  Figure  2  (Mattingly  [20]).  Figure 
2a  depicts  the  T-S  diagram  for  the  components  shown  in  Figure  2b.  The  compressor  increases  the 
pressure  of  the  air  from  step  2  to  3  without  an  increase  in  entropy  (isentropically).  As  seen  in  Figure 
2a,  the  increase  in  pressure  also  increases  the  temperature  of  the  air.  From  step  3  to  4  the  fuel  is 
burned  and  heat  is  added  to  the  system.  In  the  ideal  process,  the  addition  of  heat  to  the  system 
is  assumed  to  occur  at  a  constant  pressure.  The  hot  air  expands  from  step  4  to  9  as  the  air  exits 
through  the  turbine  blades  and  exhaust  nozzle.  The  expansion  through  the  turbine  is  also  assumed 
to  be  isentropic.  Finally,  the  heat  from  the  engine  is  rejected  externally  through  a  constant  pressure 
process  from  step  9  to  0. 
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a)  T-S  diagram 


Compressor 


On 


Figure  2.  Brayton  cycle  for  an  ideal  gas  jet  engine. 
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The  work  of  each  of  the  components  in  Figure  2b  are  depicted  as  W.  The  work  of  the  com¬ 
pressor  is  defined  as: 

Wc=m  Cp  (Tt3  -  Tm)  (4) 

where  m  is  the  mass  flow  rate  and  Cp  is  the  ratio  of  specific  heats  at  constant  pressure.  The  change 
in  temperature  across  the  compressor  is  included  as  Ttz  —  Tt 2.  The  work  of  the  turbine  is  equal  to 
the  work  required  by  the  compressor: 

Wc=Wt=m  cp  (r*4  -  Tts)  (5) 

where  Wt  is  defined  similarly  to  equation  4  for  the  turbine.  The  actual  work  output  of  a  turbine  is 
less  than  ideal  due  to  losses.  The  losses  in  an  engine  can  be  described  with  a  T-S  diagram  for  each 
part  of  the  jet  engine. 

Figure  3a  shows  the  station  numbers  for  the  stator  and  rotor  blades.  The  passage  between  the 
blades  is  from  station  2  to  2R.  The  exit  passage  to  the  next  stage  is  from  station  3Rto  3.  The  work 
done  by  the  rotor  blades  occurs  from  station  2R  to  3R.  For  discussion  purposes  it  will  be  assumed 
that  the  flow  in  the  space  between  the  blade  rows  is  isentropic.  The  T-S  diagram  for  a  typical  low 
pressure  turbine  stage  is  illustrated  in  Figure  3b. 

As  seen  in  Figure  3b,  the  generation  of  entropy  over  the  rotor  decreases  the  total  pressure  at 
station  3R.  This  decrease  in  total  pressure  is  carried  through  to  the  exit  (station  3)  of  the  stage.  The 
greater  the  drop  in  stagnation  pressure,  the  lower  the  efficiency  of  the  stage.  The  passive  techniques 
tested  reduce  this  pressure  loss  and  therefore  increase  the  stage  efficiency  of  the  low  pressure  turbine . 

The  results  of  the  experiment  discussed  in  Chapters  6  through  9  showed  the  potential  for  re¬ 
ducing  the  stagnation  pressure  drop  over  a  set  of  blades.  A  reduction  in  pressure  drop  increases  the 
efficiency  of  the  low  pressure  turbine  stage  and  hence  engine  efficiency  overall.  The  maximum  re¬ 
duction  in  stagnation  pressure  drop  is  depicted  in  Figure  3b.  For  zero  losses  the  values  at  t3R  and 

t2R  are  coincident,  and  P^r  (also  Pmd  are  increased  to  ideal  values.  For  a  given  temperature  drop 
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across  the  LPT  stage,  a  higher  exit  pressure  (Pa)  for  the  same  inlet  pressure  (Pti)  will  correlate  with 
a  decrease  in  the  production  of  entropy.  This  can  be  represented  by  the  Gibbs  equation  (Mattingly 
[20]): 

HT  HP. 

(6) 


,  dT  dPt 

dS  —  Cn— - /t— — 

^  T  Pt 


where  R  is  the  gas  constant.  For  a  constant  temperature  equation  6  reduces  to: 

ds  _  dPt 

~R  =  Pt 

for  the  case  of  constant  temperature.  The  relation  between  pressure  drop  and  efficiency  for  the 


(7) 


turbine  stage  is  represented  by  (Mattingly  [20]): 


*/t  = 


1  —  r 


(7-1  )h 


1-7T, 

where  r  =  Ttz/Tt\  and  7rt  =  Pa/Pn  ■  The  pressure  Pa  is  a  function  of  the  loss  coefficient: 


Pa  = 


Pt2R 


(  T  \7/(7-l)] 

1  +  7  loss 

M*) 

(Ik.) 

VI'aiJ 


7/(7- 1) 


(8) 


(9) 


where  7 loss  is  defined  as  the  loss  coefficient  and  7  is  the  ratio  of  specific  heats.  Equation  8  shows 
that  for  a  given  temperature  ratio  across  the  blades,  a  larger  pressure  ratio  (77)  will  increase  the 
efficiency,  77.  The  magnitude  of  typical  LPT  stage  efficiencies  can  be  determined  from  numbers 
representative  of  a  jet  engine.  Assuming  initial  values  of  r  =  0.84,  Pt 3  =  748.1  kPa,  Pt2r  =  1043.0 
kPa,  T3  =  1505.5  K,  Thfl  =  1664.4  K,  Ta  =  1560.0  K,  and  7 loss  =  0.15  for  equations  9  and  8,  the 
baseline  loss  coefficient  is  r]t  =  90.87%.  As  7ioss  is  reduced  77  increases  and  improvement  in  stage 
efficiency  with  reduction  of  loss  coefficient  can  be  determined  from  equations  8  and  9.  Figure  5 
shows  that  a  10%  reduction  in  loss  coefficient  can  improve  stage  efficiency  by  0.57%.  The  resulting 
linear  curve  shows  that: 

'A7v 


—  ^  0.057  • 

n 


(10) 


where  ^  is  the  percent  increase  in  efficiency  and  ^  is  the  percent  reduction  in  loss  coefficient. 
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percent  improvement  of  stage  efficiency  (Ax\/r\) 


Figure  4.  T-S  diagram  showing  improvement  possible  with  decrease  in  pressure  drop. 


Figure  5.  Trend  of  increasing  stage  efficiency  with  reduction  in  loss  coefficient. 
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1.4.1  Comparison  of  Experimental  Loss  Coefficients 

Installation  of  a  passive  separation  control  method  into  an  existing  engine  would  change  the 
stagnation  pressure  drop  across  the  modified  blades.  The  change  in  pressure  drop  for  a  variety 
of  passive  techniques  can  be  compared  using  a  non-dimensional  form.  The  nondimensional  loss 
coefficient  can  be  defined  using  the  nomenclature  of  Figure  3: 


_  Ptm  —  Pur 
\pU%R 

where  Pt<i  and  £/2  are  measured  downstream  of  the  blade. 


(11) 


The  increase  in  exit  pressure  within  the  engine  due  to  the  reduction  in  separation  losses  is  also 
a  consideration.  For  ideal  operation,  the  engine  exit  static  pressure,  Pg  is  equal  to  the  atmospheric 
pressure,  Pg.  Using  the  station  numbers  in  Figure  6,  the  ratio  of  actual  thrust  (Fg)  to  ideal  thrust  can 
be  written  as: 


CF g— actual  ^  ^ 
fg  =  t; -  = 


F„- 


g— ideal 


(ft) 


(7—1  )h 


v-(ft) 


(7-l)/7 


1  + 


7-1 


1  _  E* 

_ _ Ps 


27 


(12) 


where  Cd  is  the  actual  to  ideal  mass  flow  rate  ratio  and  Cy  is  the  actual  to  ideal  velocity  ratio. 
Considering  the  recent  works  of  Murawski  et  al.  [4]  and  Qiu  and  Simon  [5]  which  have  shown  that 
the  pressure  losses  due  to  separation  within  the  turbine  stages  operating  at  low  Reynolds  numbers  are 
higher  than  expected,  the  pressure  ratio  Pg/Pg  would  no  longer  equal  unity.  Equation  12  indicates 
that  with  the  increased  pressure  losses  engine  thrust  is  decreased  for  an  engine  with  fixed  exit  area. 
The  exit  area  is  at  station  9  in  Figure  6a.  For  the  case  of  a  nozzle  correctly  expanding  the  flow,  the 


T-S/R  diagram  of  Figure  6b  is  useful.  The  relationship  between  temperature  and  pressure  is: 

-(7~l)/7 


Tg 


=(») 


Ttg 


(13) 


which  can  be  substituted  into  the  energy  (assuming  no  heat  addition  and  a  calorically  perfect  gas) 
equation: 

cp(Tt9-T9)  =  ^  (14) 
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to  yield: 


CpTn  ^1 


(15) 


where  an  increase  in  engine  pressure  will  lead  to  an  increase  in  exit  velocity.  Since  thrust  is  equal 
to  the  mass  flow  rate  times  the  exit  velocity,  an  increase  in  exit  velocity  increases  the  thrust  and 
thermal  efficiency  (c.f.,  equation  3). 


Figure  6.  Station  numbers  and  T-S/R  diagram  for  exit  nozzle. 

1.5  Chapter  Summary 

Chapter  2  includes  a  discussion  on  boundary  layer  development  and  previously  conducted 
experiments,  followed  by  a  description  of  the  wind  tunnel  used  in  this  research  in  Chapter  3 .  Chapter 
4  explains  the  data  collection  equipment  used  in  the  wind  tunnel.  Simulations  for  the  baseline 
performance  using  the  2-D  VBI  code  are  shown  in  Chapter  5 .  The  baseline  performance  of  the  Pak- 
B  blades  is  then  discussed  in  Chapter  6,  illustrating  the  severe  losses  caused  by  the  flow  separation. 
The  results  of  using  a  boundary  layer  trip  are  presented  in  Chapter  7.  Chapters  8  and  9  discuss  the 
improved  performance  of  the  Pak-B  blade  with  the  addition  of  dimples  or  V -grooves  to  the  surface 
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of  the  blade.  The  conclusions  gained  from  the  experiments  and  simulations  are  discussed  last  in 
Chapter  10. 
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Chapter  2  -  Background  and  Theory 


2.1  Overview  of  Turbine  Blade  Boundary  Layer  Aerodynamics 

Low  pressure  turbine  (LPT)  blades  have  been  studied  more  in  the  last  ten  years  than  in  all  the 
decades  previously.  Halstead  et  al.  [2, 3]  documented  the  flow  consists  of  low  to  medium  levels  of 
freestream  turbulence  (1%  to  20%),  passing  wakes,  centrifugal  effects  and  vortices.  The  unsteady 
nature  of  these  interactions  makes  LPT  flowfields  a  difficult  problem  to  examine. 

Experimental  studies  by  Werle  [11],  Mayle  [12]  and  Addison  and  Hodson  [21]  of  low  pressure 
turbine  blades  operating  at  low  axial  chord  Reynolds  numbers  indicate  a  gradual  transition  process 
from  laminar  flow  to  turbulent  flow  over  the  suction  side  of  the  blade.  The  approximately  steady 
transition  process  is  interrupted  by  the  time-dependent  passing  of  wakes  due  to  the  motion  of  the 
rotor  relative  to  the  upstream  stator.  This  cyclic  influence  quickly  trips  the  flow  through  a  fast 
succession  of  transitional,  turbulent,  and  back  to  transitional  state  for  a  small  region  as  the  wake 
travels  downstream  (Halstead  et  al.  [3]).  This  unsteady  phenomena  was  not  examined  since  the 
primary  goal  in  this  research  project  was  to  reduce  the  flow  separation  which  occurs  during  non¬ 
wake  periods.  The  results  gained  from  the  baseline  blade  experiments  are  valuable  to  programmers 
who  need  experimental  data  to  increase  the  accuracy  of  their  numerical  models  (Domey  [22]). 

2.2  Experimental  and  Numerical  Studies  of  Turbine  Boundary  Layers 

Attempts  to  numerically  model  the  complex  effects  in  the  low  pressure  turbine  have  met  with 
only  limited  success  due  to  the  lack  of  experimental  data  for  low  Reynolds  numbers  (Werle  [11]). 
To  build  better  computational  transition  models  for  the  low  Reynolds  number  boundary  layer  of  a 
LPT  airfoil,  detailed  boundary  layer  velocity,  turbulence,  and  pressure  measurements  are  needed 
for  determination  of  the  simulations  accuracy.  Halstead  et  al.  [2,23]  tested  four  different  numerical 
models;  SIANX,  KEP  a  Navier-Stokes  solver,  and  the  Fan-Lakshminarayana  code.  Each  of  the 
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four  codes  used  by  Halstead  et  al.  [2,23]  have  shown  less-than-desirable  accuracy  in  the  simulation 
of  low  Reynolds  number  characteristics.  A  code  that  produces  typically  good  outputs  except  in  a 
low  Reynolds  number  environment  is  the  2-D  Allison  Vane-Blade  Interaction  (VBI)  code.  The  2-D 
VBI  code  was  used  in  this  study  to  model  the  baseline  Pak-B  blade  cassette  for  reasons  discussed 
later.  As  with  other  codes  the  VBI  results  showed  the  code’s  inability  to  correctly  model  the  region 
of  separated  flow.  Like  most  codes,  the  2-D  VBI  code  was  incapable  of  simulating  the  3-D  effects 
of  the  passive  modifications  to  the  Pak-B  blades  to  improve  efficiency.  More  accurate  numerical 
models  were  not  available  for  this  research.  Halstead  and  Solomon  [24]  have  stated  that  perhaps 
only  a  full  Laige-Eddy  Simulation  (LES)  code  or  even  a  Direct  Numerical  Simulation  (DNS)  would 
provide  accurate  results  for  blades  with  modified  three  dimensional  surface  features. 

2.2.1  Experimental  Determination  of  Laminar-to-Turbulent  Transition  Modes 

Perhaps  the  most  important  goal  of  examining  turbomachinery  blades  is  the  determination  of 
the  location  and  length  of  transition  over  the  blade  surface.  This  goal  is  difficult  due  to  the  different 
possible  causes  of  flow  turbulence.  The  complexity  ofthe  mechanical  design  and  unsteady  flowfield 
are  also  complicating  factors  for  experimenters. 

Detailed  experiments  by  Mayle  [12]  and  Walker  [25]  showed  that  flow  can  enter  transition 
over  airfoil  surfaces  through  three  basic  means  depending  on  the  Reynolds  number,  freestream  tur¬ 
bulence  and  pressure  gradients.  The  three  means  are  natural  transition,  bypass  transition,  and  sepa¬ 
ration  bubble  transition  (Walker  [25]).  For  natural  transition  to  occur,  the  environment  must  consist 
of  low  freestream  turbulence  levels  and  no  adverse  pressure  gradients.  The  first  step  of  natural  tran¬ 
sition  is  the  formation  of  2-D  Tollmien-Schlichting  waves  which  amplify  instabilities  and  create  3-D 
loop  vortices  and  large  fluctuations  (Mayle  [12]).  These  3-D  disturbances  coalesce  into  turbulent 
spots  which  then  grow  as  they  travel  downstream.  Eventually  these  spots  merge  together  in  the  cre- 
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ation  of  a  fully  turbulent  boundary  layer  concluding  the  process  identified  by  Walker  [25]  as  natural 
transition. 

For  bypass  mode  transition  to  occur,  flow  conditions  with  freestream  turbulence  levels  of  ap¬ 
proximately  20%  are  necessary.  The  initially  laminar  flow  over  the  blade  will  immediately  form 
turbulent  spots  (Mayle  [12])  under  such  high  turbulence  levels.  Since  in  this  case  the  2-D  Tollmien- 
Schlichting  waves  are  not  found  the  term  “bypass  mode”  is  used  (Walker  [25]).  Whlker  [25]  deter¬ 
mined  that  linear  stability  theory  can  still  be  valuable  as  a  guide  for  bypass  mode  in  the  determina¬ 
tion  of  the  length  of  the  transitional  flow.  The  value  of  linear  stability  theory  is  due  to  the  fact  that 
although  the  2-D  Tollmien-Schlichting  waves  are  bypassed  in  this  mode,  the  disturbances  within  the 
boundary  layer  can  still  be  predicted  using  linear  stability  theory  (Walker  [25]).  The  application  of 
linear  stability  theory  is,  however,  limited  and  unable  to  provide  relevant  results  for  separated  flow 
transition. 


Figure  7.  Laminar  separation  transition  to  turbulent  flow. 
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For  adverse  pressure  gradients  the  laminar  flow  may  separate  forming  a  separation  bubble  (see 
Figure  7).  The  fluid  undergoes  the  transitional  process,  becoming  turbulent  while  separated  from 
the  blade  surface.  Once  the  separated  flow  is  fully  turbulent,  Qiu  and  Simon  [9]  and  Mayle  [12] 
showed  the  increased  energy  of  the  fluid  enhances  re-attachment  of  the  flow  to  the  surface  of  the 
blade.  It  is  this  form  of  laminar  flow  separation  and  turbulent  flow  reattachment  that  exists  over 
the  surface  of  die  Pak-B  blade  at  low  Reynolds  numbers.  This  research  focused  on  the  prevention 
of  the  laminar  separation  bubble  to  increase  efficiency.  The  computational  modelling  of  laminar 
separation  is  an  attempt  to  identify  conditions  for  bubble  transition,  but  is  very  difficult. 

2.2,2  Limitations  of  Numerical  Code  LPT  Modelling 

The  modelling  of  transition  from  laminar  to  turbulent  flow  is  difficult  to  accomplish  accurately. 
Freestream  turbulence  level,  Reynolds  number,  and  pressure  gradient  all  influence  the  predicted  re¬ 
sults.  To  compare  the  influence  of  these  different  freestream  factors,  Reed  etal.  [26]  examined  many 
methods  of  predicting  the  linear  stability  of  fluid  flows.  Most  methods  examined  by  Reed  et  al.  [26] 
required  some  prior  knowledge  of  the  flow  to  predict  the  desired  results.  The  eN  method  (Reed  et 
al.  [26]),  for  example,  requires  the  initial  disturbance  amplitude.  Sharma  et  al.  [27]  determined  that 
part  of  the  reason  numerical  turbulence  models  have  difficulty  in  predicting  flow  behavior  is  the 
nonisotropic  nature  of  turbulence  within  some  laminar  and  transitional  boundary  layers.  Mueller 
[28]  found  that  the  accuracy  from  some  software  packages  has  been  so  poor  that  the  software  results 
would  not  only  hinder  engine  design,  but  decrease  engine  performance.  Much  of  this  inaccuracy  is 
due  to  the  difficult  nature  of  predicting  the  evolution  and  location  of  separation  bubbles  that  form 
on  the  suction  surface  of  the  blade.  These  bubbles  may  completely  separate  and  cause  high  losses 
for  the  blade  at  low  Reynolds  numbers  (Werle  [11]  and  Mayle  [12]).  Other  codes  for  predicting  the 
separation  losses  were  investigated  by  Halstead  et  al.  [23]  and  are  discussed  further  in  Appendix  K. 
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The  Allison  2-D  \hne  Blade  Interaction  (VBI)  code  was  chosen  for  use  in  this  research  because  of 
its  efficient  combination  of  airfoil  wrapping  O-grids  and  passage  filling  H-grids.  The  combination 
of  the  O  and  H  grids  in  the  VBI  code  provided  more  accurate  results  than  would  otherwise  have 
been  possible. 

The  2-D  VBI  code  is  capable  of  modelling  the  flow  of  a  low  pressure  turbine  cascade  and 
was  used  to  provide  initial  surface  pressure  coefficient  data  for  inlet  chord  Reynolds  numbers  of 
50k,  100k,  and  200k.  While  the  2-D  VBI  code  predicted  a  small  separation  bubble  at  a  Reynolds 
number  of  50k,  the  code  failed  to  predict  separation  bubbles  for  Reynolds  numbers  of  100k  and  200k. 
Detailed  information  for  the  Allison  2-D  VBI  code  such  as  user  inputs  and  simulated  boundary  layer 
velocity  profiles  are  included  in  Chapter  5  and  Appendix  E.  The  results  presented  using  the  2-D 
VBI  code  illustrate  the  shortcomings  even  the  advanced  VBI  code  has. 

Several  authors  including  Rivir  [8],Werle  [11],  and  Mayle  [12]  have  expressed  the  need  for 
more  experimental  data  to  fully  explain  the  flowfield  within  the  LPT  The  works  of  Halstead  [13] 
and  Addison  and  Hodson  [2 1  ]  have  provided  a  significant  increase  in  the  understanding  of  unsteady 
wake  phenomena  and  the  existence  of  large  regions  of  laminar  and  transitional  flow;  however  much 
remains  to  be  documented  for  low  Reynolds  number  flowfields.  Additional  experimental  data  of 
boundary  layer  growth,  turbulence  intensities  within  the  boundary  layer,  and  turbulent  length  scales 
are  all  desired  (Rivir,  [29])  in  understanding  the  problem  of  laminar  separation  bubbles. 

2.3  Separation  Control  Using  Passive  Methods 

As  shown  in  many  experiments  performed  by  Murawski  et  al.  [4]  and  Qiu  and  Simon  [9],  sep¬ 
aration  is  more  apt  to  occur  by  decreasing  either  the  axial  chord  Reynolds  number  below  300,000 
or  decreasing  the  freestream  turbulence  levels  below  10%,  conditions  under  which  low  pressure 
turbine  blades  frequently  operate  at  altitude  and  cruise  conditions.  As  the  Reynolds  number  is  de- 
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creased,  the  transition  location  moves  aft  towards  the  trailing  edge  on  the  suction  surface  of  the 
blade.  This  effect  combined  with  the  lower  freestream  turbulence  level  (which  prevents  the  flow 
from  transitioning  to  a  turbulent  state)  can  combine  to  create  a  laminar  separation  bubble  before 
transition.  The  laminar  separation  bubble  creates  losses  which  reduce  the  engine  operating  effi¬ 
ciency.  Previous  work  by  Murawski  et  al.  [4]  and  Qiu  and  Simon  [9]  with  the  Pak-B  blades  have 
exhibited  this  laminar  separation  problem.  The  Pak-B  blades  were  therefore  ideally  suited  for  fur¬ 
ther  in-depth  research  into  the  control  of  separated  flow.  Reduction  of  the  separation  problem  over 
the  Pak-B  blades  provides  the  community  with  a  better  insight  into  proper  blade  design  and  de¬ 
scribes  methods  by  which  flows  can  be  modified  and  controlled  within  the  LPT  The  experimental 
work  was  completed  in  two  parts. 

The  first  goal  of  this  study  was  to  obtain  boundary  layer  velocity  and  turbulence  profiles, 
surface  pressure  coefficients,  and  total  pressure  loss  data  for  the  unmodified  Pak-B  blade  at  low 
Reynolds  numbers  of  50k,  100k,  and  200k  and  turbulence  intensities  of  1%  and  4%.  The  experi¬ 
mental  boundary  layer  data  was  used  to  determine  the  location  of  separation  on  the  suction  surface 
of  the  Pak-B  blade.  Attention  was  given  to  verify  that  inlet  Mach  numbers  were  appropriate  and 
relevant  for  the  experiments  conducted.  The  total  pressure  loss  coefficient  (7)  and  wake  momen¬ 
tum  deficit  (©)  measurements  indicated  the  efficiency  of  the  baseline  Pak-B  blade.  A  listing  of  the 
types  of  experimental  data  collected  for  the  baseline  Pak-B  blade  is  provided  in  Table  1. 


Table  1 .  Baseline  Pak-B  data  collected 


Reynolds  number 

Cp 

B.L.  Profiles 

7 

© 

50k 

X 

X 

X 

X 

100k 

X 

X 

X 

X 

200k 

X 

X 

X 

The  second  goal  of  this  study  was  to  examine  the  baseline  Pak-B  experimental  data,  including 
the  recorded  chordwise  locations  of  separation  and  reattachment,  for  possible  means  of  reducing  the 
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losses  associated  with  laminar  separation.  Study  of  the  Pak-B  experimental  data  and  the  work  of 
several  separation  control  researchers  suggested  three  passive  design  modifications  to  the  suction 
surface  likely  to  have  the  most  beneficial  impact.  The  first  modification  previously  documented  by 
Bearman  and  Harvey  [16, 17]  involved  the  use  of  dimples  recessed  into  the  blade  surface.  These 
dimples  were  placed  spanwise  over  a  portion  of  the  blade.  As  suggested  by  Bearman  and  Harvey 
[16, 17]  these  dimples  create  vortices  within  the  boundary  layer  which  eneigize  the  flow  from  a 
baseline  laminar  state  to  a  transitional  one  before  separation  can  occur.  This  increase  in  skin  friction 
reduced  the  2-D  profile  losses  associated  with  separation. 

The  second  design  modification  tested  was  the  placement  of  streamwise  longitudinal  V-grooves 
on  the  suction  surface  of  the  Pak-B  blade.  Previously  test  by  Lin  et  al.  [18]  for  separation  control 
over  a  convex  surface  in  a  turbulent  boundary  layer,  these  grooves  are  similar  to  riblets,  but  much 
larger.  The  increased  size  of  the  V-grooves  enhanced  the  formation  of  vortices  extending  longitudi¬ 
nally  downstream  for  a  turbulent  boundary  layer  and  subsequently  mixed  higher  energy  flow  from 
the  ffeestream  with  the  lower  energy  boundary  layer  to  reduce  separation  (Yurchenko  [30]). 

The  third  modification  was  the  placement  of  a  trip  wire  on  the  surface  of  the  Pak-B  blade  at  the 
35%  axial  chord  position.  This  trip  wire  created  a  state  of  turbulent  flow  over  the  suction  side  of  the 
blade,  \ferification  was  made  that  the  flow  remains  attached  to  the  suction  surface  for  all  Reynolds 
numbers  tested  when  atrip  is  used.  The  loss  coefficient  was  greatly  reduced  at  a  Reynolds  number 
of  50k;  however,  at  higher  Reynolds  numbers  of  100k  and  200k  the  loss  coefficient  for  the  blade 
with  the  wire  trip  was  higher  than  the  baseline  blade. 

The  experiments  with  the  modified  blades  were  done  in  an  adjustable  cascade  tunnel  which 
was  capable  of  turning  the  flow  through  the  95°  necessary  to  simulate  the  operational  environment 
of  the  Pak-B  blades.  The  same  turbulence  levels  and  Reynolds  numbers  used  in  initial  tests  with  an 
unmodified  baseline  Pak-B  blade  were  repeated  for  the  three  modified  blades.  The  tests  included 


21 


measurement  of  the  surface  pressure  coefficient,  boundary  layer  profiles,  total  pressure  loss  coeffi¬ 
cient  (7),  and  the  momentum  wake  deficit  (©).  In  summary,  these  tests  show  that  for  aturbine  blade 
with  a  high  degree  of  turning  and  low  Reynolds  number  flow  separation  problems,  passive  surface 
modifications  will  reduce  the  cascade  flow  losses. 
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Chapter  3  -  Tunnel  and  Cascade  Characteristics 


3.1  Data  Collection  Locations 

The  wind  tunnel  adjustable  test  section  was  designed  to  accommodate  several  different  types 
of  data  collection.  The  tunnel  has  four  thermocouples  suspended  from  the  tunnel  ceiling  for  recov¬ 
ery  temperature  measurement  marked  T-l  through  T-4  (see  Figure  8).  The  first  thermocouple  is 
located  100cm  (39.4  inches)  up  from  the  bottom  and  61  cm  (24  inches)  from  the  inside  wall  on  the 
304.8  cm  (120  inch)  wide  by  266.7  cm  (105  inch)  high  inlet  flow  straightener  screen.  The  second 
thermocouple  is  located  45 .7  cm  ( 1 8  inches)  downstream  of  the  turbulence  generating  grid  and  1 .45 
m  (57  inches)  upstream  of  blade  #1.  Thermocouple  #3  is  located  on  the  40%  axial  chord  line  of 
blade  Ml.  The  last  thermocouple  is  76.2  cm  (30  inches)  downstream  of  blade  #1. 


AIRFLOW  EXIT  DIFFUSER 


Figure  8.  Measurement  collection  locations  for  the  adjustable  test  section. 
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The  blade  cassette  consists  of  8  Pak-B  blades  each  approximately  seven  times  engine  scale. 
Each  test  blade  has  an  axial  chord  length  of  17.78  cm  (7  inches)  with  an  axial  chord  to  spacing  ratio 
(solidity)  of  1.129.  To  enhance  the  approximately  2-D  flowfield  of  the  mid-span  of  the  blade,  each 
test  blade  is  stretched  in  the  z-axis  to  78.63  cm  (34.5  inches)  for  an  aspect  ratio  of  4.92  to  1.  The 
test  blades  were  molded  from  Ultralloy  108  white  resin  and  placed  inside  the  test  section  between 
traverse  slots  #2  and  #3. 

Four  traverse  slots,  S-l  to  S-4,  are  located  in  the  plexiglass  top  of  the  test  section.  The  slots 
are  for  hot  wire  and  pressure  probe  access  into  the  test  section  flowfield.  The  hot  wire  and  pres¬ 
sure  probe  were  used  to  record  flowfield  variations  across  the  pitch  of  the  blades.  The  first  slot  is 
upstream  of  the  blades  and  marked  S-l  in  Figure  8.  The  centerline  of  the  1.27  cm  (0.5  inch)  wide 
slot  S-l  is  approximately  one  axial  chord  length  upstream  of  the  blades  (19.69  cm)  while  slot  S-2  is 
approximately  \  chord  length  (9.53  cm)  upstream  of  the  blades.  Slots  S-3  and  S-4  are  downstream 
of  the  blade  cassette  approximately  |  chord  and  1  axial  chord  length  downstream  respectively.  The 
precise  location  of  S-3  and  S-4  are  10.8  cm  (4.25  inches)  and  21.34  cm  (8.4  inches)  downstream, 
respectively. 

A  pitot-static  probe  is  located  19.69  cm  (7.75  inches)  upstream  of  blade  #5  to  provide  pressure 
for  three  different  measurements.  The  pitot-static  probe  was  used  to  provide  inlet  total  pressure 
readings  when  collecting  data  from  the  surface  pressure  taps  on  blades  #4  and  #6.  The  pressure 
probe  also  provided  an  upstream  pressure  reference  for  the  exit  total  pressure  measurements  neces¬ 
sary  for  calculation  of  the  loss  coefficients.  Finally,  the  pressure  probe  was  used  to  determine  the 
inlet  velocity  and  Reynolds  number  when  the  hot  wire  was  being  used  for  velocity  and  turbulence 
measurements  downstream  of  the  blade  cassette.  The  inlet  and  exit  angle  of  the  flow  through  the 
blade  cassette  was  determined  through  a  pressure  probe  similar  to  the  pitot-static  probe. 
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3.2  Flow  Incidence  Angle  Measurements 

To  ensure  proper  modelling  of  the  engine  environment,  the  flow  inlet  and  exit  angles  were 
measured  using  a  3-hole  wedge  probe  (see  Figure  9).  The  two  side  ports  on  the  30°  faces  were 
connected  to  the  two  ports  of  a  Dwyer  Instruments  5.08  cm  (2  inch)  inclined  manometer.  A  reading 
of  zero  ±0.005  U1H2O  on  the  manometer  indicated  the  probe  was  facing  directly  into  the  flow.  The 
low  error  in  the  manometer  measurement  allowed  measurement  of  the  flow  angle  to  within  ±0.25°. 

flow  direction 


Figure  9.  Three  hole  wedge  probe  used  in  inlet  angle  measurements. 

The  Pak-B  low  pressure  turbine  blade  design  has  a  camber  angle  (angle  on  centerline  between 
leadine  and  trailing  edges  of  blade)  of  36°,  a  stagger  angle  (angle  between  centerline  of  blade  and 
axial  direction)  of  25°  and  an  incidence  of-l°  for  the  50,000  Reynolds  number  case.  The  measured 
experimental  inlet  and  exit  flow  angles  are  both  shown  in  Figure  10.  The  measured  inlet  angle 
values  are  near  the  design  value  of  35° .  The  exit  flow  angles  are  also  presented  in  Figure  1 0  and  are 
near  their  design  value  of  60°.  The  measured  angles  for  the  higher  turbulence  case  of  4%  created 
by  the  turbulence  grid  are  depicted  as  solid  filled  symbols.  The  traverse  distance  in  centimeters  is 
shown  on  the  lowerx-axis  ofFigure  10.  The  numbers  across  the  top  x-axis  of  Figure  10  represent 
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the  inlet  tangent  points  of  the  leading  edge  blade  circle  to  the  inlet  plane  for  each  of  the  eight  blades 
in  the  cassette.  The  exit  angles  in  Figure  10  correspond  with  the  trailing  edge  of  the  blades  labeled 
across  the  top  of  the  figure.  The  Reynolds  number  case  of  50,000  shows  a  non-uniformity  which 
was  necessary  to  get  the  Cp  distributions  for  both  blade  #4  and  #6  close  to  the  desired  operational 
setting.  The  exit  angles  for  a  Reynolds  number  of  50,000  and  freestream  turbulence  of  1%  show  a 
±3°  variation  which  was  indicative  of  the  laige  separation  regions  over  the  suction  surface  of  the 
blades  at  this  condition.  Table  2  provides  a  comparison  of  the  inlet  angles  between  test  blade  #5  and 
the  2-D  VBI  computational  model.  The  stagger  and  camber  angle  are  25°  and  36°,  respectively,  for 
both  the  experimental  and  computational  runs. 


Table  2.  Flow  inlet  angles  for  blade  5  and  2-D  VBI  software  model 


Reynolds  number 

Inlet  angle  -  Tunnel 

Inlet  angle  -  VBI  code 

50k 

35.5° 

36° 

100k 

36.5° 

35° 

3.3  Periodicity  and  Uniformity  Measurements 

Air  velocity  magnitude  and  turbulence  intensity  were  measured  across  the  blade  row  inlet  and 
exit  with  a  TSI  1210-T1.5  hot  wire  probe  to  ensure  periodicity  and  uniformity.  Periodicity  is  the 
regulamess  of  the  velocity  peaks  and  valleys  caused  by  the  air  flow  moving  around  objects  within 
the  test  section.  Good  periodicity  downstream  of  the  blade  cassette  would  have  each  peak  15.75 
cm  (6.2  inches)  apart  corresponding  to  the  distance  separating  each  blade.  Uniformity  is  defined 
as  the  consistency  of  the  mean  velocity  across  the  tunnel.  For  good  uniformity,  each  velocity  peak 
would  have  the  same  value,  and  each  velocity  valley  would  have  a  constant  and  smaller  value.  The 
measurements  were  made  across  a  60  cm  (23.62  inch)  span  from  blade  #3  to  blade  #7.  The  hot  wire 
probe  was  moved  by  a  DanTec  3 -axis  traverse  through  the  traverse  slots.  To  verify  2-D  flow  across 
the  midspan  of  the  test  blades,  velocity  and  turbulence  measurements  were  made  at  five  different 
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#1  #2  #3  #4  #5  #6  #7  #8  Outer  Wall 

65  t - : - 1 - 1 - 1 - : - 1 - 1 - 1 - 1 


traverse  (cm) 


Figure  10.  Measured  inlet  and  exit  flow  angles. 
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heights  within  the  test  section.  The  hot  wire  was  moved  across  the  blade  row  inlet  and  exit  10  cm 
below  the  center  span,  5  cm  below,  across  the  center  span,  and  at  +5  and + 1 0  cm  above .  This  method 
is  shown  in  Figure  11 . 


Test  Section  Ceiling 


Test  Section  Floor 


Figure  11.  Traverse  movement  pattern  to  verify  2-D  flow  with  a  hot  wire  probe. 

Each  velocity  measurement  was  made  three  times  with  the  data  then  ensemble  averaged.  The 
averaged  measured  velocity  results  are  presented  in  Figures  12through  15.  The  figures  each  contain 
velocities  for  all  three  Reynolds  numbers  used  in  this  experiment.  The  velocities  shown  in  Figures 
12  through  15  illustrate  the  periodicity  and  uniformity  of  the  flow  by  having  no  major  differences 
between  graphs.  The  velocity  information  collected  with  the  hot  wire  was  also  used  to  determine 
the  length  scales  of  the  turbulence  within  the  test  section  inlet. 
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3.4  Test  Section  Inlet  Turbulence  Scale  Measurements 


To  calculate  the  length  scales  of  the  turbulence  within  the  tunnel  test  section,  velocity  data  was 
collected  with  the  Auto-Correlation  VI  (see  Appendix  D.9)  software  using  the  TSI  1210-T1.5  hot 
wire  probe.  The  hot  wire  probe  was  inserted  19.69  cm  (7.75  inches)  upstream  of  the  blade  cassette. 
The  software  collected  216  velocity  sample  measurements  at  20,000  Hz  which  resulted  in  a  duration 
of  3.28  seconds.  The  high  data  collection  frequency  was  chosen  to  capture  all  of  the  frequencies 
which  exist  in  a  low  subsonic  freestream.  The  power  spectrum  of  the  collections  is  shown  in  Figures 
16  and  17  for  Reynolds  numbers  of  50k,  100k,  and  200k. 

The  y-axis  of  Figures  16  and  17  are  the  single-sided  auto  power  spectrum  of  the  rms  velocity 
component  squared  ([ti']2)  divided  by  the  period  of  the  sampling  signal  (1/20,000  Hz).  This  division 
by  the  sampling  period  results  in  the  y-axis  having  units  of  m2/sec.  When  the  curves  present  in 
Figure  17  were  integrated,  the  resulting  value  was  the  square  of  the  rms  velocity  ([u']2)  for  that 
Reynolds  number  and  turbulence  intensity.  For  example,  the  power  spectrum  for  Re  =  50k  and 
Tu  =  4%  when  integrated  (using  the  trapezoidal  rule)  resulted  in  a  rms  velocity  value  of  0.129 
m/s  compared  to  the  direct  measurement  of  0.124  m/s.  The  length  scales  for  an  inlet  freestream 
turbulence  level  of  4%  were  determined  to  be  2.87  cm  (1.13  inches),  4.50  cm  (1.77  inches),  and 
6.93  cm  (2.73  inches)  for  Reynolds  numbers  of  50k,  100k,  and  200k,  respectively.  The  software 
used  in  the  determination  of  the  length  scale  is  discussed  in  Appendix  D. 

The  lower  freestream  turbulence  curves  of  Figure  16  do  not  correlate  as  well  to  the  direct  mea¬ 
surements  of  rms  velocity  due  to  a  low  signal  to  noise  ratio  for  the  Tu  =  1%  case.  When  integrated 
using  the  trapezoidal  rule,  the  area  under  the  curves  of  Figure  16  result  in  an  inaccurate  estimate  of 
the  square  of  the  rms  velocity.  This  was  evident  with  the  measurements  for  Re  =  50k  and  Tu  =  1% 
where  the  power  spectrum  integrated  to  an  estimated  rms  velocity  of  0.044  m/s  compared  to  the  di¬ 
rect  measurement  of  0.020  m/s.  With  the  limitations  for  the  Tu  =  1%  case,  the  length  scales  of  the 
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turbulence  for  the  inlet  flowfield  at  Reynolds  numbers  of  50k,  100k,  and  200k  were  determined  to 
be  77.3  cm  (30.43  inches),  1.231  m  (48.46  inches),  and  1.533  m  (60.35  inches),  respectively.  The 
large  length  scales  are  indicative  of  the  size  of  the  test  section.  The  equipment  used  for  determina¬ 
tion  of  the  velocity  is  explained  further  in  Chapter  4. 
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Figure  16.  Power  spectrum  of  inlet  flowfield  for  Tu  =  1%. 
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Figure  17.  Power  spectrum  of  inlet  flowfield  for  Tu  =  4%. 
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Chapter  4  -  Data  Collection  Instruments 


4.1  Velocity  Measurements  Using  TSI  1210-T1.5  Hot  Wire 

A  TSI  model  1210-T1.5  single  hot  wire  was  used  for  making  velocity  magnitude  measure¬ 
ments  upstream  and  downstream  of  the  test  blades.  Moving  the  hot  wire  pitchwise  in  the  traverse 
slots  allowed  the  uniformity  and  periodicity  of  the  velocity  into  and  exiting  the  blade  cassette  to  be 
measured.  Several  pieces  of  equipment  were  used  together  to  facilitate  the  hot  wire  data  measure¬ 
ments. 

The  hot  wire  sensor  was  moved  using  a  3-axis  Dantec  model  IMC-2  traverse.  The  traverse  is 
capable  of  moving  60  cm  (23.62  inches)  in  each  axis  which  provides  access  to  the  flow  properties 
±30.0  cm  (±11.81  inches)  to  either  side  of  blade  #5.  Y-axis  (pitch  direction)  probe  location  was 
measured  with  a  Sony  SR50-075A  linear  scale  and  displayed  using  a  Sony  IY-51  digital  display. 
This  display  provided  traverse  position  information  accurate  to  ±  0.005cm  (±  0.00197  inches). 
Using  the  traverse  to  hold  the  probe,  the  hot  wire  was  calibrated  while  in  the  tunnel. 

The  hot  wire  was  calibrated  using  a  pitot  static  air  pressure  probe  co-located  approximately 
1  axial  chord  length  upstream  of  blade  #5  with  approximately  10.16  cm  (4  inches)  of  space  be¬ 
tween  the  probes.  The  pitot  probe  was  connected  to  a  5.08  cm  (2  inch)  Dwyer  Instruments  inclined 
manometer.  Software  was  used  to  match  hot  wire  voltage  output  to  tunnel  pitot-static  pressures  us¬ 
ing  the  equipment  illustrated  in  Figure  18.  The  calibration  was  completed  with  a  non-dimensional 
curve  fit  between  Nusselt  Number  and  Reynolds  number  of  the  form  (Braun  [3 1]): 

Nu  =  C{nt  +  Dsi0  ■  Re  0  45  (16) 

where  Cint  and  Dsio  are  the  intercept  and  slope,  respectively,  of  the  calibration  curve.  The  exponen¬ 
tial  value  of  0.45  was  chosen  as  it  provided  an  acceptable  curve  fit  for  measuring  low  speed  flow. 


36 


As  shown  in  Figure  18,  the  hot  wire(b)  can  move  across  the  pitch  axis  of  the  blade  row(a) 
through  the  traverse  slot(c).  The  hot  wire  was  powered  by  an  IFA-100  Intelligent  Flow  Analyzer(d) 
which  provided  a  milli-volt  output  to  a  National  Instruments  SCB-68  data  collection  card  located 
within  the  wind  tunnel  control  room  pc.  The  total-static  pressure  probe(e)  provided  Ptotai  —  Pstatic 
for  the  air  flow  inside  the  tunnel  test  section.  A  5.08  cm  (2  inch)  inclined  manometer  was  used  to 
read  the  pitot-static  probe  pressure  differential.  A  transducer(f)  could  also  have  been  used  to  read 
the  total  probe.  The  HP  3852A(g)  collects  the  temperature  readings  and  passes  the  data  from  the 
J-type  thermocouple  within  the  tunnel  to  the  computer(h).  The  calibration  software  is  discussed  in 
further  detail  in  Appendix  D.  The  method  for  calibrating  the  TSI 12 10-T1 .5  hot  wire  probe  was  also 
used  to  calibrate  the  TSI  1276-10A  subminiature  hot  film  probe  using  the  same  pitot  probe,  inclined 
manometer  and  software. 


Figure  18.  Hot  wire  data  collection  diagram. 


4.2  Boundary  Layer  Velocity  Measurements  Using  TSI  1279S-10A  Miniature 
Hot  Film  Probe 

The  TSI  1279S-10A  hot  film  probe  is  a  modified  version  of  the  TSI  1276-10A  subminiature 
probe.  Unlike  the  straight  TSI  1276-10A,  the  filament  support  prongs  of  the  TSI  1279S-10A  are 
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turned  into  the  airflow  45°as  shown  in  Figure  19.  The  hot  film  probe  was  mounted  on  a  micro¬ 
traverse  mechanism  to  provide  accurate  boundary  layer  measurements  of  the  test  blade.  Five  loca¬ 
tions  along  the  suction  surface  of  the  test  blade  were  chosen  for  boundary  layer  measurements.  To 
maintain  structural  stability  for  the  small  probe  additional  hardware  was  constructed. 

The  hot  film  probe  was  supported  using  a  3.175  mm  (0.125  inches)  diameter  tube  that  was 
10.16  cm  (4  inches)  long.  This  small  probe  holder  extends  the  length  of  the  short  probe  while 
minimizing  aerodynamic  interference  within  the  blade  passage.  The  hot  film  probe  and  the  small 
probe  holder  were  supported  by  a  9.53  mm  (0.375  inches)  diameter  tube  43.18cm  (17  inches)  long. 
This  larger  tube  provided  the  additional  length  necessary  to  reach  from  the  traverse  mechanism  to 
the  blade  passage.  The  longer  tube  was  connected  to  the  micro  traverse  through  a  rectangular  slot  cut 
through  the  outer  tail  board.  The  micro  traverse  mechanism  is  in  the  non-moving  air  space  between 
the  outer  tail  board  and  outer  test  section  wall. 


film  element 


a)  top 


b)  side 


Figure  19.  TSI  1279S-10A  subminiature  boundary  layer  probe  (not  to  scale). 


The  National  Aperture,  Inc.  motorized  5.08  cm  (2  inch)  MM-3M  micro  traverse  has  a  move¬ 
ment  resolution  of  20157.4  revolutions  per  cm  (51199.8  turns  per  inch).  The  traverse  position  en¬ 
coder  provides  a  linear  spatial  resolution  of  ±5.08  x  10-5  cm  (±0.00001  inches).  The  traverse 
position  was  controlled  and  powered  by  a  National  Aperture,  Inc.  MC-3SA  servo  amplifier  system 

connected  to  a  National  Instruments  PC-Servo-2A  motion  control  board  located  within  the  wind 

38 


tunnel  control  room  computer.  The  operating  software  is  described  in  more  detail  in  Appendix  D. 
The  software  in  combination  with  the  high  position  accuracy  allow  very  detailed  measurements  of 
the  boundary  layer. 

The  boundary  layer  probe  was  used  to  measure  velocity  and  turbulence  profiles  over  five  dif¬ 
ferent  surface  locations  on  the  baseline  smooth  blade  and  the  modified  blades.  As  seen  in  Figure  20, 
the  measurement  locations  span  the  region  of  separated  flow  over  the  suction  surface  from  67.2% 
axial  chord  to  89.8%.  The  dimple  center  and  V-groove  start  locations  were  upstream  of  the  points 
where  the  boundary  layer  was  measured. 
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Figure  20.  Locations  of  boundary  layer  measurements  over  the  suction  surface. 


4.3  Temperature  Measurements  Using  J-Type  Thermocouples 

Temperature  measurements  within  the  tunnel  were  made  using  four  J-type  thermocouples.  The 
first  thermocouple  was  placed  at  the  bell  mouth  inlet  with  the  second  thermocouple  located  1.45  m 
(57  inches)  upstream  of  blade  #1.  The  third  thermocouple  was  placed  25.4  cm  (10  inches)  down 
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from  the  top  of  blade  #1  at  the  40%  axial  chord  point.  The  final  thermocouple  was  located  76.2 
cm  (30  inches)  downstream  of  blade  #1.  A  Kaye  Instruments  Ice  Point  Thermocouple  Reference 
provided  a  0°C  (32°F)  reference  and  a  Hewlett  Packard  3852A  Data  Acquisition/Control  Unit  using 
an  HP  44708F  (20  channel  solid-state  multiplexer:  thermocouple  compensation)  card  was  used  to 
capture  the  voltages.  The  voltages  were  then  measured  by  an  HP  44701 A  (5  5  digit  voltmeter)  card 
within  the  HP  3852A  Data  Acquisition/Control  Unit.  The  voltages  were  then  relayed  to  the  wind 
tunnel  data  collection  computer  via  GP-IB  data  connection.  The  thermocouples  are  accurate  to 
within  ±0. 11  °K(±  0.2  0  F).  Air  density  was  calculated  using  the  measured  temperature  along  with 
the  atmospheric  pressure  of  the  wind  tunnel  laboratory. 

4.4  Pressure  Measurements 

Different  instruments  were  used  to  measure  the  various  pressure  readings  necessary  for  the  col¬ 
lection  of  data.  Room  atmospheric  pressure  was  read  from  a  Wallace  and  Tleman  FA- 129  pressure 
gauge  which  is  accurate  to  ±  3.5  kg/m2  (±  0.005  psi).  Pitot  pressure  readings  used  in  hot  wire  cali¬ 
bration  were  taken  using  a  5.08  cm  (2  inch)  Dwyer  Instruments,  Inc.  inclined  manometer.  Exit  total 
pressure  and  surface  static  pressure  readings  in  the  tunnel  were  measured  using  Druck  Incorporated 
351.6  kg/m2  (0.5  psi)  differential  pressure  transducers. 

Calibration  of  the  transducers  was  completed  by  using  a  Druck  Incorporated  DPI  600  Digi¬ 
tal  Pressure  Indicator  and  pressure  source  along  with  calibration  software  described  in  Appendix 
D.  After  calibration  the  transducers  have  an  accuracy  of  ±  1.25  kg/ms2  (±  0.005  in-H20).  The 
accuracy  of  the  calibration  degrades  with  changes  in  transducer  temperature.  Temperature  changes 
shift  the  zero  setting  in  the  transducer  calibration.  To  prevent  zero  shift  in  transducer  calibration, 
the  transducers  were  placed  inside  an  insulated  box.  Temperature  changes  within  the  transducer  box 
were  controlled  by  an  Omega  Engineering,  Inc.  model  6102  Temperature  Controller  which  operated 
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an  Omegalux  silicone  rubber  fiberglass  insulated  flexible  heater  placed  within  the  insulated  trans¬ 
ducer  box.  The  controller  maintained  the  temperatures  at  a  desired  temperature  setting  to  within  ± 
0.56°C  (1.0°F). 

4.4.1  Blade  Surface  Static  Pressure  Measurements 

Blade  surface  pressure  profile  measurements  were  made  using  0.76  mm  (0.03  inch)  diameter 
surface  pressure  taps  leading  to  tubes  within  the  blades.  Ihe  pressure  tubes  within  the  blade  were 
cast  in  during  the  molding  process  using  templates.  The  templates  were  cut  from  a  sample  blade 
into  0.794  cm  (0.3125  inch)  wide  segments.  Holes  were  then  drilled  into  the  templates  to  support 
the  pressure  tap  tubes.  Three  templates  were  used  in  each  pressure  tap  blade  to  provide  adequate 
support  of  the  tubes  over  the  span  of  the  blade  until  the  tubes  exit  the  blade  on  one  end.  Blade  #4 
has  37  working  taps,  while  blade  #6  has  40  usable  taps.  Problems  with  manufacturing  clogged  taps 
#0,  #17,  and  #21  on  blade  #4.  The  pressure  tap  hole  pattern  is  spread  across  the  center  22.86  cm 
(9  inches)  span  of  the  surface  of  the  blade.  The  hole  locations  were  staggered  at  an  angle  of  30°  to 
prevent  the  previous  hole  from  interfering  with  it.  The  experimental  data  was  collected  from  these 
pressure  taps  using  several  different  pieces  of  equipment  and  a  computer. 

The  pressure  tap  locations  are  shown  in  Figure  21.  The  taps  begin  with  #0  near  the  trailing 
edge  and  proceed  up  to  #12  near  the  leading  edge  for  a  total  of  12  taps  on  the  pressure  side.  The 
thirteenth  tap  was  located  near  the  leading  edge  stagnation  point,  while  taps  #14  through  #39  were 
on  the  suction  surface.  The  pressure  tap  layout  provided  detailed  measurement  of  the  region  of 
separation  in  the  aft  50%  of  the  suction  surface. 

Figure  22  illustrates  the  surface  pressure  tap  data  collection  process.  The  pressure  tubes  from 
the  surface  ports  on  the  blade(a)  are  connected  to  a  Scanivalve  Corporation  48-channel  48J9-1  mul¬ 
tiplexer^).  The  current  tap  read  by  the  transducer  is  indicated  on  a  Px  Scanner  Position  Display  (c). 
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The  multiplexer  power  is  supplied  by  a  CTRL  10P/S2-S6  Solenoid  Controller(d)  and  remotely  op¬ 
erated  using  an  HP  44728A  (8  channel  relay  actuator)  card  within  the  HP  3852A  Data  Acquisi¬ 
tion/Control  Unit(g).  Each  blade  with  pressure  taps  is  connected  to  its  own  scanivalve,  the  output 
of  each  is  connected  to  a  Druck  Incorporated  351.6  kg/m2  (0.5  psi)  pressure  transducer(f).  The 
transducers  used  in  taking  data  from  the  surface  pressure  taps  are  differential  transducers,  allowing 
the  total  port  of  the  inlet  pitot  probe(e)  to  provide  the  high  side  pressure  while  the  pressure  tap  sup¬ 
plies  the  low  side.  After  the  pressure  tap  data  is  acquired,  a  reading  of  zero  pressure  differential  is 
read  on  the  transducers  to  ensure  minimal  thermal  shift  in  the  results  with  any  result  taken  out  from 
the  acquired  data  set.  The  voltage  output  of  the  transducers  is  read  by  the  HP  3852A  Data  Acqui¬ 
sition/Control  Unit(g)  using  an  HP  4471 1A  (24  channel  high  speed  FET  multiplexer)  card  with  the 
output  relayed  to  the  data  collection  pc(h). 


Figure  21.  Surface  static  pressure  port  locations  molded  within  blades  #4  and  #6. 

4.4.2  Exit  Total  Pressure  Measurements 

A  pitot  probe  was  used  to  record  the  total  pressure  in  the  downstream  wake  of  the  test  blades. 
The  downstream  pitot  probe  was  connected  to  one  side  of  a  Druck  Incorporated  351.6  kg/m2  (0.5 
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psi)  differential  transducer  while  the  upstream  pitot  probe  was  connected  to  the  other.  The  usage 
of  two  pitot  probes  allowed  direct  measurement  of  Po, inlet  —  Po,exit  which  was  necessary  for  the 
determination  of  the  loss  coefficient.  The  hot  wire  was  also  placed  downstream  of  the  blades  to 
record  the  exit  velocity.  Together  the  differential  pressure  measurement  and  the  exit  velocity  mea¬ 
surement  allowed  easy  determination  of  the  loss  coefficient.  The  baseline  Pak-B  loss  coefficient  is 
documented  in  the  next  chapter  while  the  software  used  in  collection  of  this  data  is  described  further 
in  Appendix  D. 


Figure  22.  Blade  surface  pressure  data  collection  diagram. 


Chapter  5  -  Numerical  Simulation  of  Turbine  Cascade 


5.1  Overview  of  Allison  Software  code 

The  \hne-Blade  Interaction  (VBI)  code  developed  by  Allison  Gas  Turbines  simulates  flow 
passages  in  a  turbine  engine.  Developed  in  1992  for  2-D  flow;  the  code  can  model  a  cascade  of 
turbine  blades  using  a  combination  of  airfoil  wrapping  O-grids  and  passage  filling  H-grids.  While 
the  O-grids  provided  a  high  grid  density  near  the  surface,  the  2-D  grid  generation  program  was 
unable  to  model  the  three  dimensional  discontinuous  surfaces  tested  in  the  tunnel  for  reduction  in 
separation  losses.  The  baseline  Pak-B  shape  was  simulated  to  provide  an  initial  estimation  of  the 
surface  pressure  coefficient  for  the  blade. 

5.2  Grid  Generation  Around  Pak-B  Airfoil 

The  grid  generation  developed  by  Rao  et  al.  [32]  is  a  combination  of  airfoil  conforming  O- 
grids  embedded  within  an  H-grid  that  conforms  to  the  passage  shape.  These  two  different  grids 
overlap  at  their  boundaries  to  exchange  information.  Airfoil  boundary  layer  information  from  the 
O-grid  is  transferred  through  this  chimera  overlap  region  to  the  H-grid.  An  example  of  the  Pak-B 
blade  placed  within  this  O-H  grid  is  presented  in  Figure  23. 

5.3  VBI  Computational  Method 

The  VBI  software  used  an  explicit  5-stage  Runge-Kutta  ordinary  differential  equation  (ODE) 
solution  scheme  for  the  Navier-Stokes  equations.  Rao  et  al.,  [32]  used  a  combination  of  second  and 
fourth  order  artificial  dissipation  terms  with  this  2-D  version.  A  more  detailed  description  of  the 
software  inputs  and  the  calculation  of  turbulent  viscosity  fit  is  provided  in  Appendix  E. 
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Figure  23.  0-H  grid  surrounding  Pak-B  airfoil 


5.4  Two-Dimensional  Software  Results 

VBI  results  include  surface  pressure  contours  and  boundary  layer  velocity  profiles  calculated 
over  the  blade  surface.  The  surface  pressure  distribution  diagrams  are  steady-state  solutions.  The 
boundary  layer  velocity  plots  are  presented  for  the  same  axial  chord  locations  as  the  experimental 
data  presented  in  Chapters  6  through  9. 

A  pressure  distribution  diagram  can  reveal  areas  of  likely  separation  over  the  surface  of  the 
blade.  A  region  of  separated  flow  will  appear  as  a  horizontal  line  with  little  or  no  change  in  C'p. 
Figure  24  shows  the  pressure  distributions  obtained  using  the  VBI  code  for  the  Pak-B  blade  shape. 
As  can  be  seen  in  Figure  24,  the  code  was  unable  to  simulate  the  large  regions  of  separation  found 
on  the  actual  blade  at  Re  =  50,000,  described  in  more  detail  in  Chapter  6. 


The  VBI-2D  pressure  coefficient  results  agree  with  the  experimental  data  very  well  along  the 
pressure  surface,  and  the  first  40%  of  the  suction  surface.  Figure  24  shows  the  ideal  design  pressure 
curve  over  the  aft  60%  of  the  suction  surface.  The  code  did  predict  a  small  separation  as  Reynolds 
number  was  decreased  to  50k.  This  separation  is  also  observable  from  the  boundary  layer  velocity 
plots  presented  in  Figure  25, 


Figure  24.  Two  dimensional  VBI  code  surface  pressure  coefficient  results. 

5.4.1  VBI  2-D  Code  Results  For  Reynolds  Number  =  50k 

Figure  25a  illustrates  the  boundary  layer  velocity  profiles  for  the  simulated  case  of  Reynolds 
number  equal  to  50,000  based  upon  inlet  velocity  and  axial  chord.  The  size  of  the  separation  zone 
shown  in  Figure  25a  is  unusually  thin  due  to  the  lack  of  proper  transitional  modeling  with  the  soft¬ 
ware  code.  The  simulated  flow  transitions  at  a  single  point,  unlike  the  actual  blade  which  undergoes 
transition  over  an  extended  period  after  flow  separation.  The  turbulence  model  used  was  unable  to 
simulate  freestream  turbulence  conditions  which  are  a  significant  factor  in  determining  the  separa¬ 
tion  size  and  location. 
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5.4.2  VBI  2-D  Code  Results  For  Reynolds  Number  =  100k 

The  VBI  2-D  boundary  layer  velocity  profile  results  for  the  100k  Reynolds  number  case  are 
presented  in  Figure  25b.  As  seen  in  Figure  25b,  there  was  no  separation  bubble  present  on  the  suction 
surface  of  the  blade.  Experimentally  collected  boundary  layer  measurements  shown  in  Chapter  6 
show  separation  bubbles  for  ffeestream  turbulence  levels  of  both  1%  and  4%.  The  code  was  unable 
to  model  the  separated  transition  and  re-attachment  that  occurs  on  the  suction  surface  at  Re  =  50k 
and  Re  =  100k  at  a  freestream  turbulence  level  of  1%. 
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a)  simulated  boundary  layer  at  Re  =  50k  and  Tu  =  1% 


b)  simulated  boundary  layer  at  Re  =  100k  and  Tu  =  1% 


Figure  25.  2-D  VBI  code  simulations  of  boundary  layer  velocity. 
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Chapter  6  -  Baseline  LPT  Blade  Performance 


6.1  Baseline  Pak-B  Overview 

The  Pratt  and  Whitney  Pak-B  LPT  blade  shape  was  chosen  for  this  study  for  two  reasons.  The 
first  reason  was  that  it  was  designed  for  use  at  low  Reynolds  numbers  typical  of  U^W  flight  condi¬ 
tions.  The  second  reason  the  Pak-B  blade  shape  was  selected  was  because  the  airfoil  suffers  from 
losses  due  to  flow  separation  at  low  Reynolds  numbers  despite  being  designed  to  operate  under  this 
condition  as  a  highly  aft  loaded  airfoil.  To  illustrate  these  losses  a  linear  cascade  tunnel  was  used 
to  collect  several  types  of  experimental  data  for  the  Pak-B  blade.  Surface  static  pressure  measure¬ 
ments,  boundary  layer  velocities  and  turbulence  profiles  and  loss  coefficients  were  obtained.  The 
results  from  this  investigation  are  presented  in  this  chapter. 

6.2  Experimental  and  Simulated  Surface  Static  Pressures 

To  ensure  proper  simulation  of  the  approximately  2-D  flowfield  around  the  mid-span  surface  of 
the  blade,  surface  static  pressure  measurements  from  blade  #4  and  #6  were  compared  to  each  another 
and  the  VBI  2-D  code  results.  Similar  results  from  blade  #4  and  #6  indicated  a  good  uniformity  in 
flowfield  characteristics  for  the  blade  cassette.  The  pressure  coefficient  curves  shown  in  Figures  26 
through  3 1  were  made  using  the  following  definition: 


Cp  = 


P<0,  inlet  Pi s, 

gullet 


(17) 


where  Po, inlet  and  Uiniet  were  the  total  pressure  and  mean  velocity,  respectively,  measured  upstream 
of  the  blade  cassette.  The  inlet  pressure  and  velocity  were  measured  1  axial  chord  length  upstream 
of  the  test  blades.  The  surface  static  pressure  measurement  for  each  tap  is  denoted  Ps>x  while  p  was 
the  density  based  upon  room  pressure  and  tunnel  temperature.  No  compressibility  corrections  to  the 
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density  were  made  due  to  the  low  Mach  numbers  of  the  flow.  Pressure  curve  data  was  obtained  for 
the  baseline  Pak-B  blade  at  Re  =  50k,  100k,  200k  at  ffeestream  turbulence  levels  of  1%  and  4%. 

As  seen  in  Figure  26,  the  definition  given  in  equation  17  defines  the  pressure  coefficient  values 
as  positive,  with  the  pressure  surface  values  on  the  lower  curve  and  the  suction  surface  values  on 
the  upper  curve.  The  data  for  each  Reynolds  number  and  turbulence  intensity  was  collected  three 
times  and  the  average  values  are  shown  in  Figures  26  through  31.  The  solid  line  present  in  Figures 
26  through  31  are  the  result  from  the  VBI-2D  code.  The  VBI  code  was  explained  in  Chapter  5. 

6.2.1  Cp  Curves  for  Re  =  50k 

The  pressure  tap  blades  have  Cp  distributions  that  agree  with  the  predicted  VBI-2D  code  results 
for  regions  of  attached  flow  on  the  pressure  surface  and  the  initial  40%  of  the  suction  surface.  For 
the  aft  60%  of  the  suction  surface  the  code  fails  to  capture  the  strong  effect  separation  has  on  the 
surface  pressure.  As  seen  in  Figure  26,  the  separation  zone  extends  from  the  62%  axial  chord  line 
back  to  the  trailing  edge  of  the  blade.  The  experimental  data  and  separation  region  agree  with  that 
obtained  previously  by  Murawski  et  al.  [4], 

For  the  increased  ffeestream  turbulence  level  of  4%,  the  laige  separation  bubble  present  in  the 
laminar  flow  case  of  Tu  =  1%  shrinks  to  a  region  extending  from  the  75  to  87%  axial  chord  line  on 
the  suction  surface.  The  VBI-2D  code  is  incapable  of  modeling  the  increased  ffeestream  turbulence 
of  4%,  but  the  previously  shown  VBI  results  are  included  with  the  higher  turbulence  graphs  for 
easier  comparison. 

6.2.2  Cp  Curves  for  Re  =  100k 

The  higher  Reynolds  number  case  of  100k  shows  a  reduction  of  the  size  of  the  separation 
bubble  that  was  present  in  the  Re  =  50k  and  Tu  =  1%  case.  The  size  and  location  of  the  separation 
bubble  for  this  higher  Reynolds  number  case  were  similar  to  those  measured  with  a  higher  turbulence 


50 


Figure  26.  Experimental  and  simulated  Cp  curves  for  Re  =  50k  and  Tu  =  1%. 


Figure  27.  Experimental  and  simulated  Cp  curves  for  Re  =  50k  and  Tu  =  4%. 
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level  of  4%  for  the  Re  =  50k  case.  With  the  increased  energy  in  the  flow,  the  separation  bubble  was 
reduced  to  a  location  starting  at  74%  and  reaching  to  87%  axial  chord. 


Figure  28.  Experimental  and  simulated  Cp  curves  for  Re  =  100k  and  Tu  =  1%. 

The  effect  of  increasing  the  turbulence  level  for  the  100k  case  is  clearly  evident  when  compar¬ 
ing  Figure  29  to  Figure  28.  The  separation  bubble  was  further  reduced  to  a  smaller  range  of  80  to 
86%  axial  chord.  The  trend  in  decreasing  separation  bubble  size  with  increasing  Reynolds  number 
or  turbulence  continues  with  the  experimental  results  of  Re  =  200k. 

6.2.3  Cp  Curves  for  Re  =  200k 

Figure  30  illustrates  the  effect  that  higher  flow  velocities  have  over  the  low  pressure  turbine 
blade.  The  separation  bubble  is  reduced  to  a  region  of  75%  to  83%  axial  chord.  The  size  of  the  mea¬ 
sured  separation  bubble  was  again  dependent  upon  freestream  turbulence  level  as  in  the  previously 
discussed  cases. 

The  flow  condition  that  most  closely  approaches  ideal  fully  attached  flow  over  the  surface  of 
the  blade  is  the  high  Reynolds  number,  high  freestream  turbulence  level  case.  The  combined  effect 
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of  both  ffeestream  turbulence  and  high  flow  eneigy  help  maintain  attached  flow  for  a  longer  distance 
than  in  the  lower  Reynolds  number  and  turbulence  level  cases.  The  effect  the  increase  in  flowfield 
eneigy  is  evident  in  Figure  3 1  where  the  separation  bubble  has  shrunk  to  a  region  extending  from 
78%  to  82%  axial  chord.  The  presence  of  a  separation  bubble  under  these  high  eneigy  conditions 
indicates  there  is  room  for  reduction  in  the  loss  coefficient  for  high  Reynolds  numbers  and  turbulence 
intensities  in  addition  to  the  laige  improvements  possible  at  lower  Reynolds  numbers. 


Figure  3 1 .  Experimental  and  simulated  Cp  curves  for  Re  =  200k  and  Tu  =  4%. 

The  axial  chord  location  of  the  start  of  separation  on  the  suction  surface  was  plotted  in  Figure 
32.  Figure  32  illustrates  the  effect  both  increased  Reynolds  number  and  turbulence  intensity  have 
on  delaying  separation.  To  determine  a  cause  of  the  separation  locations  seen  in  Figures  26  through 
31,  surface  curvature  was  examined.  However,  limitations  in  the  original  coordinate  data  and  dif¬ 
ferencing  scheme  prevented  any  useful  insights  into  a  cause  of  separation  (see  Appendix  G).  The 
Cp  measurements  used  to  make  Figure  32  were  taken  with  differential  pressure  transducers.  A  de¬ 
scription  of  the  measurement  uncertainty  estimation  follows. 
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Figure  32.  Variation  of  separation  start  with  Reynolds  number  and  turbulence  intensity. 
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6.2.4  Uncertainty  in  Cp  Measurements 

The  uncertainty  in  the  Cp  measurements  was  calculated  using  the  method  described  by  Kline 
and  McClintok  [33],  Using  this  method  the  Cp  uncertainties  are  ±0.59,  ±0.18,  and  ±0.04  for 
Reynolds  numbers  of  50k,  100k,  and  200k,  respectively.  The  errors  decrease  with  increasing  Reynolds 
number  because  the  measured  pressures  increase  with  increasing  speed.  The  increased  pressure 
provides  a  larger  range  for  the  pressure  transducers  to  read,  thus  decreasing  the  effect  of  calibrator 
limitations.  The  error  calculations  for  the  Cp  curves  are  presented  in  Appendix  I.  The  velocity  of 
the  boundary  layer  downstream  of  separation  was  measured  directly  using  a  small  hot  film  probe. 

6.3  Boundary  Layer  Velocity  and  Turbulence  Profiles 

The  boundary  layer  velocity  and  turbulence  measurements  were  made  over  the  aft  portion  of 
the  suction  surface  at  67.2%,  73.0%,  79.3%,  84.8%,  and  89.8%  axial  chord.  These  five  locations 
span  the  region  most  effected  by  the  separation  bubble.  Measurements  were  made  at  Re  =  50k 
and  100k  for  turbulence  levels  of  1%  and  4%.  The  boundary  layer  measurements  of  velocity  and 
turbulence  were  made  with  a  TSI  1279S-10A  hot  film  probe.  This  probe  was  moved  by  a  National 
Aperture,  Inc.  MM-3M  5.08  cm  (2  inch)  micro  traverse  controlled  by  aNational  Aperture,  Inc.  MC- 
3SA  servo  amplifier  system.  The  software  for  this  system  is  explained  in  more  detail  in  Appendix 
D 

6.3.1  Boundary  Layer  Measurements  at  Re  =  50k 

Figure  33  illustrates  the  magnitude  of  the  suction  surface  separation  problem  typical  of  high 
altitude,  low  speed  flight.  Defining  the  boundary  layer  thickness  to  be  the  height  at  which  the 
velocity  is  99%  of  the  ffeestream  value,  the  separation  bubble  has  a  maximum  thickness  of  3 1 .0 
mm  (1 .22  inches)  at  the  89.8%  axial  chord  location.  The  x-axis  of  Figure  33  is  the  non-dimensional 
velocity  for  each  of  the  axial  chord  locations  where  the  boundary  layer  measurements  were  taken. 
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Each  velocity  reading  is  nondimensionalized  by  dividing  by  the  freestream  velocity  at  that  chordwise 
location  and  multiplied  by  100.  The  height  of  the  boundary  layer  is  determined  from  the  y-axis 
which  is  in  millimeters. 

The  growth  of  the  suction  surface  boundary  layer  shown  in  Figure  33a  appears  to  be  nonlinear 
with  the  suction  surface.  The  nonlinear  shape  was  due  to  the  convex  geometry  of  the  surface  of  the 
blade,  coupled  with  the  convex  shape  of  the  separation  bubble  itself.  The  combination  of  the  two 
convex  conditions  give  the  nonlinear  appearance  seen  in  Figure  33a  when  the  measured  velocities 
are  plotted  on  an  orthogonal  graph.  The  same  physical  conditions  which  influence  the  velocity 
graphs  also  affect  the  turbulence  graphs. 

In  Figures  33b  through  36b  the  turbulence  intensity  shown  was  calculated  using: 


where  a  is  the  rms  of  the  mean  velocity  U  recorded  at  height  y  above  the  blade  surface  thus  defin¬ 
ing  the  local  turbulence  intensity.  The  definition  provided  in  equation  1 8  was  used  to  provide  di¬ 
rect  comparison  with  the  previously  published  results  of  Murawski  et  al.  [4],  While  the  turbulence 
intensity  within  the  boundary  layer  does  match  the  levels  previously  documented  by  Murawski  et 
al.  [4],  the  high  turbulence  levels  seen  in  the  figures  are  much  smaller  when  calculated  using  the 
freestream  velocity  above  the  boundary  layer  instead  of  the  local  velocity  used  in  equation  18.  For 
example,  a  turbulence  intensity  of  50%  within  the  boundary  layer  is  only  10.3%  when  calculated 
with  the  freestream  velocity.  The  turbulence  intensity  shown  in  Figure  33b  illustrates  an  initially 
laminar  boundary  layer  quickly  separating  with  the  associated  increase  in  turbulence  intensity  to  a 
maximum  of  77.26%  at  the  73%  axial  chord  position.  Downstream  of  73%  axial  chord  the  recircu¬ 
lation  within  the  separation  bubble  is  significant  enough  to  cause  an  almost  steady  turbulence  level 
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of  35-45%.  These  high  turbulence  levels  are  consistent  with  the  previously  published  Pak-B  data 
(Murawski  et  al.  [4]). 

The  effect  of  increased  ffeestream  turbulence  is  clear  when  Figure  34  is  compared  to  Figure  33 . 
The  separation  which  continued  to  grow  to  the  trailing  edge  in  the  laminar  flow  case  was  reduced  to 
a  thin  bubble  with  a  maximum  thickness  of  only  6.38  mm  (0.252  inches)  at  the  89.8%  axial  chord 
line.  The  recirculation  region  ofthe  separation  bubble  was  delayed  from  the  73%  chord  to  the  79.3% 
axial  chord  line. 

Increasing  the  ffeestream  turbulence  level  to  4%  shows  a  delay  in  the  tripping  of  the  boundary 
layer  from  the  73%  chord  position  to  84.8%  axial  chord  line.  The  peak  turbulence  at  84.8%  chord 
is  6 1 .3%  at  a  height  of  1 .059  mm  (0.042  inches).  The  maximum  turbulence  at  the  89.8%  chord  line 
was  only  45.3%  as  compared  to  52.2%  for  the  laminar  flow  case  at  the  same  location. 

6.3.2  Boundary  Layer  Measurements  at  Re  =  100k 

The  higher  energy  in  the  flow  for  the  Re  =  100k  case  reduces  the  separation  present  over 
the  suction  surface  significantly  while  a  recirculation  region  still  exists  at  the  84.8%  chord  line. 
Comparing  the  velocity  and  turbulence  curves  in  Figure  35a  and  b,  respectively,  it  can  be  seen  that 
the  laminar  flow  separates  upstream  of  the  79.3%  axial  chord  line.  Downstream  of  the  84.8%  chord 
line,  the  flow  undergoes  transition  to  turbulent  flow  and  reattaches  upstream  of  die  89.8%  chord 
line. 

The  combined  effects  of  increased  Reynolds  number  and  ffeestream  turbulence  greatly  reduce 
the  separation  over  the  suction  surface  of  the  unmodified  Pak-B  blade.  The  completely  separated 
flow  seen  in  Figure  33  has  been  reduced  to  a  small  separation  bubble  near  the  85%  axial  chord  line. 
The  uncertainty  in  the  velocity  measurements  for  Figures  33  through  36  follows. 
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Figure  33.  Baseline  velocity  and  turbulence  profiles  over  the  suction  surface. 
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Figure  34.  Baseline  velocity  and  turbulence  profiles  with  increased  freestream  turbulence 
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Figure  35.  Baseline  velocity  and  turbulence  profiles  for  Re  =  100k. 
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Figure  36.  Baseline  velocity  and  turbulence  profiles  with  moderate  freestream  turbulence. 
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6.3.3  Uncertainty  in  Boundary  Layer  Measurements 

The  method  used  for  calculating  the  uncertainty  in  the  boundary  layer  hot  film  velocity  mea¬ 
surements  was  described  by  Kline  and  McClintok  [33],  Using  this  method  the  velocity  uncertainties 
are  ±0.24  m/s  and  ±0.15  m/s  for  Reynolds  numbers  of  50k  and  100k,  respectively.  The  errors  de¬ 
crease  with  increasing  Reynolds  number  since  the  measured  pressures  used  in  calibration  of  the  hot 
film  increase  with  increasing  speed  thus  reducing  the  effect  of  calibrator  limitations.  The  boundary 
layer  error  calculations  are  presented  in  Appendix  I  along  with  the  error  analysis  for  the  hot  wire 
measurements  for  the  exit  wake  profiles. 

6.4  Exit  Wlocity  Wake  Profiles 

The  exit  velocity  profile  for  the  Pak-B  blade  was  measured  ^  axial  chord  lengths  downstream 
of  blade  #5.  The  velocity  was  measured  using  a  TSI  1210-T1.5  single  hot  wire  probe  mounted  on 
the  3-axis  Dantec  traverse.  The  velocity  measurements  were  taken  every  0.5  cm  across  the  15.75 
cm  (6.2  inches)  wake.  The  velocity  data  was  then  used  to  calculate  the  momentum  wake  deficit  for 
the  smooth  blade. 

6.4.1  Exit  Wlocity  Wake  Data 

The  measured  exit  velocity  profiles  are  displayed  in  Figures  37  and  38  where  the  pressure  side 
of  blade  #5  is  the  left  portion  of  the  wake  from  60  cm  to  66  cm.  The  suction  side  of  the  wake  extends 
from  66  cm  to  the  right.  Figure  37  contains  the  velocity  wake  data  for  Reynolds  numbers  of  50k, 
100k,  and  200k  for  an  inlet  freestream  turbulence  of  1%.  The  x-axis  of  Figure  37  displays  the  probe 
location  within  the  tunnel  where  the  hot  wire  was  positioned  when  the  data  point  was  collected.  The 
trailing  edge  of  blade  #5  is  approximately  at  66  cm  on  Figure  37.  The  y-axis  of  Figures  37  and  38 
are  the  measured  velocities.  The  bottom  curve  of  Figure  37  was  for  an  inlet  Reynolds  number  of 
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50k.  The  middle  and  upper  curves  of  Figures  37  and  38  were  for  inlet  Reynolds  numbers  of  100k 
and  200k  respectively. 
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Figure  37.  Baseline  Pak-B  exit  velocity  profiles  for  Tu  =  1%. 


Noticeable  in  Figure  37  is  how  the  center  of  the  wake  shifts  to  the  left  with  increasing  Reynolds 
number.  The  shift  in  wake  location  was  due  to  the  decreasing  size  of  the  boundary  layer  with  in¬ 
creased  speed.  As  seen  in  Figures  33  through  36  as  the  Reynolds  number  increases,  the  thickness 
of  the  boundary  layer  decreases.  The  decrease  in  boundary  layer  thickness  leads  to  a  smaller  wake 
closer  to  the  suction  surface. 

Figure  38  contains  exit  wake  velocity  data  for  an  increased  freestream  turbulence  of  4%.  At 
the  higher  turbulence  level  the  boundary  layer  has  more  energy  and  therefore  less  separation.  The 
decrease  in  boundary  layer  thickness  reduces  the  amount  of  wake  shift  compared  to  the  laminar 
case  of  Figure  37.  All  the  velocity  wake  curves  in  Figure  38  have  less  of  a  velocity  deficit  than 
their  counterparts  for  the  laminar  freestream  case  of  Figure  37.  The  velocity  wake  data  was  used  to 
calculate  a  momentum  deficit  value  for  each  flow  condition  tested. 
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Figure  38.  Baseline  Pak-B  exit  velocity  profiles  at  Tu  =  4%. 
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6.4.2  Exit  Momentum  Deficit  Calculations 

Using  the  velocity  information  in  Figures  37  and  38  the  wake  momentum  deficits  were  calcu¬ 
lated.  The  momentum  wake  deficit  was  calculated  using: 


where  Uw  was  the  wake  velocity  measured  at  pitchwise  position  y  within  the  15.75  cm  wake.  The 
freestream  (maximum)  velocity  between  blades  in  the  downstream  wake  is  denoted  as  U^.  Figure 
39  illustrates  the  drop  in  loss  coefficient  with  increasing  Reynolds  number  and  turbulence  level. 
The  difference  between  the  two  turbulence  cases  for  a  Reynolds  number  of  50k  seen  in  Figure  39 
is  due  to  the  strong  effect  fully  separated  flow  has  upon  blade  losses.  The  problem  of  maintaining 

attached  flow  over  the  blade  at  low  Reynolds  numbers  is  the  key  goal  to  this  research. 
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Figure  39.  Baseline  Pak-B  blade  momentum  deficit  measurements. 


6.4.3  Uncertainty  in  Wake  Momentum  Deficit  Measurements 

The  uncertainty  in  the  exit  wake  momentum  deficit  measurements  were  calculated  using  the 
method  described  by  Taylor  [34].  Using  this  method  the  wake  momentum  uncertainties  were  cal- 
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culated  to  be  ±0.030  m,  ±0.009  m,  and  ±0.003  m  for  Reynolds  numbers  of  50k,  100k  and  200k, 
respectively.  The  uncertainties  decrease  with  increasing  Reynolds  number  since  the  measured  pres¬ 
sures  used  in  calibration  of  the  hot  wire  increase  with  increasing  speed,  thus  reducing  the  effect  of 
scale  limitations  from  the  manometer.  Even  with  the  uncertainty  of  the  experimental  data,  the  trend 
of  decreasing  losses  with  increasing  Reynolds  number  is  clearly  evident  from  Figure  39.  This  trend 
can  also  be  observed  from  total  pressure  coefficient  loss  measurements. 

6.5  2-D  Profile  Loss  Coefficient  Measurements 

The  performance  of  the  baseline  LPT  Pak-B  blades  were  also  characterized  by  the  total  pres¬ 
sure  loss  coefficient.  The  loss  coefficient  was  measured  using  a  total  pressure  probe  located  approx¬ 
imately  |  axial  chord  lengths  downstream  of  the  blade  cassette.  The  loss  coefficient  definition  is 
restated  here: 


Po  , inlet  1*0, exit 

7  = - — 2 -  (20) 

2  PP  exit 

where  Po,exit  and  were  averaged  over  the  exit  plane  of  blade  #5  by  dividing  the  sum  by  the 
number  of  samples.  The  inlet  total  pressure  Po, inlet  was  measured  approximately  one  axial  chord 
length  upstream  of  blade  #5.  The  exit  velocity  for  the  experiments  was  small  compared  with  the 
speed  of  sound  so  no  density  correction  for  compressibility  effects  was  used.  The  measurements 
were  made  in  0.5  cm  (0. 197  inches)  increments  across  the  blade  wake  in  the  center  span  of  the  87.63 
cm  (34.5  inches)  high  blade.  Taking  loss  coefficient  measurements  in  the  center  span  provides  the 
2-D  loss  coefficient  which  does  not  account  for  tip  losses,  comer  vortices,  or  3-D  effects. 

6.5.1  Exit  Pressure  Data 

Each  exit  pressure  data  point  was  measured  independently  three  different  times  and  then  en¬ 
semble  averaged.  The  ensemble  averaged  Pq, inlet  -  Po,y  values  for  Re  =  50k,  100k,  200k,  and  for  a 
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freestream  turbulence  level  of  1%  are  in  Figure  40.  The  x-axis  of  Figure  40  is  the  traverse  location 
of  the  probe  as  described  previously  for  Figures  37  and  38.  The  y-axis  of  Figure  40  is  the  difference 
in  total  pressure  between  Po, inlet  and  Po,exit  measured  in  units  of  kg /ms2. 
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pitchwise  traverse  (cm) 

Figure  40.  Baseline  Pak-B  exit  stagnation  pressure  loss  at  Tu  =  1%. 


6.5.2  Baseline  Performance  Loss  Coefficients 

Using  the  pressure  curves  in  Figures  40  and  41  the  loss  coefficients  were  calculated.  The 
resulting  values  for  freestream  turbulence  levels  of  1%  and  4%  are  both  in  Figure  42.  For  the  1% 
turbulence  case  the  loss  coefficient  drops  quickly  with  increasing  Reynolds  number  from  0.155 
at  50k  down  to  only  0.032  at  200k.  For  the  increased  turbulence  of  4%,  the  loss  coefficient  also 
decreases  with  increasing  Reynolds  number,  although  not  as  dramatically  as  for  the  case  of  Tu  = 
1%.  For  increased  freestream  turbulence,  the  loss  coefficient  is  0.058  at  Re  =  50k,  and  drops  to  only 
0.037  at  Re  =  200k.  The  measured  values  of  loss  coefficient  agree  with  estimations  and  curve  fits 
using  previously  recorded  data  from  other  turbine  blades.  The  early  experiments  by  Ainley  [35]  and 
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Ainley  and  Mathieson  [36]  with  the  T.6  blade  shape  provide  comparable  results  to  the  experimental 
data  of  the  Pak-B  blade.  The  curve  fits  calculated  from  Horlock  [37-39]  are  slightly  above  the 
measured  loss  coefficient  for  Tu  =  4%  due  to  the  differences  in  thickness-to-chord  (|) ratio  which 
the  curve  fits  can  not  account  for.  The  relations  of  Craig  and  Cox  [40]  go  further  to  include  the 
effect  of  blade  solidity  (|)  and  trailing  edge  thickness.  The  Craig  and  Cox  [40]  curve  fits  provide 
the  most  accurate  estimation  of  the  loss  coefficient  when  the  flow  is  attached  to  the  blade  at  the 
trailing  edge.  More  information  on  the  estimation  of  loss  coefficient  is  in  Appendix  F. 

6.5.3  Uncertainty  in  Loss  Coefficient  Measurements 

The  experimental  values  represented  in  equation  20  were  taken  using  the  same  instruments  as  in 

equation  17.  Because  of  the  similarity  between  equations  17  and  20,  the  uncertainty  analysis  for  each 

was  very  similar.  Using  the  method  described  in  Appendix  I  for  the  surface  pressure  coefficients, 

the  uncertainty  for  the  loss  coefficients  are  ±0.014,  ±0.003,  and  ±0.0006  for  Re  =  50k,  100k,  and 

200k  respectively. 
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Figure  42.  Baseline  Pak-B  exit  loss  coefficients. 
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Chapter  7  -  Performance  ofPak-B  blade  With  Boundary  Layer 

Trip  at  35%  Chord 

7.1  Trip  Configurations  Tested 

The  losses  associated  with  separated  flow  from  the  Pak-B  blade  can  also  be  controlled  with 
a  simple  boundary  layer  trip.  To  ensure  a  fully  turbulent  boundary  layer  over  the  suction  surface 
of  the  blade,  four  different  methods  were  examined  and  compared.  Each  method  was  tested  to 
determine  if  it  could  produce  an  initially  attached  turbulent  boundary  layer  for  the  region  near  the 
67%  axial  chord  line  where  separation  began  on  the  baseline  blade.  The  methods  tested  include 
surface  roughness  (sandpaper)  applied  to  the  leading  edge  of  the  blade,  roughness  extending  from 
the  leading  edge  to  the  50%  axial  chord  line,  a  small  trip  wire,  and  a  large  trip  wire.  Only  the  large 
trip  wire  produced  a  turbulent  boundary  layer  over  the  aft  portion  of  the  Pak-B  blade. 

The  height  (fccrit)  °f  the  sandpaper  roughness  necessary  to  trip  the  boundary  layer  was  deter¬ 
mined  using  (Schlichting  [41]): 


(21) 


where  ul  was  defined  as  (Schlichting  [41]): 


uh  = 


(22) 


and  r0  was  the  shearing  stress  at  the  blade  surface.  The  value  of  kcrit  determined  from  equations 
21  and  22  for  an  inlet  Reynolds  number  of  50k  was  0.87  mm  (0.0341  inches).  Sandpaper  with  this 
thickness  was  first  applied  to  the  leading  edge  of  the  Pak-B  blade,  extending  to  the  15%  chord  line. 
With  the  leading  edge  sandpaper  in  place,  a  boundary  layer  velocity  profile  was  acquired  with  the 
hot  film  probe  located  at  67%  axial  chord.  As  can  be  seen  in  Figure  43,  the  resulting  velocity  profile 
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was  still  laminar.  In  addition  to  taking  experimental  data,  typical  relations  for  laminar  and  turbulent 
flow  were  also  graphed  on  Figure  43.  The  laminar  relation  used  was  (White  [42]): 


Uy  -  I/q, 


2  y  y 


21 


(23) 


6  8\ 

where  Ua 0  was  the  freestream  velocity  and  8  was  the  boundary  layer  thickness.  The  turbulent  bound¬ 
ary  layer  estimation  used  was  (White  [42]): 


Uy  -  U( 


y 

00  1 8 


1/7 


(24) 


The  sandpaper  was  extended  to  the  50%  axial  chord  location;  however,  this  attempt  was  also 
insufficient  in  producing  a  turbulent  boundary  layer  at  the  67%  chord  line  at  a  Reynolds  number  of 
50k.  Further  attempts  were  made  to  produce  a  turbulent  boundary  layer. 

A  small  trip  wire  was  attached  to  the  suction  surface  of  the  Pak-B  blade.  The  small  wire 
diameter  was  determined  using  (Schlichting  [41]): 


^2*  >  900.  (25) 

v 

Use  of  equation  25  estimated  the  minimum  wire  diameter  to  trip  the  boundary  layer  to  be  1.83  mm 

(0.072  inches).  A  wire  of  this  minimum  diameter  failed  to  trip  the  boundary  layer  into  a  turbulent 

state  when  placed  at  the  50%  axial  chord  line.  A  larger  wire  with  a  diameter  of  2.77  mm  (0.109 

inches)  was  then  used  at  the  10%,  20%,  35%,  and  50%  axial  chord  locations  to  find  the  optimum 

location.  With  the  large  wire  at  the  10%  chord  line,  the  flow  was  tripped  to  a  turbulent  state,  but 

had  partially  re-laminarized  by  the  67%  chord  position  (see  Figures  44  and  45).  The  same  decay  in 

the  turbulence  was  also  measured  with  the  trip  wire  at  the  20%  chord  line  (see  Figures  46  and  47). 

With  the  large  wire  at  the  35%  chord  line,  the  boundary  layer  was  determined  to  be  fully  turbulent  at 

the  67%  chord  line.  However,  the  large  wire  created  a  region  of  separated  flow  when  placed  further 
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Figure  43.  Boundary  layer  trip  attempts  measured  at  67.2%  axial  chord  for  Re  =  50k. 
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Figure  44.  Boundary  layer  profiles  measured  at  67%  chord  with  trip  wire  at  10%  axial  chord. 
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Figure  45.  Boundary  layer  profiles  measured  at  84.8%  chord  with  trip  wire  at  10%  axial  chord. 
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b)  local  turbulence  profile 


Figure  46.  Boundary  layer  profiles  measured  at  67%  chord  with  trip  wire  at  20%  axial  chord. 
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Figure  47.  Boundary  layer  profiles  measured  at  84.8%  chord  with  trip  wire  at  20%  axial  chord. 
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aft  at  the  50%  axial  chord  line  (see  Figure  43).  It  was  determined  that  when  the  flow  was  tripped  to 
an  initially  attached  and  turbulent  state,  the  boundary  layer  remained  attached  to  die  suction  surface 
for  all  Reynolds  numbers  and  freestream  turbulence  intensities  tested. 

7.2  Boundary  Layer  Profiles  Using  Large  Trip  Wire 

Figures  48  through  5 1  show  the  boundary  layer  profiles  over  the  aft  portion  of  the  suction  sur¬ 
face  when  a  2.77  mm  (0.109  inches)  diameter  trip  wire  was  placed  at  the  35%  axial  chord  location. 
The  flow  remains  attached  over  the  entire  suction  surface  even  for  a  low  Reynolds  number  of  50k. 
Noteworthy  in  Figures  48  through  5 1  was  the  gradual  re-laminarization  of  the  boundary  layer.  This 
was  evident  from  the  initially  turbulent  velocity  profile  at  67.2%  chord  line  to  the  transitional  ve¬ 
locity  profile  at  89.8%  chord.  The  transitional,  attached  boundary  layer  at  the  exit  from  the  blade 
resulted  in  a  smaller  velocity  wake  profile  downstream  for  a  Reynolds  number  of  50k. 

7.3  Exit  Velocity  Wake  Measurements  Using  Large  Trip  Wire 

The  exit  velocity  profile  was  measured  using  the  TSI  1 2 1 0-T 1 . 5  hot  wire  approximately  |  ax¬ 
ial  chord  lengths  downstream  from  the  blade  with  the  wire  trip.  As  seen  in  Figure  52,  the  tripped 
velocity  profile  was  smaller  than  the  baseline  blade  for  a  Reynolds  number  of  50k  and  a  freestream 
turbulence  level  of  1%.  At  all  other  Reynolds  numbers  and  turbulence  intensities  the  velocity  wakes 
were  larger  for  the  tripped  blade  than  the  baseline  blade.  This  was  due  to  the  fact  that  the  unmodi¬ 
fied  boundary  layer  did  not  have  large  regions  of  flow  separation  at  higher  Reynolds  numbers  and 
turbulence  intensities.  The  trip  wire  created  a  turbulent  and  transitional  boundary  layer  which  in¬ 
creased  the  losses  when  the  flow  would  otherwise  be  attached  and  laminar. 

7.4  Loss  Coefficient  Measurements  Using  Large  Trip  Wire 

The  trend  in  increased  losses  for  Reynolds  numbers  of  1 00k  and  200k  was  also  evident  with  the 
exit  pressure  data  (see  Figure  53).  The  exit  total  pressure  wakes  were  measured  at  the  same  down- 
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Figure  48.  Measured  velocity  and  turbulence  profiles  for  baseline  blade  with  trip  at  35%. 
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Figure  49.  Measured  velocity  and  turbulence  profiles  for  baseline  blade  with  trip  at  35%. 
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a)  boundary  layer  velocity  measurements  at  Re  =  50k  and  Tu  =  4% 


Turbulence  Intensity 


b)  boundary  layer  turblence  measurements  at  Re  =  50k  and  Tu  =  4% 


Figure  50.  Measured  velocity  and  turbulence  profiles  for  baseline  blade  with  trip  at  35%  with  Tu  = 
4%. 
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Velocity  (U/UmaxMOO) 

a)  boundary  layer  velocity  measurements  at  Re  =  100k  and  Tu  =  4% 


Turbulence  Intensity 


b)  boundary  layer  turblence  measurements  at  Re  =  100k  and  Tu  =  4% 


Figure  5 1 .  Measured  velocity  and  turbulence  profiles  for  baseline  blade  with  trip  at  35%  with  Tu 
4%. 
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stream  locations  as  the  exit  velocities.  Only  the  case  of  a  Reynolds  number  of  50k  and  freestream 
turbulence  level  of  1%  showed  a  reduction  in  the  losses  through  its  prevention  of  the  laige  separa¬ 
tion  bubble  that  occurs  for  the  baseline  blade.  Figure  54  illustrates  the  increase  in  losses  generated 
using  the  trip  wire  at  higher  Reynolds  numbers  and  turbulence  intensities. 
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a)  Exit  pressure  data  for  baseline  blade  with  trip  wire  and  Tu  =  1% 
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b)  Exit  pressure  data  for  baseline  blade  with  trip  wire  and  Tu  =  4% 


Figure  53.  Stagnation  pressure  loss  for  baseline  blade  with  trip  wire  at  35%  axial  chord. 
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a)  Wake  momentum  values  for  baseline  blade  and  with  trip  wire 
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b)  Loss  coefficient  for  baseline  blade  and  with  trip  wire 


Figure  54.  Wake  velocity  and  total  pressure  loss  coefficients  for  baseline  blade  with  trip  wire  at  35%. 
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Chapter  8  -  Loss  Reduction  Using  Dimples 


8.1  Background  and  Separation  Control  Methodology 

The  existing  separation  problem  on  the  low  pressure  turbine  blades  could  have  been  solved  in 
a  variety  of  ways.  As  seen  in  Figure  55  two  basic  possibilities  exist.  One  possibility  was  to  redesign 
the  blades  which  allowed  the  engineer  to  retain  nothing  of  the  original  shape  if  that  was  found  to  be 
desirable .  The  other  solution  was  to  modify  the  existing  blades  in  some  way  that  simply  corrects  the 
separation  problem  while  leaving  the  other  aerodynamic  characteristics  unchanged.  Modifying  the 
existing  blades  was  the  quickest  and  most  cost  effective  solution,  and  was  examined  in  this  research. 


Figure  55.  Diagram  of  considered  separation  control  methods. 


The  modifications  can  be  further  categorized  into  two  major  divisions,  those  requiring  auxiliary 
power  and  those  that  do  not  (Chang  [43]).  Johnston  and  Nishi  [44]  have  described  methods  requiring 
auxiliary  power  such  as  blowing  jets,  boundary  layer  suction,  and  actuated  vortex  generators.  These 
methods  were  not  considered  for  this  project  due  to  their  increased  complexity,  and  the  additional 
engine  modifications  that  would  be  necessary  to  implement  them. 
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Considering  these  possibilities,  only  passive  surface  modifications  were  examined  for  the  re¬ 
duction  of  separation  losses.  The  simple  implementation  of  passive  techniques  provided  the  goal 
of  maximum  efficiency  and  lowest  implementation  cost  and  time  required  for  installation.  This 
research  provided  a  simple  and  economical  solution  to  the  problem  of  airflow  separation  on  low- 
pressure  turbine  blades  that  can  be  easily  adapted  to  a  variety  of  engines  and  conditions,  while  not 
impeding  performance  that  was  otherwise  nominal. 

An  extensively  tested  economical  means  for  reducing  separation  is  the  use  of  vortex  gener¬ 
ators.  Many  vortex  generators  are  typically  small  winglets  attached  to  the  problem  surface  (flat 
plate,  diffuser,  wing,  etc.).  The  winglet  devices  eneigize  the  slower  fluid  within  the  boundary  layer 
by  moving  higher  speed  fluid  from  the  upper  boundary  layer  or  freestream  down  into  the  slower 
boundary  layer  (Taylor  [45]).  This  mixing  is  aided  by  the  creation  of  trailing  vortices  which  ex¬ 
tend  downstream  into  the  boundary  layer.  The  efficiency  and  effect  of  this  vortex  is  dependent  on 
the  flow  conditions,  vortex  generator  design,  and  most  importantly  the  height  of  the  vortex  gen¬ 
erator.  Lin  et  al.  [18]  have  shown  that  submerged  vortex  generators  with  overall  height  less  than 
the  boundary  layer  thickness  have  less  drag  than  those  which  extend  fully  into  the  freestream,  but 
accordingly  have  a  shorter  distance  of  effectiveness  downstream.  Halstead  and  Solomon  [24]  ex¬ 
pressed  concerns  that  the  increased  height  of  vortex  generators  which  extend  into  the  freestream 
can  make  these  devices  difficult  to  manufacture  and  clean  once  in  use.  Rejected  passive  methods 
initially  considered  such  as  riblets,  ramps,  and  vanes  are  discussed  in  Appendix  J.  These  difficul¬ 
ties  have  indicated  that  the  best  options  for  modifying  the  Pak-B  blades  were  submerged  design 
concepts  that  use  vortex  generation  such  as  engraved  dimples. 
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8.2  Previous  Studies  Using  Recessed  Dimples 

The  vortex  generating  dimple  surface  treatments  used  in  this  research  were  based  upon  the  golf- 
ball  dimple  structure  previously  tested  by  Bearman  and  Harvey  [16, 17]  and  Kimura  and  Tsutahara 
[46],  Bearman  and  Harvey  [16]  first  tested  the  recessed  dimple  structures  on  a  sphere  to  simulate 
the  effects  on  a  golf  ball.  Comparison  between  a  smooth  sphere  and  the  dimpled  sphere  showed  a 
sustained  drag  reduction  for  6  x  105  <  Rep  <  3  x  106  where  R gd  is  the  Reynolds  number  based 
upon  sphere  diameter.  Bearman  and  Harvey  [  16, 17]  cut  the  dimples  to  a  specific  depth  using  a  non- 
dimensional  roughness  ratio  of  =  9  x  10-3  where  fc  is  the  maximum  depth  of  the  impression 
and  D  is  the  sphere  diameter.  Additional  tests  by  Bearman  and  Harvey  [16, 17]  using  hexagonal 
dimples  instead  of  circular  ones  showed  an  increase  in  Cl  by  20%  and  a  reduction  in  Co  by  8%. 
The  hexagonal  dimples  have  six  sharp  edges  with  which  to  generate  downstream  vortices  which  aid 
in  the  reduction  of  separation  more  efficiently  than  standard  circular  dimples.  Dimple  results  were 
later  observed  and  confirmed  to  reduce  the  size  of  the  separation  zone  by  Nakayama  et  al.  [47]  using 
a  test  sphere. 

Bearman  and  Harvey  [  1 7]  also  tested  the  dimples  on  a  cylinder  to  reduce  drag  at  low  Reynolds 
numbers.  The  dimples  were  cut  completely  around  the  cylinder,  each  row  proceeding  spanwise 
along  its  length  with  subsequent  rows  staggered  to  provide  a  honeycomb  appearance.  Each  dimple 
on  the  cylinder  was  cut  using  ^  =  9x  10-3  as  in  the  previous  spherical  experiments.  Results 
from  the  dimpled  cylinder  showed  a  significant  drag  reduction,  however,  no  attempts  were  made  to 
optimize  the  dimple  size,  location,  or  pattern. 

8.3  Dimple  Variations  Previously  Tested 

Kimura  and  Tsutahara  [46]  tested  the  effects  of  cutting  a  channel  along  the  span  of  a  cylinder 
to  simulate  a  2-D  golf  ball  dimple.  Three  cylinders  were  tested  using  three  different  channels.  The 
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channels  were  cut  to  different  depths  using  three  different  ^  values,  1.2  x  10“2, 1.7  x  10~2,  and  7 
x  10-3.  The  tests  were  performed  at  ReD  =  2200  and  a  freestream  turbulence  of  2%.  Each  cylinder 
was  rotated  to  study  the  effects  of  placing  the  single  channel  at  different  downstream  locations  from 
the  initial  stagnation  point.  Results  indicated  that  separation  was  reduced  over  the  surface  of  the 
cylinder  when  the  groove  was  place  at  83°  from  the  stagnation  line.  Separation  occurred  at  92° 
for  the  smooth  cylinder,  and  was  delayed  only  to  94.5°  for  the  best  groove  tests.  The  spanwise 
grooves  were  found  to  be  most  effective  when  recirculation  regions  formed  within  the  channel. 
The  2-D  nature  of  the  surface  modification  minimized  the  creation  of  downstream  vortices,  making 
the  2-D  channel  technique  less  effective  in  separation  control  than  the  3-D  dimples  chosen  for  this 
experiment. 

8.4  Experiment  Test  Article 

Using  a  scale  based  upon  boundary  layer  thickness,  the  dimple  structures  from  the  previous 
experiments  of  Bearman  and  Harvey  [16, 17]  were  cut  into  the  suction  surface  of  the  test  blade. 
Each  dimple  was  cut  to  a  depth  of  1.588  mm  (l/16th  of  an  inch)  using  a  5.08  cm  (2  inch)  diameter 
ball  end  mill.  Since  the  surface  of  the  blade  is  curved,  the  resulting  impressions  are  elliptical  with 
a  streamwise  length  of  15. 13  mm  (0.60  inches)  and  a  spanwise  length  of  17.53  mm  (0.69  inches). 
The  dimples  are  2.22  cm  (0.875  inches)  apart  center-to-center.  The  resulting  blade  is  in  Figure  56. 
The  axial  chord  length  of  the  blade  used  in  the  tests  is  seven  times  engine  scale,  so  the  resulting 
impressions  on  a  operational  engine  blade  would  be  only  2.16  mm  (0.086  inches)  x  2.50  mm  (0.099 
inches)  which  could  be  easily  molded  into  the  surface  of  the  blade  during  initial  manufacturing. 

As  seen  in  Figure  56,  there  are  three  sets  of  dimples  engraved  into  the  blade  at  the  65%,  55%, 
and  50%  axial  chord  lines  The  combination  of  three  locations  on  a  single  blade  allowed  for  multiple 
data  collections  to  be  performed  with  a  single  test  article.  The  benefits  of  having  dimples  at  different 
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upstream  locations  from  the  separation  bubble  was  measured  using  a  single  test  blade.  Each  set  of 
dimples  were  cut  to  the  same  size  and  depth.  The  standardization  of  dimple  size  allowed  for  easier 
comparison  of  boundary  layer  measurements  between  the  different  sets  of  dimples. 


Figure  56.  Dimples  milled  into  the  suction  surface  of  the  blade  at  three  different  locations. 


8.5  Boundary  Layer  \felocity  and  Turbulence  Measurements  Downstream  of 
Dimples 

Boundary  layer  measurements  of  velocity  and  turbulence  were  made  with  the  same  TSI 1279S- 

10A  hot  film  probe  used  previously  with  the  baseline  blade  described  in  Section  6.3.  Using  the  hot 

film  probe,  velocity  and  turbulence  measurements  were  made  over  the  aft  portion  of  the  modified 

suction  surface  at  five  different  axial  chord  locations  at  two  different  spanwise  locations.  Measure- 
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ments  were  made  both  directly  downstream  of  a  dimple  and  between  adjacent  dimples  (see  Figure 
57).  These  ten  locations  span  the  region  most  effected  by  the  dimples  cut  into  the  surface.  The 
software  for  this  system  is  explained  in  more  detail  in  Appendix  D. 

X  X  X  X  X 

X  X  X  X  X 

5  downstream  measurement  locations 


Figure  57.  Detail  of  ten  downstream  measurement  locations  for  dimples  (not  to  scale). 

8.5.1  Boundary  Layer  Measurements  at  Re  =  50k  for  Dimples  at  50%,  55%  and  65%  axial 
chord 

The  data  illustrated  in  Figures  58  and  60  present  the  dramatic  difference  in  boundary  layer 
thickness  reduction  for  the  dimple  surface  modifications.  The  boundary  layer  now  remains  attached 
with  only  a  small  separation  zone  starting  at  73.0%  and  extending  to  the  79.3%  chord  line.  Down¬ 
stream  of  the  79.3%  chord  line  the  flow  re-attaches.  Additional  figures  are  presented  in  Appendix 
A  which  show  very  little  difference  in  the  velocity  profiles  behind  a  dimple  and  between  dimples. 
The  turbulence  levels  downstream  of  the  dimples  shows  the  progression  and  dissipation  of  the  tur¬ 
bulence  near  the  surface.  The  effect  of  increased  ffeestream  turbulence  can  be  seen  in  Appendix  A 
(Figures  102  and  103)  where  a  further  reduction  in  boundary  layer  size  was  evident.  The  difference 
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between  the  baseline  and  dimpled  blade  velocity  and  turbulence  profiles  were  also  reduced  with  the 
greater  ffeestream  turbulence. 
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Height  (mm)  Height  (mm) 


Velocity  (U/UmaxMOO) 

a)  boundary  layer  velocity  measurements  at  Re  =  50k  and  Tu  =  1% 


Turbulence  Intensity 

b)  boundary  layer  turbulence  measurements  at  Re  =  50k  and  Tu  =  1% 


Figure  58.  Measured  velocity  and  turbulence  profiles  behind  a  dimple  at  50%. 
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Height  (mm)  Height  (mm) 


Velocity  (U/Umax  *  100) 

a)  boundary  layer  measurements  at  Re  =  50k  and  Tu  =  1% 


Turbulence  Intensity 

b)  boundary  layer  measurements  at  Re  =  50k  and  Tu  =  1% 


Figure  59.  Measured  velocity  and  turbulence  profiles  behind  a  dimple  at  55%. 
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Height  (mm)  Height  (mm) 


Velocity  (U/Umax  *  100) 

a)  boundary  layer  velocity  measurements  at  Re  =  50k  and  Tu  =  1% 


Turbulence  Intensity 

b)  boundary  layer  turbulence  measurements  at  Re  =  50k  and  Tu  =  1% 


Figure  60.  Measured  velocity  and  turbulence  profiles  behind  a  dimple  at  65%. 
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8.5.2  Boundary  Layer  Measurements  at  Re  =  100k  for  Dimples  at  50%,  55%  and  65%  axial 
chord 

The  data  in  Figures  61  through  63  illustrates  the  diminishing  effect  dimples  have  at  higher 
Reynolds  numbers  and  increased  freestream  turbulence.  The  profiles  still  show  an  improvement 
over  the  baseline  blade  performance;  however,  the  effect  was  not  as  strong  as  at  the  lower  Reynolds 
numbers  and  freestream  turbulence.  While  the  turbulence  profiles  were  nearly  identical  between 
dimples  and  downstream  of  a  dimple,  the  velocity  profile  at  79.3%  chord  still  shows  a  small  sepa¬ 
ration  zone  between  the  dimples  (see  Figure  110  in  Appendix  A).  The  effect  of  the  dimples  at  all 
three  locations  was  fully  mixed  downstream  of  the  84.8%  chord  line,  which  allowed  measurement 
of  the  exit  velocity  wake  at  a  single  spanwise  location. 
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Height  (mm) 


Height  (mm)  Height  (mm) 


Velocity  (U/Umax*  100) 

a)  boundaiy  layer  at  Re  =  100k  and  Tu  =  1% 


Figure  62.  Measured  velocity  and  turbulence  profiles  behind  a  dimple  at  55%. 
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Height  (mm)  Height  (mm) 


Velocity  (U/UmaxMOO) 

a)  boundary  layer  velocity  measurements  at  Re  =  100k  and  Tu  =  1% 


b)  boundaiy  layer  turbulence  measurements  at  Re  =  100k  and  Tu  =  1% 


Figure  63.  Measured  velocity  and  turbulence  profiles  behind  a  dimple  at  65%. 
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8.6  Exit  Velocity  Wake  Profiles  For  Dimpled  Blade 


The  exit  velocity  wake  data  was  collected  using  the  same  instruments  as  described  in  Section 
6.4.  The  Dantec  traverse  moved  the  hot  wire  probe  across  the  wake  of  the  dimpled  blade  at  a  location 
^  chord  length  downstream.  The  velocity  information  collected  every  0.5  cm  (0.197  inches)  was 
then  used  to  calculate  the  wake  momentum  deficit. 

Figures  64  through  66  show  the  exit  velocity  wakes  for  the  dimples  located  at  50%,  55%,  and 
65%  axial  chord  and  for  freestream  turbulence  levels  of  1%  and  4%.  The  inlet  Reynolds  numbers  of 
50k,  100k,  and  200k  are  presented  as  three  horizontal  bands  at  approximately  6  m/s,  13  m/s  and  25 
m/s,  respectively.  The  difference  between  the  graphs  for  Tu  =  1%  and  Tu  =  4%  illustrate  the  effect 
increased  turbulence  has  on  the  size  of  the  blade  wake.  The  increased  turbulence  energizes  the  flow 
at  low  Reynolds  numbers  which  reduces  the  separation  considerably. 

As  with  the  baseline  blade,  as  Reynolds  number  and  turbulence  level  increases,  the  size  of  the 
separation  zone  on  the  suction  side  of  the  blade  decreases.  As  the  separation  zone  decreases  in  size, 
the  wake  diminishes.  The  suction  side  of  the  blade  affected  the  flow  from  66cm  to  the  right  side  of 
Figure  64a  and  b.  This  shrinking  of  the  boundary  layer  thickness  with  increased  Reynolds  number 
explains  why  the  velocity  wake  shifts  to  the  left  with  increased  velocity. 

8.6.1  Exit  Momentum  Deficit  Calculations 

The  calculated  momentum  deficits  for  the  blade  modified  with  dimples  are  shown  in  Figure 
67.  The  greatest  impact  the  dimples  had  was  at  Re  =  50k,  where  even  at  higher  turbulence  levels 
of  4%,  the  wake  was  significantly  reduced  over  the  baseline  blade.  At  higher  Reynolds  numbers 
the  beneficial  influence  of  the  dimples  decreases.  The  trend  of  decreasing  benefit  with  increasing 
Reynolds  number  and  freestream  turbulence  was  also  evident  with  the  measured  loss  coefficients. 
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pitchwise  traverse  (cm) 

a)  Exit  velocity  profiles  for  Tu  =  1% 


pitchwise  traverse  (cm) 


b)  Exit  velocity  profiles  for  Tu  =  4% 

Figure  64.  Exit  velocity  profiles  for  Dimples  at  50%  chord. 


Velocity  (m/s)  Velocity  (m/s) 


pitchwise  traverse  (cm) 

a)  Exit  velocity  profiles  for  Tu  =  1% 


pitchwise  traverse  (cm) 

b)  Exit  velocity  profiles  for  Tu  =  4% 


Figure  65.  Exit  velocity  profiles  for  Dimples  at  55%  chord. 
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Velocity  (m/s)  Velocity  (m/s) 


a)  Baseline  Pak-B  blade 


b)  Dimples  at  50%  chord 


Reynolds  Number  (thousands)  Reynolds  Number  „housands) 

c)  Dimples  at  55%  chord  d)  Dimples  at  65%  chord • 

Figure  67.  Wake  momentum  values  for  baseline  blade,  and  dimples  at  50%,  55%,  and  65%  chord 


105 


8.7  Loss  Coefficient  Reduction  Using  Dimples 

The  loss  coefficient  pressure  data  was  collected  using  the  same  total  pressure  probe  as  described 
in  Section  6.5.  The  total  pressure  wake  data  was  collected  \  axial  chord  lengths  downstream  of 
the  dimpled  blade.  The  probe  was  moved  across  the  wake  of  the  blade  in  0.5  cm  (0.197  inches) 
increments  using  the  Dantec  3-axis  traverse  and  the  Pitot  Traverse  Data  Acquisition  VI.  The  software 
is  explained  in  more  detail  in  Appendix  D.  The  software  saved  to  a  file  the  exit  total  pressures, 
velocity,  turbulence  intensity  and  location  of  the  wake  measurements. 

8.7.1  Exit  Pressure  Data  For  Tu  =  1% 

The  exit  pressures  are  similar  in  form  and  magnitude  to  the  baseline  blades,  but  with  reductions 
in  width  and  depth.  These  reductions  give  the  modified  blade  improved  efficiency.  As  with  the  exit 
velocity  profiles  as  the  Reynolds  number  increases  the  width  of  the  blade  wake  decreases,  and  so  the 
wakes  appear  to  shift  to  the  left.  Since  the  separation  bubble  is  much  smaller  on  the  dimpled  blade 
surface,  the  wake  shift  effect  for  the  dimpled  blade  was  less  than  the  smooth  blade  tested  under  the 
same  circumstances.  This  shifting  was  also  reduced  with  increased  freestream  turbulence. 
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Figure  69.  Exit  stagnation  pressure  loss  for  dimples  at  55%. 
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Figure  70.  Exit  stagnation  pressure  loss  for  dimples  at  65%. 
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8.7.2  Exit  Pressure  Data  for  Tu  =  4% 


The  increased  ffeestream  turbulence  eneigized  the  boundary  layer  sufficiently  to  prevent  large 
separation  bubbles  forming  on  the  suction  surface.  The  lack  of  separation  at  all  Reynolds  numbers 
indicates  there  should  not  be  significant  variation  in  the  centerline  of  the  downstream  wake  of  the 
blade.  The  expected  alignment  of  the  blade  wakes  is  shown  in  Figures  71  through  73.  The  exit 
pressures  were  also  used  to  determine  the  loss  coefficients  for  the  different  dimple  studies. 
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Figure  71.  Exit  stagnation  pressure  loss  for  dimples  at  50%  with  Tu  =  4%. 


8.7.3  Loss  Coefficient  for  Dimpled  Blade 

The  loss  coefficients  for  the  dimpled  blade  are  illustrated  in  Figure  74.  Comparing  Figure  74 
to  the  baseline  blade  loss  coefficients  in  Section  6.5  shows  how  much  the  dimples  help.  The  loss 
coefficient  at  Re  =  50k  and  Tu  =  1%  was  reduced  by  19%  to  42%  with  the  dimples.  The  coefficient 
was  also  reduced  at  Re  =  100k  and  Re  =  200k.  At  Re  =  200k,  the  loss  coefficients  for  the  baseline 
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Figure  72.  Exit  stagnation  pressure  loss  for  dimples  at  55%  with  Tu  =  4%. 
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Figure  73.  Exit  stagnation  pressure  loss  for  dimples  at  65%  with  Tu  =  4%. 
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blade  were  slightly  more  than  the  dimpled  blade,  illustrating  the  non-damaging  nature  of  the  dimples 
at  higher  Reynolds  numbers. 
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Figure  74.  Loss  coefficient  for  baseline  blade  and  dimples  located  at  50%,  55%,  and  65%  axial 
chord. 
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Chapter  9  -  Loss  Reduction  Using  V-Grooves 


9.1  Background  and  Previous  Studies 

The  V-groove  concept  is  very  similar  to  riblets.  However,  the  V-grooves  are  much  larger  in 
size  than  the  riblet  concept  and  were  engraved  into  the  surface  with  regions  of  unmodified  surface 
separating  each  groove  (see  Figure  75).  The  V-grooves  tested  were  machined  to  a  depth  of  h  = 
0.769(50,  a  width  of  w  =  1.538#o,  and  a  spacing  of  s  =  3.077<5o  where  <5o  is  the  original  unmodified 
thickness  of  the  boundary  layer.  Tested  by  Lin  et  al.  [18]  these  V-grooves  delayed  the  onset  of 
separation  by  1.27  cm  (0.5  inches)  and  reduced  the  overall  reattachment  distance  by  up  to  66% 
using  a  convex  curved  ramp  with  no  pressure  gradient  and  a  tripped  boundary  layer.  Recommended 
placement  of  the  grooves  was  from  just  upstream  of  the  baseline  separation  point  to  the  trailing  edge 
of  the  surface  when  there  exits  a  turbulent  boundary  layer  that  is  separated. 


h 

front 


Figure  75.  Large  scale  V-groove  separation  control  design. 
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9.2  V-Groove  Experiment  Test  Article 

A  smooth  Pak-B  blade  was  modified  to  include  engraved  V-grooves  over  the  aft  portion  of  the 
suction  surface.  Two  variants  of  the  V-grooves  were  tested.  The  first  set  of  V-grooves  started  at  the 
55%  axial  chord  location  and  extended  back  to  the  95%  chord.  The  second  V-groove  application 
used  the  grooves  starting  at  the  60%  chord  location.  Testing  two  different  sets  of  V-grooves  allowed 
the  determination  start  location  had  upon  separation  loss  reduction.  The  size  of  the  V-grooves  was 
determined  from  using  the  boundary  layer  thickness  of  the  previous  experiment  of  Lin  et  al.  [18] 
and  scaling  this  to  the  Pak-B.  Each  engraved  groove  was  cut  using  a  6.35  mm  (0.25  inch)  diameter 
V-shaped  router  bit.  The  resulting  grooves  were  2.51  mm  (0.099  inches)  deep,  5.03  mm  (0.198 
inches)  wide  and  spaced  10.0  mm  (0.395  inches)  apart.  Each  set  of  V-grooves  spanned  19.05  cm 
(7.5  inches)  of  the  test  blade  z-axis.  Measurements  of  the  exit  total  pressure,  exit  velocity  and 
boundary  layer  were  completed  as  in  the  previous  case  of  the  dimpled  blade. 

9.3  Boundary  Layer  Velocity  and  Turbulence  Measurements  with  V-Grooves 
at  55%  and  60%  Chord  for  Re  =  50k 

Figures  76  and  77  show  the  velocity  and  turbulence  profiles  for  the  blade  with  V-grooves 
engraved  in  the  suction  surface  at  an  inlet  Reynolds  number  of  50k.  Boundary  layer  velocity  and 
local  turbulence  intensity  measurements  were  taken  at  two  different  spanwise  locations  for  each  of 
the  five  axial  chord  positions.  One  set  of  boundary  layer  profiles  were  taken  down  to  the  surface 
separating  the  V-grooves.  The  second  set  of  data  was  collected  over  a  V-groove  with  the  probe 
extending  0.9  mm  (0.035  inches)  into  the  V-groove  itself  (see  Appendix  B).  As  can  be  seen  by 
comparing  the  boundary  layer  in  Figure  76  and  116  (Appendix  B),  the  V-groove  has  minimal  effect 
on  the  velocity  and  turbulence  levels  above  the  surface  of  the  blade.  Boundary  layer  measurements 
were  also  taken  at  a  Reynolds  number  of  100k. 


114 


Height  (mm)  Height  (mm) 


Velocity  (U/Umax  *  100) 

a)  boundary  layer  velocity  measurements  at  Re  =  50k  and  Tu  =  1% 


Turbulence  Intensity 

b)  boundary  layer  turbulence  measurements  at  Re  =  50k  and  Tu  =  1% 


Figure  76.  Measured  velocity  and  turbulence  profiles  between  V-grooves  starting  at  55%. 
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a)  boundary  layer  velocity  measurements  at  Re  =  50k  and  Tu  =  1% 


Turbulence  Intensity 


b)  boundary  layer  turbulence  measurements  at  Re  =  50k  and  Tu  =  1% 


Figure  77.  Measured  velocity  and  turbulence  profiles  between  V-grooves  starting  at  60%. 
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9.4  Boundary  Layer  Velocity  and  Turbulence  Measurements  with  V-Grooves 
at  55%  and  60%  Chord  for  Re  =  100k 

The  boundary  layer  velocity  profiles  in  Figures  78  and  79  illustrate  attached  boundary  layer 
present  using  the  V-grooves  with  an  initially  laminar  boundary  layer.  The  attached  flow  measured 
at  Re  =  1 00k  was  very  different  from  the  separated  boundary  layer  measured  at  Re  =  50k  for  the  V- 
groove  and  baseline  blades.  Figures  78  and  79  also  highlight  how  effective  the  V-grooves  were  at 
eliminating  boundary  layer  separation  at  a  Reynolds  number  of  1 00k  compared  to  the  baseline  blade . 
The  baseline  and  V-groove  blades  were  also  tested  with  a  boundary  layer  trip  to  better  emulate  the 
turbulent  boundary  layer  conditions  of  the  previous  experiment  by  Lin  et  al.  [18], 

9.5  liirbulent  Boundary  Layer  Measurements  with  V-Grooves  at  55%  and 
60%  Chord 

A  trip  wire  was  placed  upstream  of  the  V-grooves  at  the  35%  axial  chord  line  to  create  a 
turbulent  boundary  layer.  The  wire  was  2.77  mm  (0.109  inches)  in  diameter  and  spanned  across  the 
V-groove  region  of  the  blade.  To  illustrate  the  turbulent  boundary  layer  created  over  the  V-grooves, 
approximations  for  a  laminar  and  turbulent  boundary  layer(White  [42])  are  included  in  Figures  80 
through  83. 

9.6  Exit  Velocity  Wake  Profiles  for  V-Grooves 

The  exit  velocity  wake  data  was  collected  for  the  V-grooves  using  the  same  technique  described 
in  Section  6.4.  The  hot  wire  probe  was  inserted  through  traverse  slot  #3  into  the  wake  of  the  V- 
groove  blade  approximately  |  axial  chord  lengths  downstream.  Figures  84  and  85  show  the  wake 
velocity  for  the  V-groove  blade  with  turbulence  levels  of  1%  and  4%.  In  addition  to  the  initially 
laminar  boundary  layer  results  shown  in  Figures  84  and  85,  experimental  data  was  also  collected 
for  the  V-groove  blade  with  a  boundary  layer  trip  wire. 
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Height  (mm)  Height  (mm) 


Velocity  (U/Umax  *  100) 

a)  boundary  layer  velocity  measurements  at  Re  =  100k  and  Tu  =  1% 


b)  boundary  layer  turbulence  measurements  at  Re  =  100k  and  Tu  =  1% 


Figure  78.  Measured  velocity  and  turbulence  profiles  between  V-grooves  starting  at  55%. 
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Velocity  (U/Umax  *  100) 


a)  boundary  layer  velocity  measurements  at  Re  =  100k  and  Tu  =  1% 


Turbulence  Intensity 

b)  boundary  layer  turbulence  measurements  at  Re  -  100k  and  Tu  =  1% 


Figure  79.  Measured  velocity  and  turbulence  profiles  between  V-grooves  starting  at  60%. 
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-laminar  approx. 

+  turbulent  approx. 
a  Re  =  50k,  Tu  =  1% 
o  Re  =  100k,  Tu  =  1% 
a  Re  =  50k,  Tu  =  4% 
♦  Re  =  100k,  Tu  =  4% 


Figure  80.  Tripped  boundary  layer  profiles  at  67%  chord  for  V-grooves  starting  at  55%. 
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Figure  8 1 .  Tripped  boundary  layer  profiles  at  84.8%  chord  for  V-grooves  starting  at  55%. 
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82.  Tripped  boundary  layer  profiles  at  67%  chord  for  V-grooves  starting  at  60%. 


u/Umax*100 


Figure  83.  Tripped  boundary  layer  profiles  at  84.8%  chord  for  V-grooves  starting  at  60%. 
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The  2.77  mm  (0. 109  inches)  diameter  trip  wire  was  placed  at  the  35%  chord  line  upstream  of 
the  V-grooves  which  provided  aturbulent  boundaiy  layer.  While  the  trip  wire  did  provide  a  turbulent 
boundary  layer,  the  boundary  layer  did  not  separate  from  the  baseline  smooth  blade  as  the  turbulent 
boundary  layer  did  in  the  experiments  performed  by  Lin  et  al .  [  1 8] .  The  lack  of  separation  over  the 
blade  with  aturbulent  boundary  layer  prevented  reductions  in  the  blade  wake  when  V-grooves  were 
used  compared  to  a  smooth  blade  with  a  boundary  layer  trip  (see  Figures  86  and  87). 

9.7  Exit  Momentum  Deficit  Calculations  for  V-Grooves 

Given  the  exit  velocity  profiles  in  Figures  84  through  87,  the  wake  momentum  deficit  (©)  for 
the  V-grooves  was  calculated.  Figure  88  illustrates  wake  momentum  deficits  for  the  V-groove  blade 
with  a  baseline  blade  comparison  provided.  For  the  case  of  no  boundary  layer  trip,  both  V-groove 
configurations  provided  improvements  ranging  from  3%  to  25.5%  in  the  wake  momentum  deficit 
at  Reynolds  numbers  of  100k  and  200k  with  a  freestream  turbulence  level  of  1%.  For  a  Reynolds 
number  of  50k,  only  the  V-grooves  starting  at  60%  chord  provided  a  reduction  in  0  by  30.2%  and 
10.7%,  for  freestream  turbulence  levels  of  1%  and  4%,  respectively. 

Figure  88  also  contains  ©  values  for  the  V-grooves  tested  with  a  boundary  layer  trip  at  35% 
axial  chord.  The  wake  momentum  deficit  values  for  the  V-grooves  were  higher  than  the  baseline 
tripped  blade  for  Reynolds  numbers  of  100k  and  200k.  At  a  Reynolds  number  of  50k,  both  the  V- 
grooves  starting  at  55%  and  60%  axial  chord  provided  a  small  reduction  in  ©  by  11.1%  and  5.2%, 
respectively,  for  a  freestream  turbulence  level  of  1%.  At  the  higher  freestream  turbulence  level  of 
4%,  no  reduction  in  wake  momentum  deficit  was  measured  for  the  V-grooves. 

9.8  Loss  Coefficient  Measurement  for  V-Grooves 

The  loss  coefficient  was  calculated  using  the  previously  described  exit  velocity  profiles  along 
with  the  measured  total  pressure  in  the  wake  of  the  V-groove  blade.  The  total  pressure  probe  was 
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Velocity  (m/s) 


pitchwlse  traverse  (cm) 

a)  Exit  velocities  for  Tu  =  1% 


pitchwise  traverse  (cm) 

b)  Exit  velocities  for  Tu  =  4% 


Figure  84.  Exit  velocity  profiles  for  V-grooves  starting  at  55%  chord. 
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Velocity  (m/s)  Velocity  (m/s) 


30  i 


pitchwise  traverse  (cm) 


a)  Exit  velocities  for  Tu  =  1% 


pitchwise  traverse  (cm) 


b)  Exit  velocities  for  Tu  =  4% 

Figure  85.  Exit  velocity  profiles  for  V-grooves  starting  at  60%  chord. 
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Velocity  (m/s)  Velocity  (m/s) 


pitchwise  traverse  (cm) 

a)  Exit  velocities  for  Tu  =  1% 


pitchwise  traverse  (cm) 


b)  Exit  velocities  for  Tu  =  4% 

Figure  86.  Exit  velocity  profiles  for  V-grooves  starting  at  55%  chord  with  boundary  layer  trip. 
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Velocity  (m/s)  Velocity  (m/s) 


b)  Exit  velocities  for  Tu  =  4% 


Figure  87.  Exit  velocity  profiles  for  V-grooves  starting  at  60%  chord  with  boundary  layer  trip. 
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Figure  88.  Wake  momentum  values  for  baseline  blade,  and  V-grooves  starting  at  55%  and  60% 
chord. 
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moved  across  the  V-groove  blade  wake  approximately  ^  axial  chord  lengths  downstream.  Figures 
89  and  90  show  the  exit  total  pressure  measured  for  the  V-grooves  with  an  initially  laminar  boundary 
layer.  The  size  of  the  wakes  were  reduced  by  a  small  amount  compared  to  the  baseline  smooth  blade. 
For  the  case  of  a  tripped  boundary  layer  the  total  pressure  wake  of  the  V-grooves  was  greater  than 
the  smooth  blade  with  atrip  (see  Figures  91  and  92). 

Using  the  exit  velocity  and  total  pressure  measurements  from  Figures  84  through  92,  the  loss 
coefficient  was  calculated  for  the  V-grooves.  Figures  93  and  94  illustrate  the  two  very  different 
results  V-grooves  provide  depending  upon  the  initial  state  of  the  boundary  layer.  As  seen  in  Figure 
93,  the  V-grooves  starting  at  60%  axial  chord  provided  a  beneficial  reduction  in  the  2-D  profile 
loss  coefficient  for  Reynolds  numbers  of  50k  and  100k.  The  maximum  reduction  in  loss  coefficient 
for  the  V-grooves  starting  at  60%  occurred  at  a  Reynolds  number  of  50k  and  freestream  turbulence 
level  where  the  loss  was  reduced  by  36.8%.  The  V-grooves  starting  at  55%,  however,  only  provided 
a  reduction  in  the  loss  coefficient  by  5 . 1  %  to  1 3 .2%  for  a  Reynolds  number  of  1 00k  (for  freestream 
turbulence  levels  of  1%  and  4%),  and  a  reduction  in  loss  coefficient  of  22.4%  at  a  Reynolds  number 
of  50k  when  the  inlet  freestream  turbulence  level  was  4%.  Both  V-groove  modifications  provided 
essentially  the  same  loss  coefficient  as  the  baseline  blade  at  a  Reynolds  number  of  200k  for  the 
initially  laminar  boundary  layer. 

For  the  tripped  boundary  layer,  the  V-grooves  increased  the  2-D  loss  coefficient  slightly  com¬ 
pared  to  the  baseline  smooth  blade  with  a  boundary  layer  trip.  As  seen  in  Figure  94,  the  loss  coeffi¬ 
cient  increases  with  Reynolds  number  for  the  baseline  blade  and  the  V-grooves  when  the  boundary 
is  tripped  with  a  wire. 
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pitchwise  traverse  (cm) 

a)  Exit  pressure  profiles  for  Tu  =  1  % 


pitchwise  traverse  (cm) 


b)  Exit  pressure  profiles  for  Tu  =  4% 


Figure  89.  Exit  stagnation  pressure  loss  for  V-grooves  starting  at  55%  chord. 
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Figure  91.  Exit  stagnation  pressure  loss  for  V-grooves  starting  at  55%  chord  with  boundary  layer 
trip. 


131 


□  Re  *  50k 
o  Re = 100k 
a  Re  =  200k 


X 

l 

I 

® 

c 

l 


120  1 - - - - - - - - - - 

100 - jpt - 

▲ 

80 - — - 

A 

A 

60 - 

a  A 

40 - 

A  1 

20 - — - 

A  A 

0  . . I 

60  62  64  66  68  70  72  74  76  78  80 


A 

A 

A 

A 

A 

i 

A 

A 

k 

A 

Mil 

♦ 

[■■■■ 

*♦ 

♦ 

■  ■id 

A 

.  r . 

■  Re  =  50k 
♦  Re  =  100k 
a  Re  =  200k 


pitchwise  traverse  (cm) 


b)  Exit  pressure  profiles  for  Tu  =  4% 


Figure  92.  Exit  stagnation  pressure  loss  for  V-grooves  starting  at  60%  chord  with  boundary  layer 
trip. 
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Figure  93.  Loss  coefficient  for  baseline  blade  and  V-grooves  starting  at  55%  and  60%  chord. 
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Figure  94.  Loss  coefficient  for  baseline  blade  and  V-grooves  starting  at  55%  and  60%  chord  with 
boundary  layer  trip. 
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Chapter  10  -  Summary  and  Conclusions 

Unmanned  aerial  vehicles  operating  at  high  altitudes  and  low  cruise  velocities  experience  low 
Reynolds  number  flow  over  the  aircraft  and  within  the  engine.  The  low  pressure  turbine  (LPT) 
blades  within  the  engine  may  suffer  lower  efficiencies  than  predicted  while  the  aircraft  is  performing 
its  mission  phase  at  low  Reynolds  numbers.  The  lower  engine  efficiencies  are  due  to  flow  separating 
from  the  surface  of  the  LPT  blades.  The  laminar  separation  problem  was  verified  to  exist  with  the 
Pratt  and  Whitney  Pak-B  blade. 

A  new  linear  cascade  tunnel  was  used  to  simulate  the  low  pressure  turbine  environment.  The 
test  section  consisted  of  eight  Pak-B  blades  and  nine  passages.  The  blades  were  tested  for  inlet 
Reynolds  numbers  of  50k,  1 00k  and  200k  and  for  freestream  turbulence  levels  of  1  %  and  4%.  These 
conditions  cover  the  range  under  which  the  laminar  separation  problem  is  most  prevalent  within 
the  jet  engine.  The  efficiency  of  the  Pak-B  blades  was  measured  in  terms  of  the  2-D  profile  loss 
coefficient.  The  higher  the  profile  loss  coefficient,  the  lower  the  engine  operating  efficiency.  For 
low  Reynolds  numbers  and  turbulence  level  (Re  =  50k  and  Tu  =  1%)  the  loss  coefficient  is  largest 
due  to  the  flow  separating  without  reattachment  to  the  blade. 

Previous  attempts  to  correct  for  the  decrease  in  efficiency  have  been  extensions  of  known  active 
control  methods,  such  as  vortex  generator  jets.  Passive  techniques,  however,  have  the  advantage 
of  not  requiring  additional  power  or  controls.  Installation  into  a  pre-existing  engine  would  mean 
replacement  of  the  LPT  blade  spool  with  a  set  of  modified  blades.  No  additional  controls  or  power 
would  be  required  for  installation.  No  information  has  been  previously  published  on  the  successful 
implementation  of  a  passive  separation  control  technique  for  a  jet  engine  blade  until  this  document. 

Of  the  twenty-one  different  passive  separation  control  techniques  initially  considered,  three 
were  chosen  for  further  examination.  The  first  passive  technique  was  the  attachment  of  a  wire  to 
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the  suction  surface  of  the  blade  at  the  35%  axial  chord  line.  The  initially  laminar  flow  was  tripped 
to  a  turbulent  state.  The  turbulent  flow  remained  attached  over  the  surface  of  the  blade  downstream 
of  the  wire  for  all  conditions  tested.  While  the  trip  wire  reduced  the  loss  at  a  Reynolds  number  of 
50k  and  freestream  turbulence  level  of  1%,  at  all  other  conditions  the  trip  wire  increased  the  losses 
for  the  blade. 

The  second  passive  technique  consisted  of  dimples  machined  into  the  suction  surface  of  the 
airfoil.  The  dimples  were  successful  in  reducing  the  loss  coefficient  at  low  Reynolds  numbers,  while 
not  degrading  performance  at  the  highest  Reynolds  number  tested  (see  Figure  95).  Three  variations 
of  dimples  were  tested,  the  chordwise  location  of  the  dimples  in  each  test  being  changed  to  determine 
the  effect  of  dimple  location  on  loss  reduction.  The  baseline  separation  location  moves  downstream 
with  increasing  Reynolds  number  and  turbulence  (62%  to  78%  axial  chord).  The  most  effective 
dimples  are  those  upstream  of  separations  located  near  separation.  Dimples  provided  improved 
losses  at  all  locations,  including  the  location  slightly  downstream  of  separation,  the  low  Reynolds 
number,  low  turbulence  case  (50k,  1%).  The  improved  loss  for  the  dimples  downstream  of  the 
separation  point  was  due  to  the  thin  separation  (less  than  0.16)  at  that  location,  which  allowed  for 
partial  effectiveness  of  the  dimples. 

To  take  advantage  of  this  new  knowledge,  guidance  must  be  given  for  proper  placement  of  the 
dimples  on  the  suction  surface.  The  loss  coefficient  was  plotted  for  variation  in  axial  chord  distance 
between  baseline  separation  locations  and  dimples  in  Figure  96.  Figures  96a,  b  and  c  also  show 
the  baseline  blade  loss  coefficients.  As  seen  in  Figure  96,  a  benefit  was  attained  with  the  dimples 
at  all  three  chord  locations.  Furthermore,  for  higher  Reynolds  numbers  (comparing  Figures  96a,  b 
and  c)  the  dimples  reduce  the  loss  coefficient  to  approximately  the  same  value  regardless  of  their 
location.  This  suggests  that  with  a  reduction  in  separation,  losses  are  similar  and  due  mainly  to 

viscous  effects  which  decrease  slightly  with  Reynolds  number.  For  turbulent  reattachment,  the  skin 
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friction  coefficient  is  proportional  to  1/Re1/7,  for  which  the  values  of  loss  coefficient  scale  as  seen  in 
Figure  96  (except  for  the  low  Reynolds  number  case  for  which  the  flow  a  separation  zone  remains). 

The  improvements  in  loss  coefficient  can  be  related  to  improvements  in  engine  stage  efficiency. 
For  the  reductions  in  loss  coefficient  (18.8%  to  51.7%)  made  possible  through  the  use  of  dimples, 
equation  10  of  Chapter  1  indicates  a  typical  increase  of  turbine  stage  efficiency  of  1.1%  to  2.9%. 
This  improved  efficiency  of  2.9%  can  be  expected  to  improve  the  specific  thrust  and  specific  fuel 
consumption  (SFC)  on  a  turbojet  engine  at  Mach  2, 36,000  ft.  (Mattingly  [20])  by  0.54%  and  0.54%, 
respectively.  A  decrease  of  0.54%  in  SFC  would  increase  range  by  about  1.41%  and  a  reduction  in 
fuel  consumption  by  0.59%. 
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Figure  96.  Effect  of  dimple  distance  from  baseline  separation  on  loss  coefficient. 
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Streamwise  V-grooves  were  tested  in  two  different  configurations.  The  axial  chord  starting 
location  of  the  V-grooves  was  varied  to  determine  the  effect  on  loss  reduction.  For  an  initially  lam¬ 
inar  boundary  layer,  the  V-grooves  starting  at  55%  chord  reduced  the  loss  coefficient  for  a  moder¬ 
ate  Reynolds  number  and  low  turbulence  intensity  and  low  and  moderate  Reynolds  numbers  (50k 
and  100k)  with  high  turbulence  intensity.  The  V-grooves  starting  at  60%  chord  reduced  the  loss 
coefficient  for  low  and  moderate  Reynolds  numbers  (see  Figure  97)  at  low  turbulence,  and  for  all 
Reynolds  numbers  at  high  turbulence.  Therefore  the  V-grooves  are  more  beneficial  when  placed 
closer  to  the  baseline  separation  point. 

The  two  sets  of  V-grooves  were  also  tested  with  a  boundary  layer  trip  located  at  the  35%  axial 
chord  line.  A  comparison  of  the  results  with  those  of  the  baseline  blade  with  a  boundary  layer 
tnP  showed  no  net  reduction  in  losses  with  a  freestream  turbulence  level  of  1%.  At  the  higher 
turbulence  level  of  4%,  the  V-grooves  provided  a  0.8%  reduction  in  loss  coefficient  at  the  highest 
Reynolds  number  (200k).  In  summary,  the  trip  wire  and  V-groove  experiments  showed  a  decrease 
in  loss  coefficient  for  certain  limited  cases  of  Reynolds  number  and  freestream  turbulence  level. 
Any  passive  technique  for  separation  control  must  provide  a  reduction  in  the  loss  coefficient  at  low 
Reynolds  numbers  while  not  significantly  increasing  the  losses  at  higher  Reynolds  numbers  typical 
of  take-off  operations.  Considering  the  performance  of  the  V-grooves,  only  the  V-grooves  starting 

at  60%  meet  this  criteria  while  the  trip  wire  with  V-groove  combination  should  not  be  used  in  a  jet 
engine. 

In  conclusion,  the  dimples  consistently  showed  improved  performance  over  the  baseline  blade 
and  were  the  best  passive  technique  examined  for  application  to  low  pressure  turbine  blades.  The 
dimples  unproved  the  low  Reynolds  number  performance  by  eneigizing  the  flow  sufficiently  to  pre¬ 
vent  separation.  At  higher  Reynolds  numbers  the  submeiged  dimples  did  not  increase  the  boundary 
layer  thickness  as  the  trip  wire  did.  As  with  the  baseline  blade,  the  loss  coefficient  for  the  dimples 
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also  decreased  with  increasing  Reynolds  number  and  turbulence  intensity.  The  combination  of  sep¬ 
aration  prevention  at  low  Reynolds  numbers  and  efficient  flow  at  high  Reynolds  numbers  through 
a  passive  technique  requiring  no  additional  power  or  controls  makes  the  dimple  a  significant  new 
tool  for  the  engine  designer. 
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Figure  97.  Change  in  loss  coefficient  using  V-grooves  with  and  without  atrip  wire. 
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APPENDIX  A  -  Boundary  Layer  Measurements  for  Blade  with 

Dimples 

A.1  Velocity  and  Turbulence  Measurements  at  Re  =  50k  and  Re  =  100k 


b)  boundary  layer  turbulence  measurements  at  Re  =  50k  and  Tu  =  1% 


Figure  98.  Measured  velocity  and  turbulence  profiles  between  dimples  at  50%. 
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Velocity  (U/Umax  *  100) 

a)  boundary  layer  velocity  measurements  at  Re  =  50k  and  Tu  =  4% 


b)  boundary  layer  turbulence  measurements  at  Re  =  50k  and  Tu  =  4% 


Figure  99.  Measured  velocity  and  turbulence  profiles  behind  a  dimple  at  50%  with  Tu  =  4%. 
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Velocity  (U/Umax  *  100) 

a)  boundary  layer  velocity  measurements  at  Re  =  50k  and  Tu  =  4% 


b)  boundary  layer  turbulence  measurements  at  Re  =  50k  and  Tu  =  4% 


Figure  100.  Measured  velocity  and  turbulence  profiles  between  dimples  at  50%  with  Tu  =  4%. 
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Height  (mm)  Height  (mm) 


Velocity  (U/Umax  *  100) 

a)  boundary  layer  measurements  at  Re  =  50k  and  Tu  =  1% 


Turbulence  Intensity 

b)  boundary  layer  measurements  at  Re  =  50k  and  Tu  =  1% 


Figure  101.  Measured  velocity  and  turbulence  profiles  between  dimples  at  55%. 
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Height  (mm)  Height  (mm) 


Velocity  (U/Umax  *  100) 

a)  boundary  layer  measurements  at  Re  =  50k  and  Tu  =  4% 


Turbulence  Intensity 

b)  boundary  layer  measurements  at  Re  =  50k  and  Tu  =  4% 


Figure  102.  Measured  velocity  and  turbulence  profiles  behind  a  dimple  at  55%  with  Tu  =  4%. 
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Height  (mm)  Height  (mm) 


Velocity  (U/Umax  *  100) 

a)  boundary  layer  measurements  at  Re  =  50k  and  Tu  =  4% 


Turbulence  Intensity 

b)  boundaiy  layer  measurements  at  Re  -  50k  and  Tu  =  4% 


Figure  103.  Measured  velocity  and  turbulence  profiles  between  dimples  at  55%  with  Tu  =  4%. 
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Height  (mm)  Height  (mm) 


Velocity  (U/UmaxMOO) 

a)  boundary  layer  velocity  measurements  at  Re  =  50k  and  Tu  =  1% 


Turbulence  Intensity 


b)  boundary  layer  turbulence  measurements  at  Re  =  50k  and  Tu  =  1% 


Figure  104.  Measured  velocity  and  turbulence  profiles  between  dimples  at  65%. 
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Height  (mm)  Height  (mm) 


Velocity  (U/Umax  *  100) 

a)  boundary  layer  velocity  measurements  at  Re  =  50k  and  Tu  =  4% 


b)  boundary  layer  turbulence  measurements  at  Re  =  50k  and  Tu  =  4% 


Figure  105.  Measured  velocity  and  turbulence  profiles  behind  a  dimple  at  65%  with  Tu  =  4%. 
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Height  (mm)  Height  (mm) 


40 


Velocity  (U/Umax  *  100) 

a)  boundary  layer  velocity  measurements  at  Re  =  50k  and  Tu  =  4% 


Turbulence  Intensity 

b)  boundaiy  layer  turbulence  measurements  at  Re  =  50k  and  Tu  =  4% 


Figure  106.  Measured  velocity  and  turbulence  profiles  between  dimples  at  65%  with  Tu  =  4%. 
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Height  (mm)  Height  (mm) 


Velocity  (U/UmaxMOO) 

a)  boundary  layer  velocity  measurements  at  Re  =  100k  and  Tu  =  1% 


Turbulence  Intensity 

b)  boundary  layer  turbulence  measurements  at  Re  =  100k  and  Tu  =  1% 

Figure  107.  Measured  velocity  and  turbulence  profiles  between  dimples  at  50%. 


157 


Velocity  (U/Umax  *  100) 

a)  boundary  layer  velocity  measurements  at  Re  =  100k  and  Tu  =  4% 


b)  boundary  layer  turbulence  measurements  at  Re  =  100k  and  Tu  =  4% 


Figure  108.  Measured  velocity  and  turbulence  profiles  behind  a  dimple  at  50%  with  Tu  =  4%. 
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Velocity  (U/Umax  *  100) 

a)  boundary  layer  velocity  measurements  at  Re  =  100k  and  Tu  =  4% 


b)  boundary  layer  turbulence  measurements  at  Re  =  100k  and  Tu  =  4% 


Figure  109.  Measured  velocity  and  turbulence  profiles  between  dimples  at  50%  with  Tu  =  4%. 
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Height  (mm) 


Velocity  (U/Umax  *  100) 

a)  boundary  layer  at  Re  =  100k  and  Tu  =  1% 


Turbulence  Intensity 

b)  boundary  layer  at  Re  =  100k  and  Tu  =  1% 


Figure  110.  Measured  velocity  and  turbulence  profiles  between  dimples  at  55%. 
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Height  (mm)  Height  (mm) 


Velocity  (U/Umax  *  100) 

a)  boundary  layer  at  Re  =  100k  and  Tu  =  4% 


Turbulence  Intensity 


b)  boundary  layer  at  Re  =  100k  and  Tu  =  4% 


Figure  111.  Measured  velocity  and  turbulence  profiles  behind  a  dimple  at  55%  with  Tu  =  4%. 
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Height  (mm)  Height  (mm) 


Velocity  (U/Umax  *  100) 

a)  boundary  layer  at  Re  =  100k  and  Tu  =  4% 


Turbulence  Intensity 


b)  boundary  layer  at  Re  =  100k  and  Tu  =  4% 


Figure  112.  Measured  velocity  and  turbulence  profiles  between  dimples  at  55%  with  Tu  =  4%. 


162 


Height  (mm)  Height  (mm) 


Velocity  (U/Umax  *  100) 


a)  boundary  layer  velocity  measurements  at  Re  =  100k  and  Tu  =  1% 


Turbulence  Intensity 

b)  boundary  layer  velocity  measurements  at  Re  =  100k  and  Tu  =  1% 


Figure  113.  Measured  velocity  and  turbulence  profiles  between  dimples  at  65%. 
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Height  (mm)  Height  (mm) 


Velocity  (U/Umax  *  100) 


a)  boundaiy  layer  velocity  measurements  at  Re  =  100k  and  Tu  =  1% 


Turbulence  intensity 


b)  boundary  layer  velocity  measurements  at  Re  =  100k  and  Tu  =  4% 

Figure  114.  Measured  velocity  and  turbulence  profiles  behind  a  dimple  at  65%  with  Tu  =  4%. 
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Height  (mm)  Height  (mm) 


Velocity  (U/Umax  *  100) 

a)  boundary  layer  velocity  measurements  at  Re  =  100k  and  Tu  =  4% 


Turbulence  Intensity 

b)  boundary  layer  turbulence  measurements  at  Re  =  100k  and  Tu  =  4% 


Figure  115.  Measured  velocity  and  turbulence  profiles  between  dimples  at  65%  with  Tu  =  4%. 
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APPENDIX  B  -  Boundary  Layer  Measurements  for  Blade  with 


V-Grooves 


B.l  \felocity  and  Turbulence  Measurements  at  Re  =  50k  and  Re  =  100k 


Velocity  (U/Umax  *  100) 

a)  boundary  layer  velocity  measurements  at  Re  =  50k  and  Tu  =  1% 


c  . . .  )  i  i  i  i  i  »  i  i  i  i  i  i  i  i  t  i  i  i  i  i  i  i  i  . . .  i  i  i  i  I 

0%  100%  100%  100%  100%  100% 
Turbulence  Intensity 

b)  boundary  layer  turbulence  measurements  at  Re  =  50k  and  Tu  =  1% 


Figure  116.  Measured  velocity  and  turbulence  profiles  inside  a  V-groove  starting  at  55%. 
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40 


Velocity  (U/Umax  *  100) 

a)  boundary  layer  velocity  measurements  at  Re  =  50k  and  Tu  =  4% 


b)  boundary  layer  turbulence  measurements  at  Re  =  50k  and  Tu  =  4% 


Figure  117.  Measured  velocity  and  turbulence  profiles  between  V-grooves  starting  at  55%  with  Tu 
=  4%. 
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Velocity  (U/Umax  *  100) 


a)  boundary  layer  velocity  measurements  at  Re  =  50k  and  Tu  =  4% 
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b)  boundary  layer  turbulence  measurements  at  Re  =  50k  and  Tu  =  4% 


Figure  118.  Measured  velocity  and  turbulence  profiles  inside  a  V-groove  starting  at  55%  with  Tu 
4%. 
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Height  (mm)  Height  (mm) 


Velocity  (U/Umax  *  100) 


a)  boundary  layer  velocity  measurements  at  Re  =  50k  and  Tu  =  1% 


Turbulence  Intensity 

a)  boundary  layer  velocity  measurements  at  Re  =  50k  and  Tu  =  1% 


Figure  119.  Measured  velocity  and  turbulence  profiles  inside  a  V-groove  starting  at  60%. 
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Velocity  (U/Umax  *  100) 

a)  boundary  layer  velocity  measurements  at  Re  =  50k  and  Tu  =  4% 
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b)  boundary  layer  turbulence  measurements  at  Re  =  50k  and  Tu  =  4% 


Figure  120.  Measured  velocity  and  turbulence  profiles  between  V-grooves  starting  at  60%  with  Tu 
=  4%. 
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Turbulence  Intensity 

b)  boundary  layer  turbulence  measurements  at  Re  =  50k  and  Tu  =  4% 


Figure  121.  Measured  velocity  and  turbulence  profiles  inside  a  V-groove  starting  at  60%  with  Tu  = 
4%. 
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Height  (mm)  Height  (mm) 


0%  100  100  100  100  100 

Velocity  (U/Umax*100) 

a)  boundary  layer  velocity  measurements  at  Re  =  100k  and  Tu  =  1% 


T urbulence  Intensity 

b)  boundary  layer  turbulence  measurements  at  Re  =  100k  and  Tu  =  1% 


Figure  122.  Measured  velocity  and  turbulence  profiles  inside  a  V-groove  starting  at  55%. 
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Velocity  (U/Umax  *  100) 

a)  boundary  layer  velocity  measurements  at  Re  =  100k  and  Tu  =  4% 


Turbulence  Intensity 

b)  boundary  layer  turbulence  measurements  at  Re  =  100k  and  Tu  =  4% 


Figure  123.  Measured  velocity  and  turbulence  profiles  between  V-grooves  starting  at  55%  with  Tu 
=  4%. 
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a)  boundary  layer  velocity  measurements  at  Re  =  100k  and  Tu  =  4% 
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Figure  124.  Measured  velocity  and  turbulence  profiles  inside  a  V-groove  starting  at  55%  with  Tu 

4%. 
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Height  (mm)  Height  (mm) 


0%  100  100  100  100  100 
Velocity  (U/Umax  *  100) 


a)  boundary  layer  velocity  measurements  at  Re  =  100k  and  Tu  =  1% 
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b)  boundary  layer  turbulence  measurements  at  Re  =  100k  and  Tu  =  1% 


Figure  125.  Measured  velocity  and  turbulence  profiles  inside  a  V-groove  starting  at  60%. 
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Velocity  (U/Umax  *  100) 

a)  boundary  layer  velocity  measurements  at  Re  =  100k  and  Tu  =  4% 


Turbulence  Intensity 

b)  boundary  layer  turbulence  measurements  at  Re  =  100k  and  Tu  =4% 


Figure  126.  Measured  velocity  and  turbulence  profiles  between  V-grooves  starting  at  60%  with  Tu 
=  4%. 
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Velocity  (U/Umax  *  100) 


a)  boundary  layer  velocity  measurements  at  Re  =  100k  and  Tu  =  4% 
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Figure  127.  Measured  velocity  and  turbulence  profiles  inside  a  V-groove  starting  at  60%  with  Tu 

4%. 
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APPENDIX  C  -  Adjustable  Linear  Cascade  Experimental  Facility 

C.l  Baseline  Wind  Tunnel  Design 

The  wind  tunnel  used  in  the  experiment  was  a  modified  open  loop  Aerolab  corporation  draw¬ 
down  tunnel  using  a  test  section  with  an  interior  area  121.9  cm  (48  inches)  wide  by  85.1  cm  (33.5 
inches)  high  (see  Figure  128).  Air  was  pulled  through  the  test  section  by  means  of  a  93.2  kilowatt 
(125  HP)  motor  operating  in  the  suction  mode  using  a  Joy  Technologies  Axivane  axial  flow  fan  with 
an  outer  radius  of  114.3  cm  (45  inches).  The  motor  casing  is  124.5  cm  (49  inches)  long,  containing 
the  electric  motor,  axial  fan,  exit  vanes,  and  a  wire  mesh  catcher  system  to  minimize  blade  damage 
in  the  event  of  a  test  section  structural  failure.  Air  flow  through  the  test  section  is  controlled  by  a 
Harmon/Commonwealth  Corporation  variable  frequency  motor  controller. 


The  wind  tunnel  inlet  bell-mouth  has  a  304.8  cm  (120  inch)  wide  x  266.7  cm  (105  inch)  high 
rectangular  inlet  with  a  honeycomb  flow  straightener  and  a  converging  nozzle.  The  convergent 
portion  of  the  inlet  nozzle  is  228.6  cm  (90  inches)  long  and  ends  in  a  rectangular  shape  with  the 
center  outlet  flow  121 .9  cm  (48  inches)  wide  by  85. 1  cm  (33.5  inches)  high  (see  Figure  129). 

Inserts  serve  to  minimize  comer  vortices  and  extend  into  the  test  section  10.2  cm  (4  inches) 
up  and  down  from  the  top  and  bottom  and  12.7  cm  (5  inches)  into  it  from  the  left  and  right.  The 
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standard  straight  test  section  tapers  these  inserts  down  to  6.4  cm  (2.5  inches)  high  in  the  vertical 
axis  and  8.9  cm  (3.5  inches)  wide  in  the  horizontal  axis.  With  these  inserts  in  place,  the  freestream 
turbulence  level  is  less  than  1 .0%. 

Modifications  built  by  the  AFIT  shop  for  the  wind  tunnel  did  not  include  the  comer  inserts  to 
minimize  shop  materials  and  construction  time.  To  replicate  their  effectiveness,  however,  styrofoam 
inserts  were  cut  and  placed  in  the  new  tunnel  sections  to  maintain  the  same  comer  vortex  minimizing 
effect.  With  these  inserts  in  place,  the  new  test  section  has  a  freestream  turbulence  of  approximately 
1%. 

C.2  Active/Passive  Turbulence  Generator 


Figure  130.  Turbulence  grid  capable  of  producing  high  turbulence  levels  with  and  without  blowing 
air. 


The  turbulence  generator  housing  consists  of  a  straight  channel  76.2  cm  (30  inches)  long  in 
the  flow  direction  and  has  the  same  cross-sectional  area  of  the  straight  test  section.  The  turbulator 
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grid  itself  (see  Figure  130)  is  a  removable  square-mesh  array  of  2.54  cm  (1  inch)  diameter  round 
stainless  steel  tubes  that  can  provide  both  passive  and  active  means  of  turbulence  generation.  The 
tubes  are  spaced  7.6  cm  (3  inches)  apart,  center-to-center  and  were  designed  in  accordance  with 
Roach’s  [48]  isotropic  turbulence  grid  guidelines. 

Operated  passively,  the  turbulator  disturbed  the  air  to  a  freestream  turbulence  level  of  3.5% 
to  4.5%  (see  Figure  131)  at  a  distance  of  127  to  254  cm  (50  -  100  inches)  downstream  of  the  grid. 
The  turbine  blade  cassette  is  located  within  this  downstream  distance  of  the  new  variable  angle  test 
section.  Also  included  in  Figure  131  are  the  decay  correlations  determined  from  experiments  of 
Baughn  et  al.  [10]  and  Roach  [48]. 

40 
35 
30 
25 

e£ 

3  20 
t- 

15 
10 
5 
0 

0  10  20  30  40  50  60  70  80  90  100 

Distance  (inches) 

Figure  131.  Decay  of  passive  grid  turbulence  with  downstream  position. 

Measurement  of  the  generated  turbulence  decay  and  length  scales  were  completed  inside  the 
new  adjustable  test  section.  The  recorded  turbulence  data  was  then  examined  to  ensure  proper 
isotropy  and  length  scales  in  accordance  with  Mayle  et  al.  [49],  Length  scales  were  determined  by 
taking  several  repeated  measurements  and  using  Wright  Laboratory  provided  software  to  integrate 
and  average  the  data.  The  turbulence  measurements  were  made  using  a  single  hot-wire  probe  im- 
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mersed  in  the  air  flow.  The  baseline  straight  test  section  with  the  turbulence  generator  in  place  is 
depicted  in  Figure  132. 

C.2.1  Active  Turbulence  Generation 

To  create  higher  turbulence  levels  across  the  blade  cassette,  air  can  be  blown  through  the  0.8 
mm  (1/32  inch)  diameter  holes  spaced  evenly  through  the  vertical  bars.  The  fireestream  turbulence 
level  can  be  increased  to  as  high  as  20%  at  a  distance  of  127  to  254  cm  (50  -  100  inches)  down¬ 
stream  of  the  grid.  Control  of  the  air  through  the  vertical  bars  would  be  maintained  through  the 
stainless  steel  Octopus  Air  Control  System  seen  in  Figure  133.  This  system  has  four  independent 
pressure  chambers,  each  fitted  with  a  pressure  gauge  and  four  air  exhaust  ports.  Accurate  control 
of  the  air  through  the  exhaust  ports  is  maintained  through  a  gross  adjustment  2 . 54  cm  ( 1  inch)  inner 
diameter  valve  and  a  1.27  cm  (0.5  inch)  inner  diameter  fine  control  needle  valve.  Air  is  passed  to 
the  turbulence  generating  grid  through  2.54  cm  (1  inch)  inside  diameter  high  pressure  hose  fitted 
with  quick  disconnects  at  both  ends.  Using  this  system  air  can  be  delivered  to  any  4,  8,  12,  or  all 
16  of  the  vertical  tubes  in  any  combination.  Active  turbulence  generation  using  the  blown  grid  re¬ 
quires  a  large  compressor  capable  of  producing  150+  psi  continuously.  A  large  compressor  was  not 
available  for  this  research  and  the  grid  was  only  used  passively  for  turbulence  levels  of  4%  inside 
the  new  adjustable  test  section. 

C.3  Adjustable  Inlet  and  Exit  Angle  Cascade  Test  Section 

The  test  section  contains  a  linear,  2-D  airfoil  cascade  consisting  of  8  low  pressure  turbine 
blades  with  9  airflow  passages  in  between  the  blades  and  walls.  Blade  #1  is  located  near  the  inside 
bend,  with  subsequent  blades  #2  through  #8  following  outward  towards  the  outer  comer.  The  blade 
row  is  198  cm  (78  inches)  downstream  of  the  turbulence  generator  at  its  nearest  point  and  300  cm 
(118  inches)  downstream  at  its  farthest.  The  Pak-B  blades  were  designed  by  Pratt  and  Whitney  and 
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have  axial  chord  lengths  of  17.8  cm  (7  inches).  The  axial  chord  to  spacing  ratio  (solidity,  |)  of  the 
blade  passage  is  1.129.  The  blade  spacing  is  15.75  cm  (6.2  inches)  blade  to  blade.  Each  blade  is 
molded  from  Ultralloy  108  white  resin  and  is  seven  times  engine  scale  in  the  axial  direction.  To  fur¬ 
ther  enhance  the  approximately  2-D  flowfield  of  the  midspan  of  the  actual  blades,  each  test  blade 
is  stretched  in  the  z-axis  to  78.63  cm  (34.5  inches)  high  for  an  aspect  ratio  of  4.92  to  1.  The  inner 
and  outer  walls  of  the  test  section  where  the  test  blades  are  located  were  shaped  to  approximate  the 
suction  and  pressure  sides  of  a  blade,  respectively.  The  tunnel  was  set  for  the  on-design  specifica¬ 
tions  of  the  Pak-B  blade  with  an  air  inlet  angle  of  35°  relative  to  the  engine  x-axis  (normal  to  the 
blade  inlet  plane)  and  an  air  exit  angle  of  60°. 

The  test  section  contains  a  central  pivot  joint  which  allows  adjustment  of  the  inlet  and  outlet 
channels  separately.  The  inlet  and  outlet  channels  have  the  same  cross-sectional  area  as  the  turbu¬ 
lence  generator,  and  are  connected  at  the  pivot  joint.  Both  the  inlet  and  exit  channel  may  be  turned 
through  a  range  of  0°  to  65°  independently  of  one  another.  The  maximum  total  possible  flow  turn¬ 
ing  angle  is  130°  if  both  channels  were  fully  rotated.  This  new  test  section  will  be  used  with  other 
blade  types  in  future  experiments  requiring  a  change  in  inlet  and  exit  angle.  To  facilitate  test  section 
adjustments,  the  top,  bottom,  and  outer  walls  of  the  test  section  were  built  using  clear  plexiglass. 

The  top  of  the  new  test  section  was  constmcted  from  two  sheets  of  2.54  cm  (1  inch)  thick  clear 
plexiglass.  These  two  laige  plexiglass  pieces  are  attached  to  the  top  with  a  groove  6.4  mm  (1/4 
inch)  deep  cut  around  the  edge.  This  groove  allows  the  plexiglass  to  hang  into  the  test  section  and 
create  a  flush  inner  surface.  The  good  visibility  within  the  test  section  allows  for  easy  adjustment  of 
probes  and  tailboards  when  necessary.  The  plexiglass  bottom  of  the  new  test  section  was  similarly 
constructed. 

Traverse  slots  were  cut  into  the  plexiglass  sheets,  arranged  parallel  to  the  leading  and  trailing 
edge  of  the  blade  row.  The  center  of  the  1.27  cm  (0.5  inches)  wide  slots  are  spaced  \  and  1  axial 
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Figure  134.  New  linear  cascade  test  section  with  independently  adjustable  inlet  and  exit  angle  ca¬ 
pability. 
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chord  length  away  from  the  blades.  The  slots  allowed  for  easy  measurement  of  the  inlet  and  exit 
profiles  of  the  blade  cassette  to  verify  periodicity,  a  requirement  suggested  by  Gostelow  [50],  Once 
periodicity  had  been  acquired  within  the  test  section  through  adjustment  of  the  tailboards  testing 
was  started  on  the  Pak-B  airfoils.  The  same  procedure  will  be  conducted  to  verify  periodicity  for 
other  airfoil  shapes  when  the  Pak-B  airfoils  are  replaced. 

C.4  Removable  Blade  Cassette 

The  blade  cassette  which  held  the  Pak-B  airfoil  cascade  during  testing.  Seen  in  Figure  135, 
the  blade  cassette  can  be  removed  from  the  test  section.  The  blade  cassette  was  also  constructed 
from  clear  plexiglass  which  allowed  for  improved  visibility  when  inserting  boundary  layer  probes 
into  the  test  section.  The  ease  with  which  the  cassette  can  be  installed  and  removed  facilitates  the 
testing  of  a  variety  of  blade  types  with  different  inlet  and  exit  flow  angles. 


X 

upper  clear  plate 


lower  clear  plate 


Figure  135.  Side  view  of  blade  cassette 
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The  blade  cassette  consists  of  the  eight  Pak-B  blades  and  clear  plexiglass  top  and  bottom  plates. 
To  support  the  weight  of  the  blade  cassette  steel  channels  fixed  to  the  adjustable  test  section  support 
the  leading  and  trailing  edge  of  the  cassette  plexiglass  top  and  bottom  plates.  The  blade  cassette 
design  is  compatible  with  any  other  low  pressure  turbine  blade  shape. 

C.4.1  Testing  Other  Blade  Shapes 

Once  a  new  cassette  is  created  with  the  new  blades  to  be  tested,  the  current  blade  cassette  is 
removed  by  first  removing  the  tailboards  from  the  test  section.  Then  the  current  cassette  is  lifted  out 
using  a  crane  and  the  new  cassette  lowered  into  place.  The  test  section  is  then  rotated  to  provide  the 
correct  exit  angle  for  the  blades  with  the  tail  boards  inserted  accordingly.  The  bell-mouth  inlet  is  then 
lifted  slightly  10  mm)  using  a  crane  and  rotated  to  provide  the  desired  inlet  angle.  Rotating  the 
bell-mouth  inlet  also  rotates  the  turbulence  generator  and  the  inlet  channel  of  the  new  test  section. 
Since  the  bell-mouth  is  not  bolted  to  the  floor,  it  can  remain  attached  to  the  test  section  during  this 
procedure. 
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APPENDIX  D  -  Experimental  Data  Collection  Software 


D.l  Overview  of  Software  and  Equipment 

Control  and  data  collection  for  the  wind  tunnel  was  managed  through  an  IBM-compatible  pc 
using  an  Intel  200  MHz  pentium  processor  with  a  2.0  Gb  hard  drive.  The  operating  software  was 
written  in  G,  a  visual  programming  language  used  by  the  Lab  View  5.0  software  package.  This 
software  package  included  pre-written  procedures  and  the  ability  to  use  them  with  operator  created 
software.  Each  operator  created  program  is  called  a  “VI”  (Virtual  Instrument)  and  has  a  display  on 
the  monitor  which  presents  both  the  entered  and  collected  data. 

D.2  Hot  Wire  Calibration 

The  HotWire  Calibration  VI  allowed  calibration  of  the  hot  wire  to  be  performed  within  the  new 
test  section  of  the  tunnel.  With  the  hot  wire  in  place,  measurement  of  the  probe  resistance  was  made 
using  the  TSI IFA  1 00  Intelligent  Flow  Analyzer.  This  resistance  was  entered  in  the  hot  wire  settings 
frame  on  the  left  hand  side  of  the  screen  (Figure  136)  along  with  the  factory  measured  internal 
resistance  of  the  hot  wire  probe.  The  current  atmospheric  pressure  was  read  from  a  Wallace  and 
Tieman  FA  129  pressure  gauge  and  entered  into  the  program,  while  the  temperatures  were  measured 
using  J-type  thermocouples.  To  calibrate  the  probe,  a  pitot-static  probe  was  inserted  into  the  tunnel 
1  axial  chord  upstream  of  Blade  #5.  This  pitot  probe  provided  the  Ptotai  -  Pstatic  measurement 
necessary  for  the  calibration.  With  the  pitot  and  hot  wire  probes  in  place,  the  tunnel  was  run  at  a 
high  velocity  (higher  than  that  used  during  data  collection)  and  the  resulting  pressure  from  the  pitot 
probe  and  voltage  from  the  hot  wire  entered  into  the  program. 

The  pressure  and  voltage  entry  procedure  was  repeated  for  several  different  decreasing  veloci¬ 
ties.  When  sufficient  data  was  collected,  the  “Done  Do  Fit”  button  was  pressed  on  the  screen.  This 
button  activated  a  sub-routine  to  perform  a  curve  of  the  hot  wire  data  with  the  pressure  data  using 
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a  method  described  by  Bruun  [31],  The  resulting  hot  wire  calibration  curve  was  a  nondimensional 
fit  between  the  Nusselt  number  and  Reynolds  number  of  the  form: 

Nu  =  CINT  +  Dslo  Re  0  45  (26) 

where  Cjnt  and  Dslo  were  the  intercept  and  slope,  respectively,  for  the  curve.  The  exponent  used 
for  the  hot  wire  equations  was  0.45  since  this  provided  better  agreement  with  the  low  velocities  used 
in  this  research  (Giese  [51]).  The  resulting  intercept  and  slope  were  then  entered  into  the  Current 
Tunnel  Conditions  and  other  Vi’s  where  hot  wire  velocity  information  was  collected. 

D.3  Current  Tunnel  Conditions 

The  initial  environmental  conditions  were  determined  before  every  data  collection.  This  was 
done  by  reading  the  atmospheric  pressure  from  a  Mhllace  and  Tieman  FA- 1 29  pressure  gauge  and  ob¬ 
taining  the  current  humidity  using  an  Omega  Engineering,  Inc.  RH71  Digital  Hygro-Thermometer. 
These  values  were  then  entered  into  the  Current  Tunnel  Conditions  VI  which  calculated  and  dis¬ 
played  the  current  tunnel  velocity,  air  density,  and  temperatures  throughout  the  tunnel  using  the  hot 
wire,  pressure  gauge,  and  thermocouples,  respectively.  The  desired  Reynolds  number  for  the  test 
was  then  dialed  in  by  increasing  the  tunnel  motor  speed  control  according  to  Figure  137.  \felocity 
data  was  collected  using  a  TSI  model  1210-T1 .5  hot  wire  described  previously  in  Section  4. 1 .  Each 
data  point  for  the  hot  wire  was  graphed  on  the  computer  screen  to  allow  easy  determination  of  AC 
noise  interference.  For  velocity  data  collected  behind  the  blades,  a  pitot  probe  was  used  to  deter¬ 
mine  the  inlet  velocity  and  Reynolds  number.  Figure  138  illustrates  the  Current  Tunnel  Conditions 
program. 

The  lower  left  part  of  Figure  138  shows  the  inlet  Reynolds  number  based  upon  the  pitot  probe 
measurement.  The  tunnel  pitot  probe  was  located  in  a  fixed  position  1  axial  chord  length  upstream 
of  the  blade  cassette  and  outside  of  the  boundary  layer.  The  hot  wire  calibration  information  is  just 
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.E?j  hot  wire  probe  calibration  -  a.VI ' 


Figure  136.  Software  used  for  calibra! 
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Figure  137.  Relationship  between  tunnel  inlet  velocity,  Reynolds  number,  and  drive  motor  setting 
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above  the  pitot  probe  data.  The  data  seen  in  Figure  138  was  taken  with  the  hot  wire  downstream  of 


the  blade  cassette  and  an  appropriate  velocity  increase  is  shown  between  the  upstream  pitot  velocity 
and  the  downstream  hot  wire  velocity.  All  Reynolds  number  information  presented  in  this  document 
was  based  on  inlet  velocity  and  an  axial  chord  length  of  17.78  cm  (7.0  inches).  The  calculated 


mean  velocity  and  turbulence  intensity  are  also  displayed  on  the  lower  right  portion  of  the  screen. 
Turbulence  intensity  was  based  upon  fluctuations  in  the  hot  wire  velocity  readings.  Temperature 
measurements  were  recorded  at  four  locations  within  the  tunnel  test  section. 
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Figure  138.  Current  tunnel  conditions  determination  program. 


Thermocouples  were  used  to  measure  the  temperatures  for  the  tunnel  inlet,  1 .45  m  (57  inches) 
upstream  of  blade  #1  (45.7  cm  downstream  of  the  turbulence  grid),  on  blade  #1  at  the  40%  axial 
chord  position,  76.2  cm  downstream  of  blade  #1,  and  the  insulated  pressure  transducer  box.  These 
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temperatures  were  displayed  in  the  upper  left  box.  The  calculated  air  density  was  also  displayed 
within  this  box.  Thermocouple  temperature  data  was  collected  with  J-type  thermocouples  using  a 
Kaye  Instruments  thermocouple  0°C  (32°F)  Ice  Point  Reference  and  the  HP  3852A  Data  Acquisi¬ 
tion/Control  Unit.  The  thermocouples  provided  the  current  temperature  for  the  other  data  collection 
programs,  including  the  Grid  Traverse  Data  Acquisition  VI. 

D.4  Velocity  Profiles  Across  The  Blade  Cassette  Inlet  and  Exit 

The  Grid  Traverse  Data  Acquisition  VI  was  used  to  verify  periodicity  and  uniformity  of  the 
blade  cassette  inlet  and  exit  flowfield.  The  velocity  and  turbulence  intensity  of  the  air  was  measured 
with  a  TSI  1210-T1.5  hot  wire.  The  VI  operated  the  Dantec  3-axis  traverse  to  move  the  hot  wire 
across  the  blade  row  inlet  or  exit  plane.  The  coordinate  system  for  the  traverse  is  shown  in  Figure 
139.  In  Figure  139  each  of  the  traverse  slots  is  numbered  1  through  4,  where  the  first  two  are 
for  measuring  the  inlet  flow  one  and  one-half  axial  chord  lengths  upstream  of  the  blade  row.  The 
centerline  position  of  these  1.27  cm  (0.5  inches)  wide  slots  relative  to  the  blades  is  17.78  cm  (7 
inches)  i  0.0127  cm  (0.005  inches)  and  8.89  cm  (3.5inches)  ±  0.0127  cm  (0.005  inches)  upstream 
of  the  blades  for  slots  #1  and  #2  respectively.  The  downstream  slots  #3  and  #4  are  the  same  distance 
away  from  the  blades  as  #2  and  #1,  respectively.  The  hot  wire  (symbolized  as  (b)  in  Figure  139)  was 
moved  upstream  of  the  blades  (a)  while  suspended  through  the  traverse  slot  moving  along  the  y- 
axis.  The  span  covered  by  the  hot  wire  was  6 1  cm  (24  inches)  which  included  the  passages  between 
blades  #3  through  #7.  The  hot  wire  was  moved  up  and  down  using  the  z-axis  to  get  velocity  and 
turbulence  intensity  profiles  for  different  heights  along  the  blades.  The  particular  slot  the  hot  wire 
was  in  was  changed  by  moving  the  z-axis  up  and  thus  lifting  the  probe  out  of  the  tunnel  and  then 
moving  the  x-axis  forward  or  back  to  the  next  desired  traverse  slot.  The  velocity  and  turbulence 
intensity  profiles  in  sections  3.3  and  3.4  were  taken  for  all  four  traverse  slots  at  Reynolds  numbers 
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of  50k,  100k,  and  200k.  The  Grid  Traverse  Data  Acquisition  program  was  divided  into  two  sections, 
the  left  portion  was  for  calibration  and  control  inputs,  while  the  grey  right-hand  portion  displayed 
the  velocity  and  turbulence  outputs. 


+Z  +X 


4 - 

Figure  139.  Tunnel  test  section  3-axis  coordinate  system  used  with  the  Dantec  traverse. 

D.4.1  Blade  Cassette  Velocity  Profile  Inputs 

The  program  required  the  atmospheric  pressure  and  dew  point  to  be  entered  into  the  upper 

left  frame  seen  in  Figure  140.  Other  atmospheric  conditions  such  as  temperature  and  density  were 

calculated  using  thermocouples  and  internal  calculations  based  upon  room  pressure.  The  hot  wire 

calibration  information  was  located  in  the  frame  below  the  atmospheric  data.  It  was  necessary  to 

enter  the  updated  hot  wire  calibration  information  in  this  frame.  The  number  of  data  samples  and 

the  frequency  of  their  collection  at  each  traverse  location  was  also  entered  here.  The  3 -axis  traverse 

was  controlled  through  the  lower  left  frame  seen  in  Figure  140.  The  number  of  steps  and  spacing 
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between  each  collected  data  point  across  the  pitch  axis  of  the  blades  was  entered  here.  The  same 
control  over  the  z-axis  was  maintained  through  this  frame.  Every  averaged  data  point  recorded  for 
each  y,z  location  was  displayed  on  the  right  hand  of  the  computer  screen  in  two  forms. 


Grid  Inverse  UatActj.vi ' 

'  ft 

eftfiPiiSfe  I'f . 

8  ,  PJ 

11 

1# 

at 

;  ‘  3; 

i* 

560 

9IJ: 

Bill 

iS. 

BiigSkX| 

m 

m 

B 

■ 

mmMm 

» 

B 

K1 

MS~~i 

m 

mi 

n 

■ 

m 

m 

a 

■ 

a 

■ 

m 

a 

■ 

;a  •••" '  i 

ns 

m 

IJP 

m 

m 

-a 

ai 

Bi 

P 

MjP| 

’.i  84- 

C  ■  82- 

an 

SB 

" 

RS 

FK* 

•  •  -  ..nie  tnrp  a  i  4-runs  tsM (J,  ♦>. -tl  0  cm 

•■.-v: dattsc  cotumki[iore:-:||RpM,i 

fl® 

m 

Ur  -  .  -  ■ 

□ 

i 

i 

| 

| 

1 

■’TTI" 

% 

B 

II 

Wjj 

-1 

■ 

:  V  TVaverse  .slol-.lfe  i  :  v  v  •• . ;  V::; "t :  - :  i :  e  r-: 

■ 

i 

i 

i 

1 

I 

I 

i 

te 

Bl 

il 

1 

IB 

1 

1 

■ 

:;v  fi 

m 

B 

m 

5E5pl?> 

; 

I® 

ilM 

■.  '£ 

*■ 

m 

-I 

■■ ‘i 

temperature  source' 

‘i)~I 

IT 

Si  s 

si 

l 

*  .» 

S  | 

■W 

M 

m 

N 

...  ■1; 

Keybq^rderrtryTern 

sou  eS31I 

iheee 

X  9- 

±#:i 

m 

■ii.fi 

m 

e 

a 

mm 

a 

E3 

B 

ip 

■HLiliffinB 

ietfe  Pg 

ibctij 

/Profile.'’ 

fiflr" 

km  MBPS 

_ 

* 

_J 

>5 

wm 

£ 

S' 

R 

i 

& 

■ 

■ 

■ 

■ 

m 

I 

Bd 

mi 

1  :,w- 

J*rj 

- 

£2 

9 

E8 

££ 

S 

m 

|rS; ;. 

SSEEh 

m 

1 

8 

n 

m 

51 

II 

S3 

S£ 

si- 

mSW  HImm 

s.w- 

n 

I 

■ 

i 

H 

■ 

1 

'  WMBI  - -  PCB  ■ 

t/vi 

■•■T 

O'  1 

Z  1 

4  .4 

B'l 

J  Z 

IS 

n 

s’  -,30  ¥i 

JHBHfflflffa  u&kbb  m 

&s&Ieise1 

§ii§ 

Turbulence  Mbo:F 

■ 

....... 

“  ";l" 

V.*.o. 

— -aBBs.- 

ihhhumi 

m 

fcj 

Jf 

Tjpji 

■ 

13 

iP 

11 

MB': 

«P 

■  ■■  1 

■JL, 

J 

1 

F, 

im. 

!  1 

3-AXIS  Traverse  intormation:  :  \  i-Sp^ri  fiOOrpm) Udob  oo"  1 

k  n-nnn. 

il 

■-V2 

1 

J* 

L 

w 

5*^ 

1  46  -IS 

J 

«<i 

L>.™ 

Delta  Z  (cm)  height  (cm)  ffZ  Stens  2  Step  #  NQTE^ 

av-j 

Ta(H)  |l  06195  i 

■-  a- 
gUj 

n 

in 

iill 

P) 

•  10  12  •  N 

% 

iM 

MR 

% 

il 

JA 

A 

26 

i!i 

fE-gJ  KM13  - 

DetoVM  Y-om'b  (or)  KVStgls  YSttm* 

BHHH 

Tur  Protife^  j 

WHWJ 

faa 

»w;v 

~y?? 

(EKSJ 

to 

jgpjijP 1 

1  Irrt  tinrie  (GO  VHr) 

:  Bm, 
ginwEl ; 

re- 

: 

n 

«son  raa  eh  Pi  i“®”'  5 

i1 

:* 

•;  J 

y.w 

Z*Qitectiemj  -Y-tMctionl  jews  efts eig^H 

— 

*; 

\ 

•; 

_Jj 

it 

!  \ 

J 

y 

■X- 

j:'\ 

£  •:  / 

it 

j\ 

■. 

,?4, 

rVj 

-5.. 

voll  slot  jo  1 

Kpip- 

$ 

M 

1 

Ml 

i 

5 

M 

il 

N 

*; 

-“IB 

!  . .  11-. 

i  '  -BW  .  ■  '  l 

rw 

wm 

■  ...: 

n 

A 

• 

a: 

■  fe 

MMzmwmumm 

m 

i"" : 

U 

i 

■ 

B 

n 

10  12  1 

M 

I 

m 

~ 

Si 

B 

<!.  J 

i-is 

*  »_R: 

IM 

Figure  140.  \felocity  and  turbulence  measurements  made  using  3-axis  traverse  and  hot  wire  probe 
program. 


D.4.2  Blade  Cassette  Velocity  Profile  Outputs 

The  outputs  of  the  Grid  Traverse  Data  Acquisition  VI  were  presented  in  color  to  the  user, 
but  are  printed  here  in  a  grayscale  format.  The  actual  velocity  data  taken  at  each  traverse  location 
was  displayed  in  near  real  time  in  the  upper  right-hand  comer  of  Figure  140.  Presentation  of  the 
complete  raw  data  collection  allowed  for  inspection  of  AC  noise  interference  or  unexpected  velocity 
fluctuations .  The  velocity  information  for  the  particular  x,y,z  location  was  averaged  and  added  to  the 

collection  for  the  particular  pass  in  the  y-axis  direction.  Each  y-axis  pass  was  ran  at  a  different  height 
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with  Figure  140  illustrating  5  passes,  each  5  cm  (1.97  inches)  apart  in  the  z-axis.  The  average  values 
were  presented  in  two  different  forms.  The  first  average  velocity  graph  was  a  color-scale  contour 
map.  Each  averaged  point  was  assigned  a  block  of  the  appropriate  color  depending  on  its  average 
velocity.  The  pitch  of  the  blades  spreads  across  the  x-axis  of  the  graph  (traverse  y-axis)  and  the 
different  passes  in  the  z-axis  move  up  the  maps  y-axis.  This  contour  mapping  eased  determination  of 
improper  periodicity  and  the  effects  of  height  within  the  tunnel  during  data  collection.  The  velocity 
data  was  also  graphed  on  a  regular  line  chart  just  below  the  velocity  contour  map  as  seen  in  Figure 
140.  Each  pass  along  the  pitch  axis  of  the  blades  was  given  a  different  symbol  and  line  style  to  ease 
in  differentiation.  For  example,  the  pass  along  the  center  of  the  blades  was  a  simple  line  with  no 
symbols,  while  the  pass  5  cm  (1.97  inches)  above  the  center  was  marked  by  solid-filled  squares. 
The  same  symbols  were  given  to  the  turbulence  intensity  line  chart  located  in  the  lower  right-hand 
comer  of  Figure  140.  A  color  contour  map  of  the  turbulence  intensity  data  was  also  presented.  The 
Reynolds  number  was  also  calculated  each  time  the  velocity  measured. 

The  center  of  the  display  shows  the  calculated  Reynolds  number.  Reynolds  number  was  based 
upon  axial  chord  length  and  inlet  velocity.  Also  displayed  is  the  Reynolds  number  based  upon  tme 
chord  length  and  exit  velocity.  Conversion  of  the  Reynolds  number  to  other  standards  is  shown  in 
Table  3.  As  with  the  Grid  Traverse  Data  Acquisition  VI  the  Pressure  Transducer  Calibration  VI  was 
visually  divided  into  two  separate  areas  for  the  data  entry  and  output. 

Table  3.  Reynolds  number  conversion 


Axial  Chord, 
Inlet  \blocity 

True  Chord, 
Exit  \felocity 

Tme  Chord, 
Inlet  \felocity 

Axial  Chord, 
Exit  \felocity 

50,000 

119,157 

81,500 

238,315 

476,630 

292,411 

326,000 

500,000 

1,191,574 

731,027 

815,000 

1,000,000 

2,383,147 

1,462,054 

1,630,000 
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D.5  Pressure  Transducer  Calibration 


Two  pressure  transducers  were  calibrated  simultaneously  using  the  Pressure  Transducer  Cali¬ 
bration  VI  with  the  Druck  DPI-600  Digital  Pressure  Indicator  calibrator  and  pressure  source.  The 
DPI-600  was  connected  to  the  total  pressure  side  of  the  differential  transducers  with  the  static  side 
open  to  room  air.  The  DPI-600  provided  the  pressure  to  the  transducer  as  well  as  the  resulting  refer¬ 
ence  pressure.  Using  the  HP  3852A  Data  Acquisition/Control  Unit,  the  millivolt  outputs  were  read 
from  the  transducers.  Each  transducer  output  was  read  30  times  and  then  result  averaged.  Each  cal¬ 
ibration  was  done  for  the  appropriate  range  to  be  measured  in  the  next  experiment,  with  pressures 
ranging  from  0  to  1.016  x  10~2  kg/cm2  (4  in-H20).  The  resulting  calibration  curve  coefficients 
were  displayed  in  the  right  hand  side  of  Figure  141.  These  coefficients  were  then  entered  into  the 
Wfcke  Traverse  Data  Acquisition  VI  and  other  programs  that  required  the  use  of  the  pressure  trans¬ 


ducers. 


I  it  pressure  transducer  calibration  -  a.vi " 


NOTE:  Use  home  baton  on  controller  to  read  through  portO  on  scanivalve. 

___  1:  blade  #4  rr 

#1  Transducer  to  read  jl}/  2:  grid  P-drop 

#3  Transducer  read:  IT! 1  3:  blade 

“--I  *  4  cassette  P-drop  -'  I 

5:  Pitot  P-drop 

#  of  P-Trans  Readings:  “ 

P-tap  Integration  time  (60  Hz)  ||16?| 


take  •  :••• 
.  DATA  I 


■ '  .  #1  mean  mV. |2.6031  j  .  V; 

■  '  #r Ptap' mil'KnJt  fc.59956  \  ;■ 


lidpe  f}- 1 59460i  V : 

.  ift  Ofdy intercept  : 


•  M^>»3:mdwm^:'|2.3405  j 
■  #3  Ptap  millhvofr:  ^33742  | 


'■■■.  #3  Pfop  pressure  taH2Cr.R6~497 1 j 
aentONA  t slope:  tf-V5946tjj 
• .  Qldy  intercept  ;j|-0  94S8(f  j 


Figure  141.  Pressure  transducer  calibration  software. 
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D.6  Wake  Traverse  Velocity  and  Pressure  Data  Collection 

To  measure  the  exittotal  pressure  and  velocity  downstream  ofthe  test  blade,  the  Wake  Traverse 
Data  Acquisition  VI  was  created  by  modifying  the  existing  Grid  Traverse  Data  Acquisition  VI.  The 
modifications  to  the  Grid  Traverse  Data  Acquisition  VI  included  removal  of  the  color  contour  plots 
and  addition  of  current  sample  and  mean  total  pressure  sample  plots.  The  program  inputs  were  still 
placed  on  the  left  (white)  portion  of  the  screen  with  the  outputs  on  the  right  hand  side  of  Figure  1 42 . 


Pitot  Iravvrsc  Data  Ac<i  -  with  Som.vi 


I  Header  File  (indicoton  [||c:\tAKEDATA\P^axit\smpoth\tB-1  \Smo-p1D0n3 

file  into:  fls'mooth:  Re  ■  IQDk,  Tu  -  ISi,  3  Apr  39  MS-  2{To  56  cm  ta  B8  cm  #3 
dote  column  tiffe:  |Re«1DDk 
Traverse  slot:  SO 


|  temperature  source. 


2  Keyboard  entiyTemp  (F) ptODO 
Ambient  Pressure  (psi)jtl  4273 


Channel  list  |?00-205  |  Ie7^3lj 
thermocouple:  |fl  1 


tsure  (psi)»A2m \  cflJ|  ?57m,  p  r~  25  |  I 

Tdew(F)g55.4 1  J  forupdalas  9  !  "  *  I 

Air  densrty 0<q/m3)  |71H~1  HW'cmp  I87  25  1  1 


«r-i  Probe  Type  E3  _ 

Oiannel  R(g  Q  deqC  g581  |  Mean  votoga  B980  I! 

cw^jhH2^  Rnooc-oalOn 

D  slope  |m  Rop|^g| 

|j  sample  rate  ll0000.C| _  #  of  samples/ch  1^0000  j 


Z-axi3  start 

J  Delta  Z  (cm)  height  (am)  tfZ-Steps  ZStepff  NOTE. 

■raSo]  pnN  |i|  fi  I  down?*  21.9' 

I  Delta  V  (cm)  Y-axis  (cm)  #  Y- Steps  Y  Stop  #  (55^6on) — j 

|l E3  m  ID  a 

yo> ds  offset  [33*07*1 


Figure  142.  Exit  wake  total  pressure  and  velocity  measurement  program. 


D.6.1  Wake  Traverse  Data  Collection  Inputs 

In  addition  to  the  atmospheric  conditions  and  hot  wire  calibration  information,  this  VI  also 
required  the  pressure  transducer  calibration  information.  The  pressure  transducer  calibration  infor- 
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mation  was  entered  in  the  bottom  center  frame  of  Figure  142.  Control  of  the  probe  movements  was 
maintained  in  the  lower  left  frame  where  the  spacing  between  subsequent  measurements  was  en¬ 
tered.  The  location  of  the  probes  was  recorded  to  the  output  along  with  the  pressure  and  velocity 
data. 

D.6.2  Wake  Traverse  Data  Collection  Outputs 

The  pressure,  velocity  and  turbulence  intensity  for  each  location  were  recorded  and  displayed 
on  the  screen.  The  last  sample  of  velocity  data  was  displayed  in  the  upper  right  portion  while  the 
total  pressure  sample  was  displayed  in  the  bottom  right  of  the  screen.  The  average  values  for  each 
location  for  the  velocity,  turbulence  intensity,  and  total  pressure  were  also  displayed.  The  entire  data 
collection  was  also  saved  to  a  file,  allowing  calculation  of  the  loss  coefficient.  The  transducers  were 
also  used  in  gathering  surface  pressure  coefficient  information  for  the  baseline  Pak-b  blades. 

D.7  Surface  Pressure  Coefficient  Determination 

Pressure  coefficient  profiles  for  the  blades  were  obtained  for  both  blade  #4  and  #6  simultane¬ 
ously  using  the  Scanivalve  Data  Collector  VI  shown  in  Figure  143.  This  VI  operates  the  48-channel 
scanivalves  and  records  the  pressure  results  from  two  pressure  transducers.  The  results  from  the 
transducers  were  then  converted  and  plotted  to  display  both  the  actual  pressure  readings  and  the 
normalized  coefficients.  Tunnel  velocity  was  determined  using  a  hot  wire  and  a  pitot-static  probe 
both  located  one  axial  chord  length  upstream  of  blade  #5.  Temperatures  were  obtained  using  J-type 
thermocouples  within  the  test  section.  This  VI  was  divided  into  two  categories;  inputs  on  the  upper 
and  left  white  sides,  outputs  on  the  lower  and  right  grey  sides. 

D.7.1  Surface  Pressure  Coefficient  Inputs 

Atmospheric  data  including  dew  point  and  barometric  pressure  was  entered  in  the  upper  left 
portion  while  temperatures  from  the  thermocouples  were  displayed  in  the  frame  just  below.  The 
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Figure  143.  Blade  surface  pressure  coefficient  data  collection  program. 
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hot  wire  calibration  information  was  entered  in  the  lower  left  frame.  This  information  was  used  to 
convert  the  hot  wire  voltage  into  a  velocity  and  Reynolds  number.  The  resulting  Reynolds  number 
and  velocity  were  displayed  just  above  the  hot  wire  information  frame.  The  pressure  transducer 
information  was  entered  in  the  upper  center  frame.  Transducer  information  required  for  operation 
includes  the  number  of  surface  taps  to  collect  data  from  as  well  as  the  number  of  data  points  to  take 
from  each  surface  port  for  a  good  mean  value.  The  slope  and  offset  for  each  transducer  from  the 
pressure  transducer  calibrator  VI  was  also  entered  in  this  frame. 

D.7.2  Surface  Pressure  Coefficient  Outputs 

The  software  records  a  series  of  data  points  from  zero  to  43 .  Points  zero  through  3  9  are  pressure 
tap  locations  on  blade  #4  and  #6.  The  40th  data  point  measures  Ptotai  —  Pstatic  from  the  pitot  static 
probe  which  was  then  used  to  calculate  the  tunnel  velocity  and  serves  as  a  double  check  with  the  hot 
wire.  The  41st  data  point  records  the  pressure  shift  due  to  thermal  changes  around  the  transducers 
by  recording  the  difference  between  Ptotai  readings  on  both  the  total  and  static  side  of  the  transducer. 
While  the  difference  was  zero,  any  thermal  shift  in  the  zero  offset  was  recorded.  This  value  was 
then  subtracted  from  all  of  the  previous  data  points  to  minimize  the  temperature  shift  error.  Two 
additional  measurements  were  taken  including  P0  - PwaU  (point  #42)  and  P0  -  PA  (point  #43).  Point 
#42  provides  the  difference  between  the  pitot  probe  static  pressure  and  the  tunnel  test  section  static 
pressure.  Point  #43  was  the  difference  in  pressure  between  in  the  test  section  inlet  total  pressure 
and  atmospheric  pressure  which  gave  the  operator  enough  information  to  calculate  the  actual  total 
pressure  within  the  tunnel  when  room  atmospheric  pressure  was  considered.  When  data  collection 
was  complete,  the  operator  was  asked  if  the  data  was  good  enough  to  save  to  disk.  The  output 
file  consists  of  surface  port  number  and  pressure  information  for  both  blades  #4  and  #6.  Pressure 
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information  for  the  blades  includes  the  recorded  values  in  inches-HjO,  and  Cp  based  upon  hot  wire 
velocity: 


fi  _  Ptotal  P surfaces 

V~  1  Tt2 

2  Pu  inlet 

where  p  is  determined  from  atmospheric  pressure  and  tunnel  temperature.  No  compressibility  cor¬ 
rections  were  used  for  the  density  due  to  the  very  low  Mach  numbers  of  the  test  section  airflow. 
The  inlet  velocity  used  in  equation  27  was  sampled  immediately  before  each  pressure  tap.  The  in¬ 
let  velocity  was  sampled  10,000  times  for  each  surface  port  and  then  averaged.  The  recorded  mean 
velocity  and  turbulence  intensity  for  each  measurement  were  displayed  in  the  graphs  on  the  lower 
right  of  the  screen  shown  in  Figure  143. 

The  grey  portion  of  Figure  143  contains  the  outputs  for  this  VI.  In  the  center  of  the  display  are 
the  most  recent  data  samples  for  each  blade  in  units  of  inches  of  water.  Inspection  of  these  values 
allows  determination  of  unacceptable  pressure  fluctuations  and  so  the  quality  of  the  data  can  be 
determined  for  each  data  point.  The  mean  of  these  values  were  plotted  for  each  pressure  port  in  the 
upper  right  of  the  display.  Temporary  dimensionless  Cp  values  were  also  plotted  on  the  screen  using 
equation  27. 


D.8  Boundary  Layer  Velocity  and  Turbulence  Measurements 

The  Boundary  Layer  Traverse  VI  operated  the  5.08  cm  (2  inch)  micro  traverse  which  held  a 
TSI  1279S-10A  subminiature  hot  film  probe.  The  built  in  accuracy  of  the  micro  traverse  allowed 
highly  detailed  measurements  within  the  boundary  layer,  down  to  within  0.0254  mm  (0.001  inches) 
of  the  wall.  This  VI  also  used  die  color  coded  division  of  inputs  and  outputs  to  ease  operation. 
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Figure  144. 


traverse  software. 
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D.8.1  Boundary  Layer  Traverse  Inputs 

The  inputs  for  the  Boundary  Layer  Traverse  VI  are  on  the  left,  or  white,  side  of  Figure  144. 
The  inputs  included  the  traverse  movement  controls  in  the  bottom  left  of  the  screen,  as  well  as  the 
hot  film  calibration  information  in  the  frame  just  above  it.  The  traverse  controls  allowed  the  user 
to  specify  a  step  distance  and  number  of  steps  to  complete,  \felocity,  turbulence,  and  temperature 
data  were  collected  at  every  step.  When  the  boundary  layer  profile  was  complete,  the  “Done  With 
Traverse”  button  was  pressed  and  the  output  file  examined  using  a  spreadsheet  program. 

D.8.2  Boundary  Layer  Traverse  Outputs 

The  velocity  and  turbulence  data  collected  was  displayed  in  the  right  hand  side  of  the  screen 
(Figure  144).  The  last  velocity  sample  collected  was  shown  at  the  top  right  of  the  screen,  while  the 
running  survey  was  displayed  in  the  large  box  below  it.  Both  the  velocity  and  turbulence  intensity 
profiles  were  put  on  the  same  graph  where  the  y-axis  of  the  graph  was  height  above  the  blade 
surface  and  the  x-axis  was  shared  by  the  velocity  (m/s)  and  the  turbulence  intensity.  The  displayed 
turbulence  intensity  was  the  actual  local  value  for  each  height,  defined  as: 

TU%=(i)  100  (28) 

where  o  is  the  rms  of  the  mean  velocity  U  at  height  y  above  the  surface.  In  addition  to  the  turbulence 
intensity,  the  length  scales  of  the  turbulence  were  also  calculated  for  the  inlet  and  exit  flow  from  the 
blade  cassette. 

D.9  Turbulent  Length  Scale  Determination  Using  Hot  Wire 

Turbulent  integral  length  scales  were  determined  for  the  Reynolds  numbers  and  turbulence 
intensities  under  which  the  Pak-B  blades  have  been  tested.  This  was  done  using  the  Auto-Correlation 
VI  and  the  single  hot  wire.  The  relative  humidity  and  atmospheric  pressure  were  entered  into  this 
program,  while  the  hot  wire  takes  100  complete  sets  of  data  samples.  The  100  sets  of  data  were 
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averaged  to  provide  power  spectrum  curve  data  in  addition  to  the  determined  turbulence  length 
scale. 

D.9.1  Required  Data  Inputs  for  Length  Scale  Software 

The  software  package  developed  by  Wright  Laboratory  for  the  calculation  of  turbulent  length 
scales  records  65536  velocity  data  points  at  a  rate  of  20,000  Hz  (Giese  [5 1])  from  a  single  hot  wire 
placed  in  the  wind  tunnel.  The  user  controlled  software  settings  allowed  the  maximum  frequency 
monitored  to  be  10,000  Hz.  The  high  sampling  rate  was  required  to  ensure  all  the  turbulent  fre¬ 
quencies  present  within  the  flowfield  were  captured  by  the  hot  wire  sensor  and  not  masked  by  a  low 
sampling  rate  (Halstead  and  Solomon  [24]).  Most  flowfield  frequency  fluctuations  taper  off  around 
5,000  Hz  (Hinze  [52]).  The  eneigy  of  the  various  frequencies  present  in  the  tunnel  flowfield  were 
presented  to  the  user  on  the  computer  monitor.  The  control  screen  of  the  Auto-Correlation  VI  is  pre¬ 
sented  in  Figure  145,  where  the  control  inputs  are  in  the  upper  left  white  sections  and  the  outputs  in 
the  right-hand  side  grey  sections.  Inputs  for  the  software  included  hot  wire  calibration  information 
and  the  atmospheric  properties  of  pressure  and  relative  humidity.  The  outputs  included  the  power 
spectrum  of  the  measured  flowfield,  the  auto-correlation  of  the  power  spectrum,  as  well  as  running 
strip  charts  of  velocity,  turbulence  level,  and  temperature.  The  hot  wire  velocity  values  were  used 
further  in  the  calculation  of  the  turbulent  length  scale. 

D.9.2  Calculation  Process  of  Length  Scale  Software 

The  Auto-Correlation  VI  used  the  recorded  velocity  data  to  find  the  root-mean-square  velocity 
and  average  velocity  U .  Other  parameters  used  were  the  time  step  (At  =  ~^)  equal  to  the  inverse 
of  the  sampling  rate  (SR)  and  the  frequency  resolution  equal  to  the  sampling  rate  divided  by  the 
number  of  data  points  (n)  taken  in  a  single  collection  (A/  =  §£).  The  eneigy  power  spectrum 
curve  was  then  created  using  a  Fourier  transform  of  the  rms  velocity  data.  The  power  spectrum 
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Figure  145.  Turbulent  length  scale  determination  software. 


207 


curve  was  a  single-sided  auto  spectrum  plot  of  G(f)  against  frequency  (/).  A  single  set  of  data 

contains  unsteady  fluctuations  which  could  result  in  a  jagged  curve  with  high  error.  To  reduce  the 

error  in  the  energy  power  spectrum  curve  calculation,  several  data  collections  were  required  using 

an  averaging  process.  The  averaging  process  may  be  done  in  real  time  by  dividing  the  values  in  the 

next  data  set  by  the  total  number  of  data  sets  collected  so  far  plus  one  (Bendat  and  Piersol  [53]): 

[running  average]  — -I - - —  [new  data  set]  (29) 

n  +  1  n  +  1 

and  adding  this  to  the  current  average  multiplied  by  the  ratio  A  typical  resulting  power  spec¬ 
trum  is  presented  in  Figure  146. 


D.9.3  Turbulence  Intensity  Verification  Using  Power  Spectrum  Curve 

A  verification  was  performed  to  ensure  accuracy  of  the  obtained  spectral  curve.  First  the  nor¬ 
malized  velocity  was  calculated  by  integrating  the  power  spectrum  curve  (Bruun  [31]): 

—  u  r°° 

<  =  2 -JQ  GU)df  (30) 

which  gives  a  result  in  units  of  .  Assuming  the  v!  and  v'  velocity  components  are  of  the  same 


using  other  methods. 


D.9.4  Integra]  Length  Scale  Determination  and  Program  Outputs 

The  auto-correlation  curve  (px(t)  versus  r)  was  created  using  an  inverse  Fourier  transform 
on  the  power  spectrum  curve.  This  required  the  time  step,  At,  and  the  frequency  resolution.  A/. 
Normalization  of  the  resulting  auto-correlation  curve  was  obtained  by  dividing  the  data  set  by  the 
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maximum  value  within  it.  The  integral  time  scale,  Tj,  was  determined  by  locating  the  zero  crossing 
of  the  curve  (t*  -  see  Figure  147).  The  length  scale  was  calculated  by  first  integrating  the  auto¬ 
correlation  curve  (Bruun  [31]): 


(33) 


and  then  multiplying  it  by  the  mean  velocity: 


L  =  TiU  (34) 

which  resulted  in  the  desired  integral  length  scale  (L).  Graphs  of  the  power  spectrum  and  auto¬ 
correlation  curves  may  be  produced  on  any  spreadsheet  program  using  the  output  data  file. 
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G(f)  [m  ^  2/s] 


n  Re=50k 


frequency  (Hz) 


Figure  146.  Power  spectrum  for  inlet  conditions  of  Re  =  50k  and  Tu  =  4%. 
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APPENDIX  E  -  2-D  Vane  Blade  Interaction  Code 


E.l  VBI  Software  Inputs 

The  VBI  2-D  code  was  capable  of  modeling  a  complete  turbine  stage  consisting  of  both  rotors 
and  stators.  While  the  capability  to  model  a  moving  blade  row  and  the  associated  wakes  was  not 
used,  the  inputs  for  both  sets  of  blades  are  shown  in  Table  4.  The  required  units  for  the  program  are 
also  given  in  the  table.  Values  for  user  defined  parameters  such  as  inlet  relative  Mach  number  and 
inlet  total  pressure  were  determined  using  the  17.8  cm  (7  inch)  axial  chord  length  Pak-B  blades  and 
the  environmental  conditions  of  the  experimental  facility.  All  of  the  initial  conditions  were  applied 
to  both  the  airfoil  wrapping  O  grid  and  the  passage  filling  H  grid. 

Table  4.  2-D  VBI  code  aerodynamic  inputs 


Parameter 

Stator 

Rotor 

Rotor  speed 

rpm 

Stage  equivalent  work 

Btu/lb 

stage  expansion  ratio 

ratio 

Inlet  relative  Mach  number 

user  defined 

user  defined 

Exit  relative  Mach  number 

user  defined 

user  defined 

Inlet  relative  flow  angle 

degrees 

degrees 

Exit  relative  flow  angle 

degrees 

degrees 

Inlet  total  temperature 

°R 

°R 

Inlet  total  pressure 

psi 

Corrected  flow 

lb/s 

Reynolds  number 

user  defined 

Inlet  turbulence  intensity 

% 

Reduced  frequency 

user  defined 

user  defined 

\hne  setting  angle 

degrees 

\hne-Blade  spacing 

inches 

E.2  The  k  —  Two-Equation  Model  for  Turbulent  Viscosity 

The  overlapping  nature  of  the  O  and  H  grids  used  in  the  VBI  code  presented  a  problem  for  typ¬ 
ical  turbulence  models  (e.g.,  Baldwin-Barth)  where  calculation  of  the  dimensionless  wall  distance 
y+  was  required.  The  turbulent  viscosity  fxt  for  the  VBI  code  was  modeled  without  the  calculation 
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of  y+  by  using  the  k  -  3ft  two-equation  model  (Kelecy  and  Delaney  [54]).  The  turbulent  viscosity 

Pt  was  calculated  by  a  multistep  process.  The  first  step  was  solving  for  the  transport  variable  9?: 

k2 

&  =  —  (35) 

ev 

and  then  using  this  to  find  the  turbulent  Reynolds  number: 

Re  y  =  — .  (36) 

Once  and  Re  t  were  calculated,  the  turbulent  viscosity  was  solved  for  using: 


Vt  =  C^pdt  (37) 

where 


,  _  1  -  exp(-AM  Re  f.) 

M  1  —  exp(— A€  Re  ^) 

The  constants  involved  with  equation  37  and  38  were  =  0.09,  A p  =  2.5  x  10-6 


(38) 

,  and  Ae=  0.2. 
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APPENDIX  F  -  Baseline  Performance  Loss  Coefficient  Estimation 


The  loss  coefficient  has  been  measured  extensively  by  several  researchers  with  an  effort  to 
create  a  universal  curve  fit  for  various  blade  shapes  with  different  turning  angles.  The  original  ex¬ 
perimental  work  of  Ainley  [35]  has  been  expanded  upon  and  re-iterated  through  the  texts  of  Horlock 
[39]  and  Mattingly  [20,55], 

F.1  Extrapolation  of  Losses  from  Previous  Experiments 

Ainley  [35]  first  published  experimental  loss  coefficient  data  for  turbine  blades  in  1948.  Mea¬ 
surements  of  2-D  profile  loss  and  total  loss  were  made  for  variations  of  the  T6  turbine  blade  profile . 
The  blade  shape  was  tested  at  different  inlet  angles,  axial  chord  to  pitch  ratios  (solidity,  a),  Reynolds 
numbers,  and  thickness  to  chord  ratios  provide  experimental  data  for  a  wide  range  of  practical  con¬ 
ditions.  Ainley’s  measurements  were  made  using  a  modified  T6  airfoil  base  profile.  The  T.6  was 
placed  in  a  cascade  and  tested  at  a  range  of  camber  angles  and  spacings.  Measurement  of  the  loss 
coefficient  (7)  for  the  T.6  turbine  blade  in  a  reaction  blade  configuration  was  7  =  0.04  and  0.03  for 
attached  flow  and  inlet  Reynolds  numbers  of  50k  and  100k,  respectively,  and  for  an  incidence  of  0° 
(Ainley  [35]).  The  reaction  blade  results  of  Ainley  [35]  agree  very  well  with  the  measured  Pak-B 
loss  coefficient  of  7  =  0.038  for  an  inlet  Reynolds  number  of  100k. 

F.2  Estimation  of  Losses  Using  Curve  Fits 

Ainley  and  Mathieson  [36]  continued  to  experiment  with  turbine  blades  providing  a  series  of 
loss  coefficient  estimation  curves  for  a  range  of  solidities  and  gas  flow  exit  angles.  For  the  Pak-B 
design  exit  angle  of  60°,  the  loss  coefficient  was  0.115  for  impulse  blade  relations  (inlet  angle,  f31 
=  -exit  angle,  012)  and  0.023  for  nozzle  blade  relations  (p1  =  0°).  The  experiments  for  these  curves 
were  completed  at  an  approximate  inlet  Reynolds  number  of  123k  with  an  incidence  of  0°.  Ainley 
and  Mathieson  [36]  also  provided  a  curve  fit  for  their  data  so  the  profile  loss  coefficient  other  blade 
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shapes  could  be  estimated: 


7  = 


7(^=0)  + 


2 


7(/3l=-a2)  ~  7(^=0) 


(39) 


where  t/c  is  the  thickness  to  chord  ratio.  Using  equation  39  for  the  Pak-B  design  angles  and  thickness 
to  chord  ratio,  the  loss  coefficient  was  estimated  to  be  0.069.  The  loss  coefficient  can  be  adjusted 
to  account  for  changes  in  Reynolds  number  using  a  curve  fit.  The  Reynolds  number  corrected  loss 
coefficients  are  0.07  and  0.12  at  50k  and  100k,  respectively.  The  difference  in  the  estimation  of  loss 
coefficient  using  equation  39  and  experiment  can  be  attributed  to  the  advised  restriction  that  the 
equation  only  be  used  for  0. 15  <  t/c  <  0.25  (Ainley  and  Mathieson  [36])  where  the  Pak-B  has  a  t/c 
=  0.134  and  the  inability  to  account  for  differences  in  aspect  ratio. 

Horlock  [37, 38]  added  to  the  work  of  Ainley  [35]  and  Ainley  and  Mathieson  [36]  by  incor¬ 
porating  the  unpublished  data  of  C.  R.  Soderbeig.  Horlock  [37]  presented  the  previously  unpub¬ 
lished  data  next  to  the  work  of  the  published  available  data.  Soderbeig ’s  data  had  loss  coefficients 
of  0.082,  0.07  and  0.059  for  Reynolds  numbers  of  50k,  100k  and  200k,  respectively,  using  blades 
with  a  higher  thickness  to  chord  ratio.  Horlock  [39]  also  presented  a  curve  fit  for  Soderbeig ’s  data: 


7  =  °.°25[l  +  (A)2 


(40) 


where  e  is  the  deflection  angle  in  degrees.  Using  the  Pak-B  total  turning  (deflection)  angle  of  95°, 
the  loss  coefficient  was  predicted  to  be  0.053  with  equation  40. 

Craig  and  Cox  [40]  published  a  series  of  curve  fits  to  account  for  several  aspects  of  blade 
geometry.  Flow  angles,  blade  spacing,  contraction  ratio,  and  trailing  edge  thickness  were  accounted 
for  in  the  figures  published  by  Craig  and  Cox  [40] .  Using  the  figures  with  the  geometry  data  for  the 
Pak-B  blade,  loss  coefficients  were  estimated  to  be  0.052,  0.043  and  0.034  for  Reynolds  numbers 
of  50k,  100k,  and  200k.  These  values  match  the  experimental  data  obtained  for  the  Pak-B  blade 
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when  flow  exits  the  blade  passage  attached  to  the  suction  surface.  For  example,  for  small  separation 
zones  at  a  freestream  turbulence  level  of  4%,  the  estimated  losses  are  only  10.3%,  10.3%,  and  8.1% 
off  from  the  experimental  losses  for  Reynolds  numbers  of  50k,  100k,  and  200k,  respectively. 
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APPENDIX  G  -  Blade  Surface  Curvature  Distribution 


Wilson  and  Korkianitis  [56]  have  shown  that  the  boundary  layer  is  directly  affected  by  the 
continuity  of  the  slopes  that  make  up  the  surface  of  the  blade.  The  continuity  of  the  blade  surface 
slope  can  be  examined  using  the  radius  of  surface  curvature  and  the  slope  of  the  radius  of  curvature. 

The  radius  of  surface  curvature  can  be  determined  using  (Wilson  and  Korkianitis  [56]): 

1  ( d2y/dx 2) 

r  =  /  x3/2  (41) 

(l  +  (dy/dx) 2) 

where  dy  and  dx  are  the  differences  between  sequential  y  and  x  coordinates,  respectively,  using  a 
three  point  differencing  scheme.  The  resulting  values  of  1/r  for  the  suction  and  pressure  surfaces 
are  shown  in  Figure  148.  The  curves  in  Figure  148  show  several  discontinuities  due  to  limitations 
in  the  original  coordinate  data.  These  limitations  prevented  useful  insight  into  what  effects  surface 
curvature  may  have  had  on  the  baseline  separation  position.  The  slope  of  the  radius  of  curvature 
could  not  be  calculated  due  to  the  same  limitations. 
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Figure  148.  Radius  of  curvature  plots  for  Pak-B  blade. 
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APPENDIX  H  -  Hot  Wire  Control  Equations 


H.l  Hot  Wire  Physical  Characteristics 

The  hot  wire  used  in  ffeestream  velocity  and  turbulence  measurements  was  a  single  wire  TSI 
1210-T1.5.  The  wire  element  had  a  diameter  of  3.87  x  10'3  mm  (1.52  x  1(T 4  inches)  and  a  length 
of  1.27  mm  (0.05  inches).  The  hot  wire  was  powered  by  an  IFA-100  Intelligent  Flow  Analyzer 
which  also  provided  a  voltage  output  which  was  calibrated  to  provide  the  flow  velocity. 

H.2  Velocity  Determination  Equations 

The  velocity  was  calculated  through  a  series  of  equations  which  accounted  for  changes  in  the 
flow  condition  including  temperature,  pressure,  and  humidity.  The  temperature  within  the  tunnel 
was  read  through  thermocouples  while  the  current  room  pressure  was  entered  into  the  computer 
through  the  keyboard.  The  humidity  was  obtained  from  an  Omega  Engineering,  Inc.  RH-71A 
Digital  Hygro-thermometer  and  entered  into  the  computer.  The  following  equations  were  used  for 
every  data  point  taken  by  the  hot  wire.  The  control  room  computer  processor  speed  was  capable 
of  providing  all  the  velocity  and  turbulence  level  information  in  near  real  time.  To  calculate  the  air 
velocity,  several  factors  were  first  calculated  including  viscosity,  density,  and  the  current  molecular 
weight  of  air. 

H.2.1  Atmospheric  Molecular  Weight 

The  humidity  was  entered  into  the  computer  in  terms  of  its  dew  point  in  degrees  Fahrenheit. 
This  was  quickly  converted  to  Celsius  and  the  following  exponent  calculated  (Reid,  et  al.  [57]): 

Z  =  a$  +  Tp  ■  [a\  +  Tq  ■  [a,2  +  Tp  •  (03  +  To  •  0.4)]]  (42) 

where  To  is  the  dewpoint  in  Celsius  and  ao  through  a4  are  constants.  The  humidity  mass  fraction 
is  defined  as: 
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TUf  =  10Z. 


(43) 


The  humidity  corrected  molecular  weight  of  air  is  therefore: 


-  r 


1 


\  (i-TTC/)  ,  mt 
[  Mwa ir  '  Mwhio  j 

where  Mwair  is  the  molecular  weight  of  dry  air  and  Mwh20  is  the  molecular  weight  of  water. 


(44) 


H.2.2  Atmospheric  Constants 


The  mole  fraction  of  water  in  the  atmosphere  was  calculated  with: 

X  =  (  MWair  A _ mf 


\Mwji20  J 

where  rrif  is  the  previously  calculated  mass  fraction.  The  mole  fraction  of  air  is  simply: 


1—1—  ( Mmtutn  \  .  m  , 
\  Mwair.  )  ml\ 


(45) 


Xa  =  l-Xv. 


The  gas  constant  was  computed  using  the  corrected  molecular  weight  as: 


(46) 


which  is  in  units  of  J/(kg  ■  °K). 


8314.3 

XIw 


(47) 


To  calculate  the  specific  heat  at  constant  pressure  two  initial  variables  were  calculated: 


Cpa  =  bo  +  Tk  ■  [&i  +  Tk  ■  (p2  +  Tk  ■  63)] ,  (48) 

where  TK  was  the  tunnel  flow  recovery  temperature  in  Kelvin  and  60  through  63  are  constants.  The 
other  variable  was: 


CpV  =  1797.5  +  0.23  •  Tk- 


(49) 
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The  specific  heat  was  then  calculated  using' 


Cp  =  mfcpv  +  (l-mf)-cpa  (50) 

The  ratio  of  specific  heats  was  determined  from: 


7  = 


Op  —  R 


and  the  speed  of  sound  was  calculated  using: 


(51) 


a  —  VlR-  Tk ■ 


(52) 


H.2.3  Atmospheric  Properties 

To  calculate  the  viscosity  of  the  tunnel  airflow,  three  coefficients  were  first  computed.  The 
first  was: 


Hua  —  °o  +  Tk  (ci  +  Tk  ■  C2)  (53) 

where  co  through  c2  are  constants  and  TK  is  the  tunnel  temperature  in  degrees  Kelvin.  The  second 
was  (Reid,  et  al.  [57]): 


where  nuv 


was  a  constant.  The  final  coefficient  was: 


(54) 


4>b  =  4>a- 


{ f^uv\  {  \ 

UuJ  '  \MWH20)  ‘ 


The  viscosity  was  then  computed  using  (Reid,  et  al.  [57]): 


(55) 


+ 


-^a  *  Muq 

Xa  +  Xy  ■  (f>a  XV  +  Xa  ■  <f)b 


■Mu 


(56) 
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where  Xv  and  Xa  were  the  mole  fractions  calculated  previously. 

The  thermal  conductivity  was  determined  from  four  coefficients  based  upon  the  tunnel  test 
section  temperature  in  degrees  Kelvin.  The  first  coefficient  was  calculated  using: 


ka  —  do  +  Tk  ( d\  4-  Tk  ■  (I2) 
where  do  through  are  constants.  The  second  coefficient  was: 


(57) 


kv  =  e 0  +  Tk  (ei  +  Tk  ■  +  Tk  •  e 3))  (58) 

where  eo  through  e$  are  constants.  The  third  and  fourth  coefficients  were  based  upon  ka  and  kv. 
The  third  was  (Reid,  et  al.  [57]): 


The  final  coefficient  was: 


(59) 


V>21  — 


1  + 


( Mmua  \ 

\  MWair  ) 


0.25 


1  2 


The  thermal  conductivity  was  then  calculated  as  (Reid,  et  al.  [57]): 


(60) 


k  = 


Xa  *  ka 


+ 


Xy  *  ky 


Xa  +  Xv  •  2  Xv  +  Xa  ■  1p21 

which  allowed  quick  determination  of  the  Prandtl  number: 


(61) 


P  -  c2-tL 

r  k 

which  used  the  previously  calculated  specific  heat  at  constant  pressure. 


(62) 


222 


The  density  of  the  air  within  the  tunnel  was  simply  calculated  using: 


P  R'Tk  (63) 

where  P  was  the  atmospheric  pressure  of  the  wind  tunnel  facility  converted  to  metric  units.  No 
corrections  were  made  for  changes  in  density  with  increasing  speed  or  blockage  due  to  the  low 
velocities  used  during  testing. 

H.2.4  Velocity  Determination  Using  Atmospheric  Properties  and  Hot  Wire  Voltage 

The  velocity  of  the  air  within  the  test  section  was  calculated  using  the  previously  determined 
atmospheric  properties  of  density,  thermal  conductivity,  and  Prandtl  number.  The  calculated  flow 
recovery  and  hot  wire  gas  temperatures  were  also  used  along  with  some  hot  wire  physical  charac¬ 
teristics.  The  first  of  four  coefficients  calculated  was: 


Fi  = 


P_du 

P 


(64) 


where  dw  is  the  hot  wire  diameter.  The  second  coefficient  was: 


f2=  y2'Qw 


(«w  +  10)2  (65) 

where  Qw  was  the  hot  wire  probe  operating  resistance  and  V  was  the  measured  voltage.  The  third 
factor  was: 


F3 


(  Tm  \ 


Tr 


(66) 


where  Tm  was  the  gas  temperature  determined  as  an  average  of  the  flow  recovery  temperature  Tg 
and  the  probe  temperature  Tw.  The  exponent  Tr  was  the  temperature  overheat  ratio  used  by  the 
IFA-100.  The  final  coefficient  was: 
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F4  =  7T  •  k  ■  lw  •  (Tw  -  Tg)  (67) 

where  lw  was  the  length  of  the  hot  wire  and  7r  is  3.14159.  The  temperature  and  humidity  corrected 


velocity  was  then  calculated  first  using: 


which  was  then  used  in: 


1 


(69) 


where  Reexp  was  the  exponent  used  in  the  calibration  of  the  hot  wire.  The  intercept  for  the  hot  wire 
calibration  is  Cint  and  the  slope  is  Dsio. 

To  determine  the  flow  stagnation  temperature  the  mean  velocity  U  was  first  calculated  and 
then  the  following  equation  solved: 


Tg  =TW  ~ 


0.85  •  if 

2  •  Cp 


(70) 


which  was  then  compared  to  the  initially  used  Tg.  These  equations  were  all  then  recalculated  to 
provide  a  more  accurate  velocity  until  the  previously  used  and  the  current  calculation  for  Tg  were 
within  0.005  Kelvin. 


H.3  Turbulence  Determination  Equations 

Turbulence  intensity  was  determined  after  the  velocity  had  been  accurately  resolved.  The  col¬ 
lection  of  velocity  data  involved  sampling  the  voltage  of  the  hot  wire  40,000  times  at  a  rate  of  1 0,000 
samples  per  second.  The  velocity  was  calculated  for  each  of  the  40,000  data  points.  With  the  40,000 
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velocity  values,  a  mean  velocity  was  calculated  as: 

E  ?£4 


V  =  —  -  (71) 

where  Ux  was  the  individual  velocity  data  point  stored  in  an  array  Also  calculated  was  the  deviation, 
a,  from  the  mean  and  the  division  of  these  two  provided  the  turbulence  intensity.  The  turbulence 
intensity  was  multiplied  by  100  to  give  the  resulting  value  in  terms  of  percent.  The  calculation  of 
turbulence  was: 


Tu  = 


•100. 


(72) 
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APPENDIX  I  -  Uncertainty  Analysis 


1.1  Surface  Static  Pressure  Cp  Measurements 

Using  the  method  described  by  Kline  and  McClintok  [33],  the  first  step  in  determining  the 
uncertainty  is  the  calculation  of  the  partial  derivatives.  The  surface  static  pressure  measurement 
was  calculated  using: 

Cp  =  P°’in^r  Ps'x  (73) 

hpUinlet 

where  Potiniet  and  U inlet  were  the  total  pressure  and  velocity,  respectively,  measured  upstream  of 
the  blade  cassette.  This  equation  can  be  re-written  to  account  for  the  calculation  of  density  as: 


CV: 


2  (A P)  RT 


(74) 


PAUhet 

where  A P  was  the  pressure  difference  measured  across  the  transducer.  The  density  term  was  broken 
into  its  components  where  Pa  was  the  atmospheric  pressure  within  the  wind  tunnel  building,  R  is 
the  gas  constant  and  T  was  the  tunnel  test  section  temperature. 

There  are  four  partial  derivatives  that  need  to  be  calculated  from  equation  74,  including  mea¬ 
surements  of  A P ,  T,  Pa,  and  Ihnict  -  There  is  no  need  to  calculate  the  partial  derivative  with  respect 

to  R  since  it  is  a  constant.  The  four  partial  derivatives  are: 

dCp  2  RT 


dAP  pAulla 


(75) 


dCp 

dT 


2  (AP)  R 
PAUiraet 


dCp  _  -2  (AP)  RT 
9  Pa  p2TJ2 

rAU  inlet 


dCp  _  -4  (AP)  RT 
dU  PaTJ 3 

P'4  U  inlet 
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(76) 


(77) 


(78) 


The  second  step  in  determining  the  uncertainty  is  to  calculate  the  square  of  each  of  the  partial 
derivatives  and  divide  each  by  the  original  expression  squared  (equation  73).  Each  of  the  partial 
derivatives  after  squaring  and  dividing  becomes: 


fdC£_\2 
\dAP  ) 

i 

cv 

“  A  P2 
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The  uncertainty  in  the  surface  pressure  coefficient  equation  can  now  be  determined  using: 


r  .  .  o  i  1/2 

-P  =  Cp.[(^)2+(f)2+(^)  +(2.^)2]  (S3, 

where  zn  was  the  uncertainty  of  each  measurement  made  based  upon  the  equipment  used. 

1.2  Loss  Coefficient  Measurements 

Since  the  loss  coefficient  and  pressure  coefficient  use  the  same  equation  but  with  different 
symbols,  the  uncertainty  analysis  for  the  loss  coefficient  was  determined  using  the  method  described 
previously  for  the  surface  pressure  coefficient. 
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1.3  'Vfelocity  Measurements 


The  velocity  from  the  hot  wire  was  calibrated  using  a  water  manometer.  The  velocity  can  be 
determined  from  the  manometer  through: 


rr  1 2  (AP)  ~RT 

<84) 

where  A  P  was  the  difference  measured  between  the  total  and  static  ports  of  a  pitot  tube.  Following 
the  same  steps  as  done  with  the  surface  pressure  coefficient,  the  uncertainty  for  the  velocity  can  be 
written  as: 


wu  =  U  ■ 


(85) 


1.4  Wake  Momentum  Deficit  Measurements 


The  uncertainty  estimation  for  the  wake  momentum  deficit  was  calculated  using  the  same 
method  described  previously.  The  momentum  deficit  was  defined  as: 

(86) 


where  the  three  partial  derivatives  were: 

a© 

dUw 


dy  _  Uwdy 
U0 o  ui 


(87) 


d®  _  —Uwdy_  U^dy 
9Uoo  Ul  +  U'i 


(88) 


a©  =  -uw  ul 

ddy  t/oo  Ul’ 

The  next  steps  were  then  completed  using  a  spreadsheet  program.  Each  of  the 


(89) 

partial  deriva¬ 


tives  were  calculated  and  then  squared.  The  square  of  the  three  uncertainties  (ujf-  ,  ufj  ,  uj^y)  was 

also  calculated.  The  spreadsheet  then  multiplied  the  square  of  the  partial  derivative  with  its  respec- 
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tive  uncertainty  squared.  These  products  were  then  added  and  the  square  root  taken.  This  can  be 
written  as: 


+ 


1/2 


(90) 
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APPENDIX  J  -  Rejected  Separation  Control  Methods 


J.l  Separation  Control  Methods  Unsuitable  for  LPT  Use 

Several  methods  initially  considered  for  this  project  are  presented  in  Figures  149  through  155. 
These  methods  have  been  documented  in  research  for  the  reduction  of  separated  flow  and  reduced 
drag.  These  methods  were  rejected;  however,  since  the  reported  findings  have  shown  various  unde¬ 
sirable  qualities. 

J.1.1  Triangular  Plow  Vortex  Generator  Device 

The  plow  is  a  submerged  vortex  generator  (see  Figure  149a)  which  has  been  suggested  by 
many  authors  as  a  means  for  reducing  separation.  Chang  [43]  documents  this  design  which  was  to 
be  used  in  a  compressor  cascade  to  alleviate  the  shock-induced  boundary  layer  separation  (Austin 
[58]).  Time  constraints  prevented  Austin  [58]  from  completing  that  testing  and  it  was  not  until  1996 
that  Gamerdinger  and  Shreeve  [59]  managed  to  complete  the  work.  The  unexpected  result  for  the 
Mach  1.4  flow  indicated  an  improved  turning  angle  of  0.94  degrees,  but  also  an  increase  in  flow 
losses  by  approximately  8%. 


a)  Plow  b)  Semi-Circular  c)  Simple  Vane 

Figure  149.  Plow,  Semi-Circle  and  Vane  concepts  eliminated  by  available  data. 


230 


J.1.2  Semi-Circle  Vortex  Generator  Device 


The  semi-circle  is  constructed  from  a  solid  rod  cut  in  half  longitudinally  and  attached  flat  side 
down  to  the  surface  which  suffers  separation  problems  (see  Figure  149b).  These  cut  rods  were  then 
placed  in  either  a  parallel  or  V  pattern.  Rao  and  Kariya  [60]  conducted  semi-circle  experiments 
with  different  device  sweep  angles,  spacing  and  height  on  a  flat  plate  but  with  results  poorer  than 
the  other  designs  investigated  by  the  same  authors  which  included  typical  vane  vortex  generation 
devices.  The  unfavorable  pressure  recovery  results  found  by  Rao  and  Kariya  [60]  indicate  this 
design  was  a  poor  candidate  for  the  conditions  found  in  a  LPT 

J.1.3  Simple  Vane  Vortex  Generator  Device 

The  simple  vane  was  another  concept  considered  by  Rao  and  Kariya  [60]  that  consisted  of  a 
simple  vertical  triangle  with  height  equal  to  the  boundary  layer  thickness.  This  design  was  consid¬ 
ered  as  a  baseline  case  for  comparison  and  illustrated  the  capability  for  retaining  attached  flow  over 
a  flat  plate  (see  Figure  149c).  The  vanes  also  showed  a  substantial  reduction  in  the  total  pressure 
loss,  but  with  a  highly  distorted  profile  for  Cp  over  the  baseline  case.  Sweep  angles  tested  were  60° 
and  70°  relative  to  the  flow  velocity.  Considering  the  height  of  this  design,  a  substantially  higher 
drag  penalty  would  result  from  its  use  compared  to  a  submerged  vortex  generator  design. 

J.1.4  Passive  Cavity  Using  Adverse  Pressure  Gradient 

McCormick  [6 1  ]  examined  the  passive  cavity  method  on  a  flat  plate  as  a  means  for  reducing  the 
large  separation  bubble  due  to  shock  interaction  at  speeds  of  Mach  1.56  to  1.65.  The  passive  cavity 
utilizes  the  adverse  pressure  gradient  behind  the  separation  location  to  provide  some  flow  through 
a  cavity  within  the  plate  or  airfoil  and  into  the  boundary  layer  ahead  of  the  separation  line.  Results 
indicated  the  benefit  of  reducing  the  total  pressure  loss  through  the  shock  system,  but  significantly 
increasing  the  boundary  layer  losses  downstream.  Possible  implementation  of  the  passive  cavity  is 
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displayed  in  Figure  150a.  This  method  was  again  attempted  by  Lin  et  al.  [18]  with  aconvex  curved 
surface,  simulating  one  side  of  an  airfoil  with  turning.  Results  indicated  that  the  passive  porous 
surface  increased  the  reattachment  distance  instead  of  reducing  it.  For  the  low  Reynolds  numbers 
of  the  LPT,  this  method  was  considered  unsuitable. 
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blade 
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a)  Passive  Cavity  (side  view) 

b)  Passive  Injection  (side  view) 

Figure  150.  Porous  suction  surface  methods. 

J.1.5  Passive  Injection  Using  Adverse  Pressure  Gradient 

The  passive  injection  method  uses  a  porous  wall  similar  to  the  passive  cavity  method  described 
by  McCormick  [61],  but  with  the  air  now  only  allowed  to  flow  through  a  narrow  slot  ahead  of  the 
separation  line.  This  provides  a  powered  jet  effect  in  the  downstream  flow  direction  as  seen  in 
Figure  150b.  This  method  was  rejected  because  Lin  et  al.  [18]  discovered  that  passive  injection 
does  not  provide  any  reduction  in  separation  losses  over  the  unmodified  smooth  surface. 

J.1.6  ±45°  Small-Scale  Riblets  Located  On  the  Suction  Surface 

The  ±45°  riblet  technique  is  similar  to  conventional  riblets  but  with  the  ribs  in  a  ±45°  pattern 
in  the  area  where  separation  occurs.  Attempted  by  Lin  et  al.  [18]  over  a  convex  surface,  the  results 
for  a  variety  of  flow  conditions  always  increased  the  reattachment  distance.  The  reduced  pressure 
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recovery  of  the  model  with  this  modified  surface  (see  Figure  151a)  indicated  increased  losses  when 
compared  to  the  baseline  case. 

J.1.7  Leading  Edge  Vane  Vortex  Generators 

Chang  [43]  designed  the  leading  edge  vane  vortex  generators  to  extend  into  the  freestream 
airflow  at  the  front  of  the  airfoil  to  provide  trailing  longitudinal  vortices.  Chang  [43]  theorized 
that  the  vortices  generated  would  be  of  sufficient  strength  to  maintain  attached  flow  over  the  entire 
surface.  The  application  of  this  design  (see  Figure  151b)  would  create  several  problems  in  turbine 
engine  applications.  Addition  of  the  leading  edge  vane  vortex  generators  would  effectively  decrease 
the  spacing  between  the  blade  rows.  In  addition  to  the  spacing  and  heat  transfer  problems  involved, 
the  lack  of  data  on  this  concept  means  that  implementation  would  have  unpredictable  results  with 
the  Pak-B  blade. 

J.1.8  Wheeler  Channel  Groove  Cut  Into  Surface 

The  Wheeler  channel  is  a  V-shaped  groove  cut  into  the  surface  (see  Figure  151c).  The  air 
within  the  boundary  layer  spills  into  the  channel  and  vortices  are  thus  formed.  Wheeler  [62]  suggests 
using  multiple  rows  of  these  channels  if  necessary  to  maintain  attached  flow  all  the  way  to  the 
trailing  edge.  Implementation  of  the  complex  Wheeler  channel  design  involves  cutting  rows  of 
channels  into  the  turbine  blade  which  presents  a  formidable  machining  cost.  This  design  concept 
was  eliminated  due  to  the  complexity  and  cost  of  machining  the  wheeler  channels. 

J.1.9  Afterbody  Groove  Cut  Into  Surface 

The  afterbody  groove  is  a  cut-out  V-shaped  notch  with  a  valley  angle  of  90°  placed  over  the 
trailing  portion  of  the  surface  where  separation  normally  occurs.  Tested  by  Howard  et  al.  [63]  on 
an  axisymmetric  bluff  body  for  Reynolds  numbers  (based  on  diameter)  from  20  x  103  to  200  x 
103,  these  grooves  reduced  drag  by  16%  to  33%  by  causing  the  flow  within  the  grooves  to  remain 
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Figure  151.  Rejected  methods  with  difficult  application. 
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attached.  Howard  et  al.  [63]  speculated  that  the  groove  provides  a  pumping  action  for  the  three- 
dimensional  separated  flow  regions  similar  to  a  fluid  ejector.  The  lack  of  available  documentation 
for  this  experiment  does  not  support  application  to  a  LPT  blade.  The  more  extensively  tested  V- 
grooves  of  Lin  et  al.  [18]  discussed  in  Chapter  9  are  similar  to  the  afterbody  groove  and  have  been 
chosen  for  implementation  on  the  Pak-B  blade  since  they  provide  the  same  effect  as  the  afterbody 
groove  but  had  better  documentation  on  the  implementation  of  the  technique. 

J.1.10  Small-Scale  V-Groove  Riblets 

Riblets  are  small  scale  V-grooves  (see  Figure  152a)  which  have  been  extensively  investigated 
by  Whlsh  [64-66]  and  by  Bacher  and  Smith  [67,68].  The  riblet  sizes  are  determined  by  the  non- 
dimensionalized  height  and  spacing  based  upon  law  of  the  wall  coordinates  using  the  skin  friction 
coefficient  C /  and  kinematic  viscosity  v.  Nondimensional  height  of  the  riblets  is  defined  as: 


where  h  is  the  physical  height.  This  allows  easy  graphing  of  results  on  a  h+  vs.  s+  plot  for  a  variety 
of  test  cases.  The  nondimensional  spacing  between  riblet  peaks  is  similarly  defined  as: 


Walsh  [64]  has  attained  drag  reductions  for  a  flat  plate  with  a  turbulent  boundary  layer  of  up  to 
8%  for  h+  «  10  and  s+  «  15  while  Gallagher  and  Thomas  [69]  found  no  drag  reduction  and  only  a 
local  20%  reduction  in  skin  friction.  The  scope  of  possible  drag  reduction  using  riblets  is  narrowed 
by  the  results  of  Choi  et  al.  [70]  who  found  that  for  laminar  flows  no  net  drag  reduction  is  possible 
with  riblets. 

Many  attempts  have  been  made  to  use  the  apparent  drag  reduction  and  flow  attachment  quali¬ 
ties  of  the  riblets.  These  include  the  works  of  Wreck  [71]  on  cylinders  and  airfoils,  Reagan  [72]  on 
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Figure  152.  Other  submerged  design  modifications  considered  and  rejected. 
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diffusers  and  Rothenflue  [73,74]  on  concave  surfaces.  Wieck  [71]  demonstrated  a  separation  delay 
from  2%  to  1 1%  with  riblets  for  an  airfoil  in  an  adverse  pressure  gradient,  depending  on  freestream 
velocity.  Rothenflue  [73]  used  riblets  in  a  compressor  cascade  with  the  riblets  giving  the  maximum 
benefit  on  the  pressure  surface.  Suction  surface  riblet  placement  results  were  summarized  as  in¬ 
creasing  the  flow  turning  angle  by  1%.  Rothenflue  [74, 75]  also  showed  a  delay  in  transitioning  due 
to  Gortler  vortices  forming  above  the  riblets  for  a  concave  surface. 

Application  of  the  riblet  concept  to  turbomachinery  has  been  attempted  by  Chen  et  al.  [14]. 
Using  the  ideal  riblet  parameters  h+  <  30  and  s+  <  30  as  determined  by  Wblsh  [64],  a  10%  loss 
reduction  benefit  for  a  compressor  cascade  was  found  by  Chen  et  al.  [14]  when  the  riblets  were 
applied  to  the  pressure  surface  and  not  the  suction  surface.  When  riblets  were  applied  to  the  suc¬ 
tion  surface  Chen  et  al.  [14]  found  no  drag  reduction  or  delay  in  separation  from  the  suction  side. 
While  small-scale  riblets  have  shown  reductions  in  drag  for  flat  plate  flows,  the  extreme  turning  and 
pressure  gradients  of  the  LPT  require  larger  flow  mixing  devices  if  the  technique  is  to  be  applied 
successfully.  Considering  the  disappointing  results  of  small-scale  riblets  they  will  not  be  applied  to 
the  Pratt  and  Whitney  Pak-B  blades. 

J.1.11  Triangular  Ramp  Vortex  Generator 

Similar  to  the  plow  design  (see  Figure  152b),  this  device  was  found  by  Barber  et  al.  [76]  to 
have  an  optimal  height  of  h  «  ^  with  an  ideal  placement  ahead  of  the  original  separation  line.  The 
optimal  length  to  half-width  to  height  ratio  for  the  ramp  was  found  to  be  6:5: 1 .  Flow  spilling  over 
the  ramp  produces  trailing  edge  vortices  which  energize  the  lower  momentum  fluid  near  the  surface, 
encouraging  separation  delay  or  prohibiting  it  entirely.  Barber  et  al.  [76]  found  that  this  submerged 
concept  has  one-tenth  of  the  parasitic  drag  than  larger  devices  which  extend  into  the  freestream.  The 
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same  author  also  found  that  the  effective  range  downstream  of  the  ramp  was  only  ™<!>o,  which  for 
the  turbine  case  is  insufficient  to  manage  the  entire  separation  region. 

J.1.12  Parallel  Concave  Slat  Vortex  Generators 

Concave  slats  used  by  Rao  and  Kariya  [60]  were  constructed  from  standard  window  Venetian 
blinds  and  affixed  to  a  flat  plate  surface.  Depicted  in  Figure  152c,  these  slats  were  tested  both  in 
parallel  and  V-shaped  configurations.  The  V-shaped  test  used  the  slats  connected  at  their  leading 
edges  facing  into  the  flow,  with  the  trailing  edges  spreading  apart  at  a  30°  angle .  The  V  configuration 
however  produced  a  low-eneigy  wake  behind  the  attached  devices,  giving  the  parallel  configuration 
better  performance  and  a  higher  CPtmax  recovery.  The  design  and  implementation  of  this  device  over 
the  span  of  a  curved  turbine  airfoil  would  prove  difficult  and  perhaps  retard  performance  at  higher 
Reynolds  numbers.  The  complexity  of  installing  this  device  with  the  concerns  of  keeping  the  slats 
clean  in  an  operable  engine  eliminate  this  concept. 

J.1.13  Wheeler  Vortex  Generators 

Designed  as  a  low  speed  vortex  generator,  this  device  has  an  overall  height  of  only  h  =  0.10<5o. 
Wheeler  vortex  generators  are  placed  in  an  overlapping  fashion  on  the  surface  (see  Figure  152d), 
forming  a  “Wheeler  Doublet.  ”  Optimum  dimensions  determined  from  subsonic  testing  with  Uoo 
=  40.23  m/s  (132  ft/s)  by  Lin  et  al.  [18]  showed  placement  should  be  26q  ahead  of  the  original 
unmodified  separation  location  with  a  length  of  1.28<50.  Wheeler  Doublets  tested  by  Lin,  et  al.  [18] 
on  a  curved  ramp  simulating  the  suction  surface  of  a  blade  showed  a  reduction  in  reattachment 
distance  of  up  to  66%  over  the  baseline  case,  with  a  delay  in  separation  of  12.7  millimeters  (0.5 
inches).  Drag  forthis  device  was  estimated  at  1/44  ofthe  0.80<5o  tall  vane  type  vortex  generator  also 
tested  by  Lin,  etal.  [18],  McCormick  [61]  tested  slightly  taller  Wheeler  Doublets  with  an  h  =  0.36<50 
for  supersonic  flow  (M^  =  1 .56- 1 .65)  with  similar  beneficial  results.  This  good  performance  along 
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with  the  small  height  and  less  drag  than  many  other  devices  makes  this  a  desirable  option  for  this 
study;  however,  implementation  of  such  small  devices  over  the  highly  curved  surface  of  a  turbine 
blade  with  a  severe  adverse  pressure  gradient  would  be  very  difficult  and  so  this  concept  was  not 
tested. 

J.1.14  Ogee  Submerged  Vortex  Generators 

The  Ogee  (a  reverse  S-shaped  curve)  is  a  submerged  vortex  generator  that  is  quite  similar 
to  the  geometry  of  the  Wheeler  doublet  described  previously.  The  difference  between  the  Ogee 
submerged  vortex  generator  and  the  Wheeler  vortex  generator  is  that  the  Ogee  is  engraved  into  the 
surface  rather  than  sticking  out  of  it  as  the  Wheeler  vortex  generator  does  (see  Figure  152e).  Also 
designed  by  Wheeler  [62],  the  Ogee  has  a  depth  equal  to  9%  of  the  axial  chord  length  and  is  doubled 
up  by  placing  the  apex  of  the  next  row  of  Ogee’s  at  the  50%  location  of  the  current  row.  Optimum 
placement  from  the  baseline  separation  location  all  the  way  to  the  end  has  maintained  attached  flow 
through  turning  angles  of  65°  on  wind  tunnel  models  of  automobiles  (Wheeler  [62]).  While  quite 
capable  of  maintaining  attached  flow,  the  drastic  modification  necessary  to  the  low-pressure  turbine 
blades  would  be  far  too  extensive  for  practical  application.  Implementation  of  this  concept  would  be 
much  more  difficult  than  the  simple,  passive  fix  this  investigation  has  found  with  the  use  of  dimples. 

J.1.15  Small-Scale  Square  Ribs 

Simple  inspection  of  Figure  152,  parts  a  and  f,  shows  the  similarity  of  the  square  rib  con¬ 
cept  to  those  of  riblets.  While  not  tested  as  thoroughly  as  V-groove  riblets,  the  same  principles  of 
nondimensional  height  (h+)  and  spacing  (s+)  have  been  applied.  Maciejewski  and  Rivir  [77]  tested 
square  ribs  to  see  the  effects  on  the  heat  transfer  coefficient  in  a  linear  turbine  cascade  and  found 
an  average  reduction  of  7%.  Their  experiments  were  conducted  for  44,000  <  Re  <  410,000  and 
freestream  turbulence  of  0.8%  <  Tu  <  12.3%;  however,  no  measurements  were  made  concerning 
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the  effect  on  separation.  Maciejewski  and  Rivir  [77]  do  provide  the  insight  though  that  for  applica¬ 
tion  of  such  a  device  over  an  airfoil,  the  average  value  of  h+  across  the  blade  must  be  less  than  5  in 
order  for  the  local  values  of  h+  to  take  on  desirable  values  between  5  and  25.  Since  this  research 
will  focus  only  on  concepts  which  will  provide  a  profile  drag  reduction  for  Pak-B  LPT  blades,  the 
lack  of  any  drag  data  for  this  device  eliminated  this  concept. 

J.1.16  Large-Scale  Vortex  Generators 

Based  upon  the  early  work  of  Taylor  [45]  the  large-scale  vortex  generator  is  a  small  NACA  64- 
812  airfoil  sticking  through  the  boundary  layer  and  into  the  ffeestream  flow.  While  inherently  having 
more  drag  than  submerged  designs,  it  is  also  much  more  effective.  The  design  consists  of  a  tapered 
airfoil  of  height  h  =  1.2<5o,  with  the  quarter  chord  point  unswept.  The  root  chord  is  optimally  3. 846o 
while  the  tip  chord  is  just  1.92<5o.  The  airfoils  are  placed  spanwise  across  the  blade  in  an  alternating 
fashion  (to  provide  counter-rotating  vortices)  with  a  spacing  of  2.4<50  <  s  <  6<5o  between  vortex 
generators,  and  ahead  of  the  baseline  separation  a  distance  of  10<5o  to  30<$o-  Adjustments  might 
have  been  necessary  to  find  the  correct  angle  of  attack  for  Pak-B  airfoils,  but  1 6°  was  suggested  as  a 
starting  point  by  Taylor  [45]  for  any  airfoil  shape.  Taylor  [45]  found  through  experiments  with  these 
large  devices  that  drag  reductions  are  obtained,  but  only  when  the  initial  conditions  were  quite  poor. 
Considering  the  magnitude  of  the  separation  that  was  found  over  the  Pak-B  suction  surface  at  Re  = 
50k  and  Tu  =  1%,  these  devices  might  have  performed  well.  For  example,  testing  by  Taylor  [45]  of 
a  NACA  16-325  airfoil  at  Rec=  1.8  x  106  with  the  large-scale  vortex  generators  (see  Figure  153) 
showed  a  33%  increase  in  maximum  lift  coefficient  (C l)  and  a  general  decrease  in  drag  except  when 
Cl  dipped  below  0.4  where  flow  was  fully  attached  anyway.  Drag  increases  for  this  benign  flow 
of  low  Cl  was  approximately  72%  from  0.0 1 8  to  0.03 1 .  Considering  the  drag  increase  for  attached 


240 


flows  at  high  Reynolds  numbers  in  the  LPT  and  the  difficulty  in  implementing  and  maintaining  the 
concept  in  a  real  engine,  this  idea  was  rejected. 


Figure  153.  Large-Scale  vortex  generator  using  NACA  airfoil. 

J.1.17  Transverse  Grooves  Placed  at  Baseline  Separation  Point 

As  seen  in  Figure  1 54,  the  transverse  groove  concept  includes  square  ribs  engraved  transversely 
to  the  flow  direction  along  the  span  of  the  blade.  The  method  by  which  this  reduces  the  separation 
is  described  by  Lin  et  al.  [18]  as  a  substitution  of  several  smaller  separation  regions  which  provide 
a  wall  slip  boundary  condition.  These  smaller  regions  of  separation  replace  the  larger  separated- 
flow  region.  Optimal  placement  of  the  grooves  is  one  boundary  layer  thickness  («50)  ahead  of  the 
baseline  separation  extending  past  the  original  separation  point  to  one  boundary  layer  thickness 
behind  it,  coinciding  with  the  region  of  the  greatest  adverse  pressure  gradient.  Experiments  were 
conducted  by  Lin  et  al.  [18]  on  a  curved  ramp  with  a  turning  of  25°  and  a  radius  of  20.3  cm  (8 
inches).  The  experiments  were  conducted  at  low  velocities  of  Uoo  =40.23  m/s  (132  ft/s)  and  showed 
these  grooves  reduced  the  reattachment  distance  by  approximately  50%.  To  achieve  this  desirable 
result,  the  width  of  the  groove  valleys  was  b  =  0.144<$o,  the  rib  height  was  h  =  0.385<50,  and  the  rib 
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thickness,  t  =  0.048<So  using  the  nomenclature  of  Figure  154.  No  other  data  concerning  skin  friction 
coefficient  was  given  for  this  experiment. 


mjTJTJijnjixLr  >h 


Figure  154.  Transverse  groove  separation  control  device. 

Upon  further  consideration  this  method  was  rejected  due  to  the  higher  loss  coefficients  which 
may  result  from  use  at  high  Reynolds  numbers.  High  Reynolds  number  flows  such  as  take-off  and 
sea-level  operations  do  not  require  separation  control  and  increasing  the  boundary  layer  thickness 
with  a  boundary  layer  trip  would  be  detrimental  to  performance. 

J.1.18  Increased  Surface  Roughness  Over  Surface  of  Pak-B  Blade 

The  intentional  roughing  of  the  airfoil  surface  (see  Figure  155)  eneigizes  the  boundary  layer 
into  a  turbulent  state.  Organized  surface  roughness  as  a  means  to  reduce  separation  losses  has 
been  suggested  by  many  authors  including  Radeztsky  et  al.  [78],  Mueller  [28],  and  Yurchenko  [30], 
Previous  experiments  have  placed  roughness  strips  at  two  different  locations:  the  leading  edge,  and 
the  point  of  suction  side  peak  curvature .  Bloch  and  Mueller  [79]  conducted  experiments  for  both  the 
leading  edge  and  suction  peak  cases.  Bloch  and  Mueller’s  results  indicated  a  prevention  of  trailing 
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edge  separation  with  roughness  applied  to  the  suction  peak  of  a  Wortmann  FX-63-137  airfoil  for 
low  Reynolds  numbers.  Minimum  drag  decreased  slightly  over  the  baseline  smooth  airfoil  while 
maximum  Cl  decreased  by  21%  for  a  range  of  Reynolds  numbers  from  100  x  103  to  200  x  103. 
Halstead  and  Solomon  [24]  recommended  covering  the  entire  airfoil  with  increased  roughness  which 
may  prevent  the  loss  in  lift  coefficient  by  tripping  the  flow  into  a  turbulent  state  early  on. 


side 

- ► 


t 

roughness  strip 

Figure  155.  Surface  roughness  strip  applied  to  Pak-B  airfoil. 

The  moderate  success  of  Bloch  and  Mueller’s  (Bloch  and  Mueller  [79])  work  may  be  increased 
by  using  a  different  roughness  height  for  the  given  conditions.  Williams  [80]  found  a  range  of 
roughness  values  that  least  adversely  affected  the  loss  coefficient  for  his  cascade.  Roughness  values 
below  or  above  this  range  resulted  in  higher  loss  coefficients.  An  example  of  how  grit  roughness 
could  have  been  used  on  the  Pak-B  LPT  airfoil  would  be  a  2  mm  wide  strip  placed  5  mm  downstream 
of  the  leading  edge.  The  mean  height  of  the  roughness  elements  would  be  selected  based  on  the 
optimal  roughness  Reynolds  number: 
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where  k  is  the  roughness  height  and  ujt  is  the  velocity  at  the  top  of  the  roughness  element.  Rough¬ 
ness  studies  by  Braslow  et  al .  [8 1  ]  and  Bragg  et  al .  [  82]  on  a  NACA  00 1 2  airfoil  (for  chord  Reynolds 
number  greater  than  3  x  1 06)  have  found  an  ideal  value  for  Rck  of  600  was  useful  for  tripping  a  lam¬ 
inar  boundary  layer.  The  associated  increase  in  drag  coefficient  was  15%  for  adding  the  roughness 
strip.  Achenbach  [83]  performed  research  with  sand-roughened  spheres  and  found  a  drag  reduction 
for  very  narrow  ranges  of  Reynolds  numbers,  but  drag  increases  (compared  to  smooth  spheres)  for 
Reynolds  numbers  (based  upon  diameter)  above  3.8  x  105.  The  intentional  tripping  of  the  bound¬ 
ary  layer  using  a  roughness  strip  was  also  rejected  because  of  the  hindered  performance  at  high 
Reynolds  numbers  typical  of  take-off  conditions. 


APPENDIX  K  -  Limitations  of  LPT  Computational  Modelling 


K.1  Research  Conducted  by  Halstead  et  al.  [23] 

Halstead  et  al.  [23]  studied  different  types  of  software  codes  that  included  comparisons  to 
experimental  values  of  displacement  thickness,  <5*,  momentum  thickness,  9,  momentum  Reynolds 
number,  Ree,  and  shape  factor,  H  =  ^ .  The  software  used  by  Halstead  et  al.  [23]  can  be  divided  into 
three  basic  categories:  steady  boundary  layer  analyses,  steady  Navier-Stokes  analyses,  and  unsteady 
boundary  layer  analysis.  All  of  the  codes  failed  to  accurately  predict  transition  location  and  length. 
Many  of  the  codes  predicted  laminar  separation  where  none  existed.  Results  from  all  the  codes  used 
by  Halstead  et  al.  [23]  gave  reasonable  answers  for  momentum  and  displacement  thickness  as  well 
as  momentum  Reynolds  number.  One  steady  two-dimensional  boundary  layer  analysis  code  used 
by  Halstead  et  al.  [23]  was  STANX,  a  variation  of  the  STAN5  program  developed  by  Crawford  and 
Kays  [84],  This  program  uses  the  compressible  form  of  the  boundary  layer  equations  while  incor¬ 
porating  a  mixing  length  turbulence  model  adjusted  for  pressure  gradients.  The  program  must  have 
a  correlation  supplied  to  estimate  transition  location,  and  cannot  determine  the  length  of  transition. 
With  this  constraint,  the  flow  is  assumed  to  transition  to  turbulent  flow  at  a  single  point.  In  addition 
to  error  in  transition  prediction,  the  STANX  code  results  showed  laminar  separation  prior  to  tran¬ 
sition  on  a  second-stage  turbine  nozzle,  even  though  no  separation  was  measured  by  Halstead  [13] 
for  the  same  environmental  conditions.  This  is  in  part  due  to  the  inability  of  the  code  to  include  the 
unsteady  calming  effect  of  the  passage  of  wakes  through  the  boundary  layer  (Halstead  et  al.  [23]). 
For  the  low  Reynolds  number  cases,  the  STANX  code  overpredicted  the  profile  loss  coefficient  by 
up  to  56%  (Halstead  et  al.  [23]). 

Halstead  et  al.  [23]  also  examined  a  code  similar  to  STANX,  the  two-dimensional  KEP  program 
for  study  of  a  low  pressure  turbine  blade.  KEP  is  a  modified  version  of  STANX  which  includes  a  k- 
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e  turbulence  model  for  low  Reynolds  numbers.  This  program  uses  internal  correlations  to  estimate 
transition  onset,  eliminating  the  need  for  user  inputs.  KEP  estimated  transition  locations  far  too  aft 
along  the  blade  suction  surface  compared  to  the  measurements.  Predictions  of  loss  coefficient  by 
the  KEP  program  were  also  overpredicted  by  6.2%  to  59. 1%  depending  upon  Reynolds  number  and 
turbulence  intensity. 
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